WO2023024639A1 - 超声空化参数调整方法和超声空化装置 - Google Patents
超声空化参数调整方法和超声空化装置 Download PDFInfo
- Publication number
- WO2023024639A1 WO2023024639A1 PCT/CN2022/097141 CN2022097141W WO2023024639A1 WO 2023024639 A1 WO2023024639 A1 WO 2023024639A1 CN 2022097141 W CN2022097141 W CN 2022097141W WO 2023024639 A1 WO2023024639 A1 WO 2023024639A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- cavitation
- data
- intensity distribution
- intensity
- ultrasonic
- Prior art date
- Legal status (The legal status 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 status listed.)
- Ceased
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
- B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
- B01J19/08—Processes employing the direct application of electric or wave energy, or particle radiation; Apparatus therefor
- B01J19/10—Processes employing the direct application of electric or wave energy, or particle radiation; Apparatus therefor employing sonic or ultrasonic vibrations
-
- 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 agents, 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
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
- B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
- B01J19/0006—Controlling or regulating processes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods
- A61B17/22—Implements for squeezing-off ulcers or the like on inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; for invasive removal or destruction of calculus using mechanical vibrations; for removing obstructions in blood vessels, not otherwise provided for
- A61B17/22004—Implements for squeezing-off ulcers or the like on inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; for invasive removal or destruction of calculus using mechanical vibrations; for removing obstructions in blood vessels, not otherwise provided for using mechanical vibrations, e.g. ultrasonic shock waves
- A61B2017/22005—Effects, e.g. on tissue
- A61B2017/22007—Cavitation or pseudocavitation, i.e. creation of gas bubbles generating a secondary shock wave when collapsing
- A61B2017/22008—Cavitation or pseudocavitation, i.e. creation of gas bubbles generating a secondary shock wave when collapsing used or promoted
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N2007/0039—Ultrasound therapy using microbubbles
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
- B01J10/00—Chemical processes in general for reacting liquid with gaseous media other than in the presence of solid particles, or apparatus specially adapted therefor
- B01J10/002—Chemical processes in general for reacting liquid with gaseous media other than in the presence of solid particles, or apparatus specially adapted therefor carried out in foam, aerosol or bubbles
Definitions
- the invention relates to the technical field of ultrasonic detection and ultrasonic imaging, in particular to an ultrasonic cavitation parameter adjustment method and an ultrasonic cavitation device.
- Ultrasonic cavitation technology is widely used in medical or chemical fields. By applying ultrasonic energy, the tiny bubbles in the liquid vibrate and grow until a certain energy threshold is reached, and the bubbles collapse and close sharply to promote liquid emulsification, gel liquefaction, solid Dispersion, and the effect of enhancing drug absorption capacity.
- the technical solution provided by the prior art is to assist the operator to adjust artificially by showing the strongest MI value (Mechanical Index, mechanical index) under the current sound field, but the above solutions can only It serves as a reference for the control of steady-state cavitation and transient cavitation. Due to the large number of parameters affecting cavitation and the uneven energy distribution of the ultrasonic sound field, operators often cannot manually adjust the ultrasonic cavitation parameters to Best state.
- One of the objectives of the present invention is to provide a method for adjusting ultrasonic cavitation parameters to solve the technical problems in the prior art that cavitation parameter adjustment effects are poor and optimal parameters cannot be automatically adapted according to operator needs.
- One of the objectives of the present invention is to provide an ultrasonic cavitation device.
- an embodiment of the present invention provides a method for adjusting ultrasonic cavitation parameters, including: receiving a detection frame selection signal, analyzing the detection frame position, and obtaining the marker position; retrieving the preset intensity distribution map Collecting the cavitation intensity data at the marked position, screening the cavitation intensity distribution map with the largest cavitation intensity data, to obtain the optimal intensity distribution map; extracting the cavitation output parameters corresponding to the optimal intensity distribution map, to obtain the optimal Output parameters and output.
- the method further includes: outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters; Two sets of radio frequency data, at least two cavitation intensity distribution diagrams are obtained according to the analysis of the radio frequency data, and an intensity distribution atlas is established.
- the method specifically includes: calculating and obtaining corresponding cavitation intensity data according to the radio frequency data; respectively performing mapping on the cavitation intensity data to obtain false color data, and
- the cavitation intensity distribution map was generated from the pseudo-color data described above.
- the method specifically includes: receiving echo data corresponding to the cavitation trigger signal; performing processing on the echo data to obtain corresponding radio frequency data; extracting radio frequency data from the radio frequency data
- the sub-harmonic signal of the two-dimensional scanning plane is calculated to obtain the total pressure data of at least the first test point; according to the total pressure data of the first test point and the liquid pressure data of the first test point, the calculation is obtained. Cavitation intensity data for the first test point.
- the cavitation intensity data is the difference between the total pressure data and the liquid pressure data.
- the method specifically includes: traversing all the cavitation intensity data to obtain maximum intensity data and minimum intensity data; According to the grayscale data and the RGB mapping curve, the pseudo-color data corresponding to the grayscale data is obtained by mapping, and according to the pseudo-color data Generate a cavitation intensity distribution map; wherein, the grayscale data satisfies: x is the cavitation intensity data, g is the pseudo-color data corresponding to the cavitation intensity data, x max is the maximum intensity data, and x min is the minimum intensity data.
- the method specifically includes: retrieving and calculating the average value of the grayscale data of each cavitation intensity distribution map at the marked position in the intensity distribution atlas; A cavitation intensity distribution diagram of the data mean value to obtain an optimal intensity distribution diagram; wherein, the grayscale data mean value is used to characterize the cavitation intensity at the position of the mark.
- the method further includes: receiving an ultrasonic detection image; and superimposing the optimal intensity distribution map translucently on the ultrasonic detection image, generating a cavitation display image and outputting .
- the method specifically includes: selecting a region of interest in the ultrasonic detection image; setting the optimal intensity distribution map to have a first weight, and setting the region of interest With a second weight, the optimal intensity distribution map and the region of interest are weighted and mixed to generate and output a cavitation display image; wherein the value range of the first weight and the second weight is 0 to 1 , and the sum of the first weight and the second weight is 1.
- an embodiment of the present invention provides an ultrasonic cavitation device, including a main control module, an ultrasonic probe, an image processing module, a storage module, and a parameter setting module.
- the ultrasonic cavitation device performs any of the above tasks.
- a method for adjusting ultrasonic cavitation parameters described in the technical proposal performs parameter adjustment.
- the ultrasonic cavitation parameter adjustment method selects the cavitation intensity distribution map corresponding to the largest cavitation intensity data by analyzing the cavitation intensity data at the marked position in the preset intensity distribution map, and Adjust the output parameters to the optimal output parameters corresponding to the cavitation intensity distribution map, so that regardless of the position of the selected frame, the most suitable preset cavitation output parameters can be retrieved, avoiding manual adjustments by the operator, and achieving improved efficiency , Realize overall automation and ensure the best technical effect of output cavitation.
- Fig. 1 is a structural principle diagram of an ultrasonic cavitation device in an embodiment of the present invention
- Fig. 2 is a schematic diagram of the steps of the ultrasonic cavitation parameter adjustment method in an embodiment of the present invention
- Fig. 3 is a schematic diagram of the steps of the ultrasonic cavitation parameter adjustment method in another embodiment of the present invention.
- Fig. 4 is a schematic diagram of the steps of the cavitation intensity distribution map preset process of the first embodiment of the ultrasonic cavitation parameter adjustment method in another embodiment of the present invention
- Fig. 5 is a schematic diagram of the steps of step 221 of the second embodiment of the ultrasonic cavitation parameter adjustment method in another embodiment of the present invention.
- FIG. 6 is a schematic diagram of steps in step 222 of the third embodiment of the method for adjusting ultrasonic cavitation parameters in another embodiment of the present invention.
- Fig. 7 is a schematic diagram of the steps of the ultrasonic cavitation parameter adjustment method in another embodiment of the present invention.
- Fig. 8 is a schematic diagram of the steps of the first embodiment of the ultrasonic cavitation parameter adjustment method in another embodiment of the present invention.
- Fig. 9 is a schematic diagram of the steps of the second embodiment of the method for adjusting ultrasonic cavitation parameters in another embodiment of the present invention.
- Ultrasonic cavitation is used in medical and chemical fields. By applying ultrasonic energy, the cavitation nucleus vibrates, grows and continuously gathers energy under the action of the sound field, and finally collapses and closes sharply, thereby improving the drug absorption capacity of the target lesion to assist treatment, or Accelerate the mixing between different media to accelerate physical or chemical reactions.
- cavitation nuclei are defined as micron-sized or smaller air bubbles existing in living organisms, and specifically refer to tiny air bubbles existing in liquids in one embodiment of the present invention.
- the ultrasonic cavitation device applies ultrasonic energy through the ultrasonic probe.
- the sound field formed by the ultrasonic probe emitting ultrasonic waves is often of uneven intensity, resulting in the sound field intensity in different areas within the scanning range of the ultrasonic probe.
- the above MI value Mechanical Index, mechanical index
- the cavitation effect is also different. In this way, if the sound field is always established with the same parameters, it will be impossible to concentrate the sound field on the target to be treated, and the auxiliary effect on treatment and response will be greatly weakened.
- an embodiment of the present invention provides an ultrasonic cavitation device to solve the problems in the prior art and improve the auxiliary effect of ultrasonic cavitation.
- the ultrasonic cavitation device includes a main control module 11, an ultrasonic probe 12, an image processing module 13, a storage module 14, and a parameter setting module 15. Adjustments to enhance the effects of cavitation.
- the main control module 11 is connected with the peripheral display device and the operating device, receives the detection frame selected by the operator on the display device through the operating device, receives the mark position analyzed by the image processing module 13 and executes it in the storage module 14. Retrieval, put the retrieved qualified images into the image processing module 13 for processing, extract the corresponding parameters and output them to the parameter setting module 15 to perform parameter setting; in some specific embodiments, the main control module 11 also It can be further configured to perform steps such as outputting cavitation trigger signals, receiving echo data, receiving radio frequency data, forwarding false color data, receiving and outputting cavitation display images, etc. This point can be adjusted according to the needs of those skilled in the art, and the present invention does not This is limited.
- the display device can specifically refer to the display screen of the host computer or other computers, and the operating device can specifically refer to devices such as the host computer, keyboard and mouse set, etc., and the display device is at least used to display the .
- the image received and output by the ultrasonic cavitation device, and the operating device is at least used for inputting a detection frame selection signal.
- the ultrasonic probe 12 is connected to the main control module 11 and the parameter setting module 15. In one embodiment, it is configured to integrate the output function of cavitation, the feedback function of cavitation and the function of ultrasonic imaging. Of course, in other embodiments, the above functions can be separately Set in different devices, at this time, at least the part of the ultrasonic probe 12 that is used to output cavitation is configured to connect to the parameter setting module 15, and the part of the ultrasonic probe 12 that is used for feedback of cavitation and for ultrasonic imaging Part is configured to be connected to the main control module 11; further, the ultrasonic probe 12 is at least configured to accept the cavitation output parameter setting of the parameter setting module 15, output the cavitation effect, and receive the cavitation effect feedback data, in a specific In the embodiment, it can also be configured to output and receive ultrasound signals for imaging, which can be adjusted according to the needs of those skilled in the art.
- the image processing module 13 is connected to the main control module 11, configured to analyze and output the position of the detection frame in the image (or display screen), analyze and output the cavitation output parameters corresponding to the filtered cavitation intensity distribution map; in one implementation In the mode, it can also be configured to perform steps such as calculating the cavitation intensity distribution map according to the feedback data, establishing and outputting the intensity distribution atlas, generating and outputting the cavitation display image; of course, the image processing module 13 and the main control module 11 can be divided into The body settings can also be integrated, and the functions corresponding to the two modules can also be configured to be undertaken by the other one.
- the storage module 14 is connected to the main control module 11, configured to accept the retrieval of the main control module 11, store the intensity distribution atlas, and store the data and process images involved in the method of adjusting ultrasonic cavitation parameters; the parameter setting module 15 is connected to the main control module 15 respectively.
- the control module 11 and the ultrasonic probe 12 are configured to receive and output cavitation output parameters.
- the ultrasonic cavitation device can also specifically include an echo processing module 16, a bandpass filter 17, and a data processing module 18 connected in sequence, and the echo processing module 16 and the data processing module 18 are connected with the main control module respectively. 11 connection; wherein, the echo processing module 16 is configured to perform processing on the echo data forwarded by the main control module 11 to obtain radio frequency data, and the bandpass filter 17 is configured to extract and output subharmonic signals in the radio frequency data, the data The processing module 18 is configured to calculate cavitation intensity data according to the sub-harmonic signal and/or radio frequency data, perform mapping on the cavitation intensity data to obtain and output pseudo-color data; of course, the above-mentioned modules are not necessary technical features for realizing the expected technical effect of the present invention, Modules and functional configurations can be adjusted as needed. Some of the above-mentioned data and images are process data content, and some are final output data content. The output data content may also contain other process data content before. Of course, those skilled in the art can also use Process data content
- an embodiment of the present invention provides a method for adjusting ultrasonic cavitation parameters, as shown in Figure 2, including:
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the detection frame is an area selected by the operator in the display screen, which can be a frame with a certain area, or a point or point set, etc.
- the display screen is defined as an ultrasonic detection image in one embodiment, and the operation The operator selects an area on the ultrasonic detection image as a selected frame as required, and generates a corresponding selected signal.
- the process of analyzing the position of the detection frame can refer to the two diagonal points of the detection frame or the position of the geometric center.
- a plane Cartesian coordinate system can be established on the display screen and the position can be expressed in the form of coordinates.
- the reference point can also be set in the display screen to match the position of the geometric center. The relative position of the reference point determines the position of the detection frame, and the embodiments formed according to the usual means of those skilled in the art are within the protection scope of the present invention.
- the marked position generated by executing the above process is stored, and the main control module 11 uses the marked position as a reference to retrieve the cavitation intensity data at the marked position of each cavitation intensity distribution map in the intensity distribution atlas.
- the marked position is displayed as a box or point set with an area, scan the area area in the cavitation intensity distribution map to obtain the average or weighted cavitation intensity data and compare them; if the marked position is displayed as For a single point, locate the corresponding position in the cavitation intensity distribution map, and extract the cavitation intensity data for comparison.
- the cavitation intensity distribution diagram in this embodiment is defined as the distribution of the cavitation effect output by the ultrasonic probe 12 in the display screen under a certain parameter combination (ie, cavitation output parameters).
- the cavitation intensity near the center of the ultrasonic probe 12 is stronger than the cavitation intensity far from the center, but under different parameter settings, the regional distribution of stronger cavitation intensity and weaker cavitation intensity There is a difference with the size of the area.
- the cavitation intensity distribution map generated under different setting parameter states is preset in the device for query and call.
- the intensity distribution atlas is defined in this embodiment as a collection of multiple cavitation intensity distribution atlases under different parameter setting states, and the specific form can be adjusted according to the needs of those skilled in the art.
- the cavitation intensity data is the data at different regions in the cavitation intensity distribution diagram, and is used to represent the cavitation intensity in this region in the cavitation intensity distribution diagram;
- the cavitation output parameter is the parameter setting condition for generating each cavitation intensity distribution map.
- the cavitation output parameters correspond to the cavitation intensity distribution map, and the cavitation intensity data correspond to a certain area in the cavitation intensity distribution map.
- the above-mentioned intensity distribution atlas can be obtained by pre-computing with reference to big data, or it can be obtained by using the ultrasonic cavitation device provided by the present invention to perform detection and calculation.
- another embodiment of the present invention provides A method for adjusting ultrasonic cavitation parameters is proposed, as shown in Figure 3, which specifically includes:
- Step 21 outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters
- Step 22 analyzing and obtaining at least two sets of radio frequency data corresponding to at least two sets of cavitation trigger signals, obtaining at least two cavitation intensity distribution maps according to the analysis of the radio frequency data, and establishing an intensity distribution atlas;
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the cavitation trigger signal contains at least the data affecting the effect of cavitation in the cavitation output parameters, such as transmission frequency, transmission period, sound power and PRF (Pulse Repetition Frequency, pulse repetition frequency) and other information, two sets of cavitation output parameters and corresponding Among the cavitation trigger signals, there is a difference in at least one of the above parameters, so that at least two cavitation intensity distribution maps with different intensity distributions are obtained.
- PRF Pulse Repetition Frequency, pulse repetition frequency
- the feedback form of cavitation is mainly radio frequency signals in this embodiment, so the cavitation intensity distribution diagram is analyzed and calculated according to the radio frequency data corresponding to the cavitation trigger signal in this embodiment, and the feedback form has other embodiments
- the above steps 21 to 22 can of course be adjusted accordingly, as long as it is enough to generate two different cavitation intensity distribution maps.
- the order of steps 21 and 22 and steps 31 to 33 can be adjusted, as long as the intensity distribution atlas is preset before step 32 is performed, the expected technical effect of the present invention can be achieved.
- step 21 and step 22 can be configured to be performed in an environment different from step 31 to step 32, and the environment is One embodiment is defined as a pure water environment injected with contrast medium and free of air bubbles, which can be set up in a sink.
- the detection surface of the ultrasonic probe 12 is set parallel to the water surface of the water tank and immersed in the water.
- the cavitation intensity distribution diagram mentioned above can be a simple grayscale image, and the distribution of cavitation intensity is represented by the size of the grayscale data, but in this embodiment, in order to improve the display effect of the cavitation intensity distribution diagram, further
- the first embodiment shown in Figure 4 is provided, specifically including:
- Step 21 outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters
- Step 221 analyzing and obtaining at least two sets of radio frequency data corresponding to at least two sets of cavitation trigger signals, and calculating corresponding cavitation intensity data;
- Step 222 respectively perform mapping on the cavitation intensity data to obtain pseudo-color data, and generate a cavitation intensity distribution map based on the pseudo-color data, and establish an intensity distribution atlas;
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the mapping process of pseudo-color data is to use the cavitation intensity data as the code value, correspondingly look up the R, G, and B intensity values corresponding to the cavitation intensity data in the CLUT (Color Look-Up Table), and then execute map. Since the mapped pseudo-color data still contains cavitation intensity data representing cavitation intensity, and the human eye is more sensitive to color than grayscale, the corresponding generated cavitation intensity distribution map can more significantly represent the cavitation intensity. Distribution.
- the obtained cavitation intensity distribution diagram or intensity distribution atlas can not only be stored in the storage module 14, but also be output to a display device for display.
- the cavitation trigger signal is output to the ultrasonic probe 12.
- the ultrasonic probe 12 performs cavitation correspondingly and receives the echo signal representing the cavitation effect.
- the collection of the echo signal is performed according to the time course, and the establishment of the cavitation intensity distribution map
- the relative positional relationship between depth and lateral direction is required.
- a second implementation of the above-mentioned process is further provided. For example, as shown in Figure 5, it specifically includes:
- Step 21 outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters
- Step 2211 receiving echo data corresponding to the cavitation trigger signal
- Step 2212 perform processing on the echo data to obtain corresponding radio frequency data
- Step 2213 extract the sub-harmonic signal in the radio frequency data, and calculate the total pressure data of at least the first test point in the two-dimensional scanning plane;
- Step 2214 according to the total pressure data of the first test point and the liquid pressure data of the first test point, calculate the cavitation intensity data of the first test point;
- Step 222 respectively perform mapping on the cavitation intensity data to obtain pseudo-color data, and generate a cavitation intensity distribution map based on the pseudo-color data, and establish an intensity distribution atlas;
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the process of processing the echo data to obtain the radio frequency data may have a specific configuration manner.
- the echo data may be used to perform beamforming, and the formed receiving beam may be used as the radio frequency data.
- the radio frequency data is extracted by the band-pass filter 17 to obtain a sub-harmonic signal.
- the sub-harmonic signal has a strong correlation with the pressure in the water body. Therefore, the present embodiment utilizes the amplitude variation of the sub-harmonic signal, The total pressure data of at least the first test point in the water body can be calculated.
- the sub-harmonic signal is distinguished from the first harmonic component (the first harmonic component) or the second harmonic component (the second harmonic component), and in this embodiment specifically refers to the sub-harmonic component (subharmonic component, or subharmonic components), which have a better correlation with pressure, and the estimated total pressure data is more accurate.
- the cavitation intensity data is set as the difference between the total pressure data and the liquid pressure data, that is:
- P final is the equivalent cavitation intensity data
- P ultrasound is the sound field pressure emitted by the ultrasonic probe 12
- P bubble is the pressure generated by cavitation nucleus (micro bubble) resonance, that is, the actual cavitation intensity data
- P cal is the all The total pressure data
- P water is the liquid pressure data.
- the principle of the above derivation formula is that the actual cavitation intensity data P bubble has a strong correlation with the parameters affecting the sound field pressure P ultrasound , especially with the sparse pressure peak (PRP, Peak Rarefactional Pressure), pulse repetition frequency (PRF, Pulse Repetition Frequency), and pulse duration (PD, Pulse Duration), so it can be considered that there is a corresponding relationship between the cavitation intensity data P bubble and the sound field pressure P ultrasound , and the sound field pressure P ultrasound can be converted into the cavitation intensity data P
- the parameters of the bubble can also define an equivalent cavitation intensity data P final to represent the current cavitation intensity.
- the first test point will be simultaneously affected by the liquid pressure P water , the actual cavitation intensity P bubble and the sound field pressure P ultrasound during operation, that is, it will be affected by the liquid pressure P water and the equivalent cavitation intensity
- the combined effect of the data P final constitutes the total pressure P cal detected by the ultrasonic probe 12 .
- the above relational formula can be used to simply calculate the cavitation intensity data that can represent the cavitation intensity, that is, the equivalent cavitation intensity data.
- step 2213 and step 2214 can be split into two stages.
- the ultrasonic probe 12 No cavitation is applied to the water body, and only the liquid pressure P water of the first test point is measured.
- the ultrasonic probe 12 applies cavitation to the water body, and the total pressure data P of the first test point is measured cal , so that the data measured in the two stages are used to make a difference, and the equivalent cavitation intensity data representing the cavitation effect is obtained.
- the present invention simply calls the cavitation intensity data for comparison, and then determines the cavitation output parameters. In fact, using the original (equivalent) cavitation intensity data is enough to achieve the purpose of the invention. Further, in order to improve the display quality of the cavitation intensity distribution diagram, the third embodiment of this embodiment provides a technical solution for generating the cavitation intensity distribution diagram by mapping the pseudo-color data from the grayscale data, as shown in FIG. 6 , Specifically include:
- Step 21 outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters
- Step 221 analyzing and obtaining at least two sets of radio frequency data corresponding to at least two sets of cavitation trigger signals, and calculating corresponding cavitation intensity data;
- Step 2221 traversing all cavitation intensity data to obtain maximum intensity data and minimum intensity data
- Step 2222 according to the maximum intensity data, the minimum intensity data and the cavitation intensity data, calculate the grayscale data corresponding to the cavitation intensity data;
- Step 2223 according to the grayscale data and the RGB mapping curve, map to obtain the pseudo-color data corresponding to the grayscale data, and generate a cavitation intensity distribution map according to the pseudo-color data;
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the grayscale data satisfies: x is the cavitation intensity data, g is the pseudo-color data corresponding to the cavitation intensity data, x max is the maximum intensity data, and x min is the minimum intensity data.
- the cavitation intensity data with a large data range can be first mapped to the grayscale data in the range of 0 to 255, and then the grayscale data can be mapped to pseudo-color data by means of table lookup or curve mapping, and finally according to the current space Generate a cavitation intensity distribution map for all the pseudo-color data corresponding to the output parameters.
- the cavitation intensity data will not be lost during the process of mapping and generating the cavitation intensity distribution map, or set an independent step before step 32 to convert the false color data
- the grayscale data is re-converted into cavitation intensity data, or the retrieved cavitation intensity data in step 32 is changed to retrieved grayscale data or retrieved pseudo-color data, all of which can achieve the desired technical effect.
- the cavitation intensity distribution map generated by the above steps is displayed in color, the area with large cavitation intensity data value (strong cavitation effect) is displayed in dark red or red, and the area with small cavitation intensity data value (weak cavitation effect) It is shown as purple or blue-purple, and the areas in between are shown as yellow, green, etc.
- FIG. 7 For the technical solution of using grayscale data instead of cavitation intensity data to perform retrieval and screening, another embodiment of the present invention provides an ultrasonic cavitation parameter adjustment method as shown in Figure 7, which specifically includes:
- Step 21 outputting at least two sets of cavitation trigger signals with at least two sets of preset cavitation output parameters
- Step 221 analyzing and obtaining at least two sets of radio frequency data corresponding to at least two sets of cavitation trigger signals, and calculating corresponding cavitation intensity data;
- Step 2221 traversing all cavitation intensity data to obtain maximum intensity data and minimum intensity data
- Step 2222 according to the maximum intensity data, the minimum intensity data and the cavitation intensity data, calculate the grayscale data corresponding to the cavitation intensity data;
- Step 2223 according to the grayscale data and the RGB mapping curve, map to obtain the pseudo-color data corresponding to the grayscale data, and generate a cavitation intensity distribution map according to the pseudo-color data;
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 321 retrieve and calculate the mean value of the gray data at the marked position of each cavitation intensity distribution map in the intensity distribution map set;
- Step 322 screening the cavitation intensity distribution diagram with the largest gray data mean value to obtain the optimal intensity distribution diagram
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution graph, obtaining the optimal output parameters and outputting them.
- the mean value of the gray data is used to represent the cavitation intensity data at the marked position.
- the present invention provides another embodiment, which specifically includes the first embodiment as shown in Figure 8 and the second embodiment as shown in Figure 9 example.
- Step 30 receiving an ultrasonic detection image
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieve the cavitation intensity data at the marked position in the preset intensity distribution atlas set, and filter the cavitation intensity distribution diagram with the largest cavitation intensity data to obtain the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution diagram, obtaining the optimal output parameters and outputting them;
- Step 34 translucently superimposing the optimal intensity distribution diagram on the ultrasonic detection image, generating and outputting a cavitation display image.
- Ultrasonic detection images are usually black and white images, and the optimal intensity distribution map that is processed into translucency and displayed as color is superimposed on the ultrasonic detection image, which can clearly show the current detection target and space.
- the relationship between the distribution of cavitation intensity so that the operator understands the current situation of detection (or imaging) and cavitation application.
- other methods can be used to generate a cavitation display image that can integrate the characteristics of both the optimal intensity distribution map and the ultrasonic detection image, and the expected technical effect can also be achieved.
- step 30 and step 34 is also not limited, as long as step 30 is executed before step 34 and step 34 is executed after step 32, the expected technical effect can be achieved.
- Step 30 receiving an ultrasonic detection image
- Step 31 receiving the detection frame selection signal, analyzing the detection frame position, and obtaining the marker position;
- Step 32 retrieving the cavitation intensity data at the marked position in the preset intensity distribution atlas, screening the cavitation intensity distribution diagram with the largest cavitation intensity data, and obtaining the optimal intensity distribution diagram;
- Step 33 extracting the cavitation output parameters corresponding to the optimal intensity distribution diagram, obtaining the optimal output parameters and outputting them;
- Step 341 selecting a region of interest in the ultrasonic testing image
- Step 342 setting the optimal intensity distribution map to have a first weight, setting the region of interest to have a second weight, performing weighted mixing of the optimal intensity distribution map and the region of interest, generating a cavitation display image and outputting it.
- the first weight and the second weight range from 0 to 1, and the sum of the first weight and the second weight is 1. In this way, by adjusting the values of the first weight and the second weight, the transparency of the optimal intensity distribution map can be adjusted, and the display effect of the cavitation display image can be changed.
- the region of interest (ROI, Region of Interest) described here is applied in the field of image processing and machine vision, and is defined as a frame, circle, ellipse, etc. from the processed image (ultrasonic detection image in this embodiment) Outline the area that needs to be processed, often used in Halcon, OpenCV, Matlab and other software.
- the ultrasonic detection image may be globally selected as the region of interest, or a detection frame or other regions input by the operator may be selected as the region of interest.
- the ultrasonic detection image, the cavitation intensity distribution map, and the cavitation display image have the same size, preferably the same size as the display screen.
- the order of the steps can be adjusted according to the needs of those skilled in the art without affecting the realization of the technical effect. It should be noted that the multiple ultrasonic cavitation parameter adjustment methods provided by the present invention cannot be viewed in isolation, and the steps of each embodiment or example can of course be combined and/or replaced, and the new embodiment thus generated should be included in the present invention within the scope of protection.
- the ultrasonic cavitation parameter adjustment method selects the cavitation intensity distribution map corresponding to the largest cavitation intensity data by analyzing the cavitation intensity data at the marked position in the preset intensity distribution map set, and adjusts the output parameters as The optimal output parameters corresponding to the cavitation intensity distribution map, so that regardless of the position of the selected frame, the most suitable preset cavitation output parameters can be retrieved, avoiding manual adjustments by the operator, achieving improved efficiency and global automation , To ensure the best technical effect of output cavitation.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Organic Chemistry (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Toxicology (AREA)
- Biophysics (AREA)
- Hematology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
Claims (10)
- 一种超声空化参数调整方法,其特征在于,包括:接收检测框选定信号,分析所述检测框位置,得到标记位置;检索预设的强度分布图集中所述标记位置的空化强度数据,筛选具有最大空化强度数据的空化强度分布图,得到最优强度分布图;提取所述最优强度分布图对应的空化输出参数,得到最优输出参数并输出。
- 根据权利要求1所述的超声空化参数调整方法,其特征在于,所述方法还包括:以至少两组预设空化输出参数输出至少两组空化触发信号;分析得到与所述至少两组空化触发信号对应的至少两组射频数据,根据所述射频数据分析得到至少两张空化强度分布图,建立强度分布图集。
- 根据权利要求2所述的超声空化参数调整方法,其特征在于,所述方法具体包括:根据所述射频数据,计算得到对应的空化强度数据;分别对所述空化强度数据执行映射,得到伪彩数据,并根据所述伪彩数据生成空化强度分布图。
- 根据权利要求3所述的超声空化参数调整方法,其特征在于,所述方法具体包括:接收与所述空化触发信号对应的回波数据;对所述回波数据执行处理,对应得到射频数据;提取所述射频数据中的次谐波信号,计算得到二维扫描平面内至少第一测试点的总压强数据;根据所述第一测试点的总压强数据和所述第一测试点的液体压强数据,计算得到所述第一测试点的空化强度数据。
- 根据权利要求4所述的超声空化参数调整方法,其特征在于,所述空化强度数据为所述总压强数据和所述液体压强数据的差。
- 根据权利要求6所述的超声空化参数调整方法,其特征在于,所述方法具体包括:检索并计算所述强度分布图集中,每张空化强度分布图在所述标记位置处的灰度数据均值;筛选具有最大灰度数据均值的空化强度分布图,得到最优强度分布图;其中,所述灰度数据均值用以表征所述标记位置处的所述空化强度。
- 根据权利要求1所述的超声空化参数调整方法,其特征在于,所述方法还包括:接收超声检测图像;以及,将所述最优强度分布图半透明叠加于所述超声检测图像之上,生成空化显示图像并输出。
- 根据权利要求8所述的超声空化参数调整方法,其特征在于,所述方法具体包括:在所述超声检测图像中选定感兴趣区域;设定所述最优强度分布图具有第一权重,设定所述感兴趣区域具有第二权重,将所述最优强度分布图和所述感兴趣区域进行加权混合,生成空化显示图像并输出;其中所述第一权重和所述第二权重取值范围为0至1,且所述第一权重和所述第二权重之和为1。
- 一种超声空化装置,其特征在于,包括主控模块、超声探头、图像处理模块、存储模块和参数设定模块,所述超声空化装置执行权利要求1-9任一项所述的超声空化参数调整方法进行参数调整。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP22859986.6A EP4393582A4 (en) | 2021-08-23 | 2022-06-06 | ULTRASONIC CAVITATION PARAMETER ADJUSTMENT METHOD AND ULTRASONIC CAVITATION APPARATUS |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110970195.2 | 2021-08-23 | ||
| CN202110970195.2A CN113663622B (zh) | 2021-08-23 | 2021-08-23 | 超声空化参数调整方法和超声空化装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2023024639A1 true WO2023024639A1 (zh) | 2023-03-02 |
Family
ID=78545189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2022/097141 Ceased WO2023024639A1 (zh) | 2021-08-23 | 2022-06-06 | 超声空化参数调整方法和超声空化装置 |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP4393582A4 (zh) |
| CN (1) | CN113663622B (zh) |
| WO (1) | WO2023024639A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117038530A (zh) * | 2023-10-07 | 2023-11-10 | 东莞市楷德精密机械有限公司 | 一种半导体清洗设备的智能控制方法及系统 |
| CN119224115A (zh) * | 2024-12-02 | 2024-12-31 | 东北石油大学三亚海洋油气研究院 | 一种基于超声作用的破乳及乳化分析系统及方法 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113663622B (zh) * | 2021-08-23 | 2022-12-23 | 飞依诺科技股份有限公司 | 超声空化参数调整方法和超声空化装置 |
| CN114813037B (zh) * | 2022-04-21 | 2023-06-20 | 中国船舶科学研究中心 | 一种空化流动结构频率分布特征分析方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102249367A (zh) * | 2011-04-02 | 2011-11-23 | 河海大学常州校区 | 水处理中超声声强检测与控制方法 |
| CN103235041A (zh) * | 2013-04-26 | 2013-08-07 | 西安交通大学 | 基于超声主动空化成像的空化起始阈值分布重建方法 |
| US20150141734A1 (en) * | 2012-03-08 | 2015-05-21 | (Inserm) Institut National De La Sante Et De La Re Cherche Medicale | Method for determining optimized parameters of a device generating a plurality of ultrasound beams focused in a region of interest |
| CN105496455A (zh) * | 2015-12-10 | 2016-04-20 | 飞依诺科技(苏州)有限公司 | 一种超声空化强度的调节方法和装置 |
| CN108827876A (zh) * | 2018-04-03 | 2018-11-16 | 中国科学院声学研究所 | 一种空化云空化强度的测量方法 |
| CN111220700A (zh) * | 2019-12-09 | 2020-06-02 | 中北大学 | 超声空化泡运动矢量估计方法 |
| CN113663622A (zh) * | 2021-08-23 | 2021-11-19 | 飞依诺科技(苏州)有限公司 | 超声空化参数调整方法和超声空化装置 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7238085B2 (en) * | 2003-06-06 | 2007-07-03 | P.C.T. Systems, Inc. | Method and apparatus to process substrates with megasonic energy |
| US20100056924A1 (en) * | 2006-11-20 | 2010-03-04 | Koninklijke Philips Electronics N.V. | Control and display of ultrasonic microbubble cavitation |
| JP2015185813A (ja) * | 2014-03-26 | 2015-10-22 | 株式会社Screenホールディングス | 基板洗浄方法および基板洗浄装置 |
| CN104215581B (zh) * | 2014-07-30 | 2017-08-25 | 中国科学院声学研究所 | 一种检测超声空化强度的装置及方法 |
| CN105548006B (zh) * | 2015-11-27 | 2019-07-16 | 西安科技大学 | 一种检测超声空化强度的系统 |
| WO2018071551A1 (en) * | 2016-10-11 | 2018-04-19 | Thomas Jefferson University | Non-invasive method for pressure measurement |
| US11259781B2 (en) * | 2016-12-16 | 2022-03-01 | Koninklijke Philips N.V. | Adaptive pulsing for sonothrombolysis treatment |
| CN106950832B (zh) * | 2017-03-08 | 2020-01-31 | 杭州电子科技大学 | 一种利用空化强度反馈的超声分散控制装置 |
| KR102123072B1 (ko) * | 2019-01-29 | 2020-06-16 | 제주대학교 산학협력단 | 수동 공동 영상을 표시하는 방법 및 장치 |
| CN113040814B (zh) * | 2021-03-11 | 2023-12-22 | 飞依诺科技股份有限公司 | 超声微泡空化设备的成像处理方法及成像处理系统 |
-
2021
- 2021-08-23 CN CN202110970195.2A patent/CN113663622B/zh active Active
-
2022
- 2022-06-06 WO PCT/CN2022/097141 patent/WO2023024639A1/zh not_active Ceased
- 2022-06-06 EP EP22859986.6A patent/EP4393582A4/en not_active Withdrawn
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102249367A (zh) * | 2011-04-02 | 2011-11-23 | 河海大学常州校区 | 水处理中超声声强检测与控制方法 |
| US20150141734A1 (en) * | 2012-03-08 | 2015-05-21 | (Inserm) Institut National De La Sante Et De La Re Cherche Medicale | Method for determining optimized parameters of a device generating a plurality of ultrasound beams focused in a region of interest |
| CN103235041A (zh) * | 2013-04-26 | 2013-08-07 | 西安交通大学 | 基于超声主动空化成像的空化起始阈值分布重建方法 |
| CN105496455A (zh) * | 2015-12-10 | 2016-04-20 | 飞依诺科技(苏州)有限公司 | 一种超声空化强度的调节方法和装置 |
| CN108827876A (zh) * | 2018-04-03 | 2018-11-16 | 中国科学院声学研究所 | 一种空化云空化强度的测量方法 |
| CN111220700A (zh) * | 2019-12-09 | 2020-06-02 | 中北大学 | 超声空化泡运动矢量估计方法 |
| CN113663622A (zh) * | 2021-08-23 | 2021-11-19 | 飞依诺科技(苏州)有限公司 | 超声空化参数调整方法和超声空化装置 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117038530A (zh) * | 2023-10-07 | 2023-11-10 | 东莞市楷德精密机械有限公司 | 一种半导体清洗设备的智能控制方法及系统 |
| CN117038530B (zh) * | 2023-10-07 | 2024-01-16 | 东莞市楷德精密机械有限公司 | 一种半导体清洗设备的智能控制方法及系统 |
| CN119224115A (zh) * | 2024-12-02 | 2024-12-31 | 东北石油大学三亚海洋油气研究院 | 一种基于超声作用的破乳及乳化分析系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113663622B (zh) | 2022-12-23 |
| EP4393582A1 (en) | 2024-07-03 |
| CN113663622A (zh) | 2021-11-19 |
| EP4393582A4 (en) | 2024-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113663622B (zh) | 超声空化参数调整方法和超声空化装置 | |
| US20220354467A1 (en) | Methods and apparatus for configuring an ultrasound system with imaging parameter values | |
| US8047992B2 (en) | Brightness adjustment method and system for 3D ultrasound | |
| CN100448406C (zh) | 超声波诊断装置和超声波诊断方法 | |
| US20210369241A1 (en) | Imaging system and method with live examination completeness monitor | |
| EP4033983B1 (en) | Methods and apparatus for configuring an ultrasound system with imaging parameter values | |
| CN1644170A (zh) | 超声波诊断装置、超声波图像数据处理方法以及图像处理装置 | |
| JP4831539B2 (ja) | C面ボリューム複合イメージングのための方法及び装置 | |
| US20260114847A1 (en) | Methods and Systems for Evaluating Echo Data Contemporaneous with an Electrodiagnostic Study | |
| US11099263B2 (en) | Automatic time gain compensation (TGC) in ultrasound imaging | |
| CN106691510A (zh) | 优化超声波图像的方法和装置 | |
| JP2005074226A (ja) | 多数の強調画像を提示する方法及び装置 | |
| CN113509209B (zh) | 一种眼科超声成像方法及装置 | |
| CN111714157A (zh) | 一种多普勒超声血流自动识别方法及装置 | |
| US10716545B2 (en) | Ultrasound system for imaging and protecting ophthalmic or other sensitive tissues | |
| JP6733445B2 (ja) | 超音波診断装置、超音波画像生成方法及びプログラム | |
| CN107874782A (zh) | 一种多普勒超声自动定位取样框的方法 | |
| JP2021159696A (ja) | 超音波診断装置、およびプログラム | |
| US20180146954A1 (en) | Method of ultrasound apparatus parameters configuration and an ultrasound apparatus of using the same | |
| WO2021061957A1 (en) | Methods and apparatuses for providing feedback for positioning an ultrasound device | |
| KR101097578B1 (ko) | 3차원 컬러 도플러 영상 제공 방법 및 그를 위한 초음파 시스템 | |
| US20170000463A1 (en) | Ultrasonic diagnostic apparatus | |
| CN108230250A (zh) | 超声图像对比度优化系统及方法 | |
| CN120285438B (zh) | 一种辅助治疗关节炎的干扰电设备 | |
| CN111281423A (zh) | 一种超声图像优化方法和超声成像设备 |
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: 22859986 Country of ref document: EP Kind code of ref document: A1 |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 2022859986 Country of ref document: EP |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| ENP | Entry into the national phase |
Ref document number: 2022859986 Country of ref document: EP Effective date: 20240325 |
|
| WWW | Wipo information: withdrawn in national office |
Ref document number: 2022859986 Country of ref document: EP |