CN117637173A - Method and device for determining microwave ablation parameters and electronic equipment - Google Patents
Method and device for determining microwave ablation parameters and electronic equipment Download PDFInfo
- Publication number
- CN117637173A CN117637173A CN202311618626.4A CN202311618626A CN117637173A CN 117637173 A CN117637173 A CN 117637173A CN 202311618626 A CN202311618626 A CN 202311618626A CN 117637173 A CN117637173 A CN 117637173A
- Authority
- CN
- China
- Prior art keywords
- ablation
- initial
- simulation
- parameters
- model
- 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.)
- Pending
Links
- 238000002679 ablation Methods 0.000 title claims abstract description 393
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000004088 simulation Methods 0.000 claims description 202
- 230000000694 effects Effects 0.000 description 22
- 238000010586 diagram Methods 0.000 description 18
- 206010028980 Neoplasm Diseases 0.000 description 14
- 210000004185 liver Anatomy 0.000 description 12
- 230000006870 function Effects 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- 208000024770 Thyroid neoplasm Diseases 0.000 description 1
- 206010046798 Uterine leiomyoma Diseases 0.000 description 1
- 238000010317 ablation therapy Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000000415 inactivating effect Effects 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 208000014018 liver neoplasm Diseases 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 235000015277 pork Nutrition 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 201000002510 thyroid cancer Diseases 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Abstract
The invention provides a method and a device for determining microwave ablation parameters and electronic equipment, and size data of an area to be ablated is obtained; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time. According to the mode, the microwave ablation power and the microwave ablation time matched with the area to be ablated can be obtained by inputting the size data of the area to be ablated into the ablation prediction model, and therefore the ablation efficiency can be improved.
Description
Technical Field
The present invention relates to the field of medical technologies, and in particular, to a method and an apparatus for determining microwave ablation parameters, and an electronic device.
Background
The tumor microwave ablation therapy is connected with an ablation needle through microwave ablation equipment, and a microwave source emits microwave to generate heat so as to solidify protein to achieve the purpose of inactivating cells of the tumor, and is widely applied to tumor operations such as liver cancer, lung cancer, thyroid cancer, hysteromyoma and the like. In the related art, a doctor scans by imaging equipment such as CT (Computed Tomography) or MRI (nuclear Magnetic Resonance Imaging, magnetic resonance imaging) to obtain a tumor size before an operation, and a commercial microwave ablation device can provide an ablation effect table for the doctor, and displays the size of an ablation area of 55w,60w,70w corresponding to 5min,8min,10min,12min and 15 min. Therefore, the related art cannot perform efficient ablation because a doctor cannot determine the corresponding relationship between the size of the area to be ablated and the power and time before operation.
Disclosure of Invention
The invention aims to provide a method and a device for determining microwave ablation parameters and electronic equipment so as to improve the ablation efficiency.
The invention provides a method for determining microwave ablation parameters, which comprises the following steps: acquiring size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
Further, the ablation prediction model includes: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; the step of inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model comprises the following steps of: inputting the long-diameter data into a long-diameter submodel to obtain a first submodel; inputting the short-diameter data into a short-diameter submodel to obtain a second submodel; inputting the volume data into a volume sub-model to obtain a third sub-model; and calculating microwave ablation parameters matched with the region to be ablated based on the first sub-model, the second sub-model and the third sub-model.
Further, the ablation predictive model is determined by: obtaining a plurality of groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration; inputting the group of simulation parameters into a pre-verified simulation model aiming at each group of simulation parameters, and reading simulation cross section sizes of designated cross sections of simulation ablation areas corresponding to a plurality of first time points in simulation ablation time under simulation power; wherein the simulated ablation area is ellipsoidal; the simulated cross-sectional dimensions include: simulating a long diameter dimension and a short diameter dimension; determining the simulation volume size corresponding to the simulation ablation area according to the simulation cross section size; performing curve fitting on a plurality of simulated long diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a long diameter sub-model; performing curve fitting on a plurality of simulated short diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a short diameter sub-model; and performing curve fitting on a plurality of simulation volume sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a volume sub-model.
Further, the pre-validated simulation model is determined by: acquiring at least one set of initial parameters; wherein each set of initial parameters comprises: an initial power and an initial ablation duration; inputting the initial parameters into a preset initial simulation model aiming at each initial parameter group, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation duration under initial power; acquiring an actual ablation result of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power; and comparing the initial simulation result with the actual ablation result to obtain a comparison result so as to determine a pre-verified simulation model based on the comparison result.
Further, the initial simulation result includes: an initial cross-sectional dimension and an initial volumetric dimension; inputting the initial parameters into a preset initial simulation model aiming at each initial parameter, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation duration under initial power, wherein the steps comprise: inputting the initial parameters into a preset initial simulation model aiming at each initial parameter, and reading initial section sizes of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation time under initial power; wherein the initial ablation region is ellipsoidal; the initial cross-sectional dimensions include: an initial long diameter dimension and an initial short diameter dimension; and determining the initial volume size corresponding to the initial ablation region according to the initial cross-sectional size.
Further, the actual ablation results include: actual cross-sectional dimensions and actual volumetric dimensions; the step of obtaining the actual ablation result of the designated section of the actual ablation area of the preset sample to be ablated corresponding to each second time point under the initial power comprises the following steps: acquiring the actual section size of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power according to each group of initial parameters; wherein, actual ablation area is oval, and actual cross-section size includes: an actual long diameter dimension and an actual short diameter dimension; and determining the actual volume size corresponding to the actual ablation region according to the actual cross-sectional size.
Further, the step of comparing the initial simulation result with the actual ablation result to obtain a comparison result, so as to determine a pre-verified simulation model based on the comparison result comprises the following steps: calculating the difference value between the initial long diameter size and the actual long diameter size corresponding to each second time point to obtain a plurality of first difference value results; calculating the difference value between the initial short diameter size and the actual short diameter size corresponding to each second time point to obtain a plurality of second difference value results; comparing the difference value between the initial volume size and the actual volume size corresponding to each second time point to obtain a plurality of third difference value results; if at least one first difference result exceeds a preset first threshold range, and/or at least one second difference result exceeds a preset second threshold range, and/or at least one third difference result exceeds a preset third threshold range, receiving a model parameter adjustment instruction for the initial simulation model to adjust the initial simulation model to obtain a second simulation model; and taking the second simulation model as a new initial simulation model, and repeatedly executing the step of acquiring at least one group of initial parameters until each first difference result belongs to a first threshold range, each second difference result belongs to a second threshold range, and each third difference result belongs to a third threshold range, so as to obtain a pre-verified simulation model.
The invention provides a device for determining microwave ablation parameters, which comprises: the acquisition module is used for acquiring the size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; the output module is used for inputting the size data into a preset ablation prediction model so as to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
The invention provides an electronic device, which comprises a processor and a memory, wherein the memory stores machine executable instructions which can be executed by the processor, and the processor executes the machine executable instructions to realize the method for determining the microwave ablation parameters in any one of the above steps.
The invention provides a machine-readable storage medium storing machine-executable instructions that, when invoked and executed by a processor, cause the processor to implement a method of determining microwave ablation parameters of any of the above.
The invention provides a method and a device for determining microwave ablation parameters and electronic equipment, and the method and the device are used for acquiring size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time. According to the mode, the microwave ablation power and the microwave ablation time matched with the area to be ablated can be obtained by inputting the size data of the area to be ablated into the ablation prediction model, and therefore the ablation efficiency can be improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for determining microwave ablation parameters according to an embodiment of the present invention;
FIG. 2 is a diagram of a system for microwave ablation and simulation data acquisition provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram showing a comparison of an actual ablation section and a simulated section according to an embodiment of the present invention;
FIG. 4 (a) is a schematic diagram showing a comparison of actual ablation data and simulation data provided by an embodiment of the present invention;
FIG. 4 (b) is a schematic diagram showing a comparison of actual ablation data and simulation data provided by an embodiment of the present invention;
FIG. 5 (a) is a schematic diagram showing a comparison of short-diameter fitting functions according to an embodiment of the present invention;
FIG. 5 (b) is a schematic diagram showing a comparison of long-diameter fitting functions according to an embodiment of the present invention;
FIG. 5 (c) is a schematic diagram showing a comparison of a volume fitting function according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a final predictive model of ablation according to an embodiment of the invention;
FIG. 7 is a schematic diagram of comparison of ablation data according to an embodiment of the present invention;
fig. 8 is a schematic structural diagram of a device for determining microwave ablation parameters according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
The technical solutions of the present invention will be clearly and completely described in connection with the embodiments, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
At present, however, there are many problems in practical clinical use, and development of microwave ablation is limited, wherein one of the most important pain points is that a doctor can scan the tumor size by imaging equipment such as CT or MRI before operation, but does not know how much power and time of microwave ablation can ablate how large tumor area. At present, commercial microwave ablation equipment can give a table of ablation effects to doctors, however, due to the influence of various factors, the power can be different even at the same time, so that only a approximate range can be given for the doctors to refer, and the doctors cannot acquire ablation effect parameters of other powers and times. Therefore, in the actual use process, a doctor needs to select a relatively proper ablation power and time according to own operation experience, but the situation that the tumor is not completely ablated at one time and needs secondary ablation is often caused by insufficient setting of the ablation power or time; or because of excessive power or timing, many surrounding normal tissue cells are killed during ablation. Secondly, the doctor does not know the corresponding relation between the size of the ablated region, the power and the time before operation, so that the efficient ablation cannot be performed. The ablation effect which can be achieved after 5min of ablation is possible, and a doctor cannot know the exact ablation size, so that the time of 2-3min of ablation can be increased to ensure successful operation, and therefore, the corresponding relationship between the size of the area to be ablated and the power and time cannot be determined before operation by the doctor in the related technology, so that the efficient ablation cannot be achieved.
For the understanding of this embodiment, a method for determining microwave ablation parameters disclosed in the embodiment of the present invention is first described, as shown in fig. 1, and the method includes the following steps:
step S102, acquiring size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated.
The area to be ablated may be an area including a tumor or the like which needs to be actually ablated; in actual implementation, the actual size of the tumor can be obtained by scanning by imaging equipment such as CT or MRI before operation, an area to be ablated is determined according to the actual size of the tumor, the area to be ablated is usually an ellipsoidal area containing the tumor, long diameter data, short diameter data and volume data of the area to be ablated can be determined, and the long diameter data can be understood as the longest diameter of the ellipsoidal area to be ablated; short diameter data can be understood as the shortest diameter of an ellipsoidal region to be ablated.
Step S104, inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
In actual implementation, after the long diameter data, the short diameter data and the volume data of the area to be ablated are obtained, the data can be input into a preset ablation prediction model, proper microwave ablation power and proper microwave ablation time can be predicted through the ablation prediction model, and the problem that the area to be ablated is ablated according to the microwave ablation power and the microwave ablation time, an ablation result matched with the area to be ablated can be obtained, and secondary ablation is required or a plurality of surrounding normal tissue cells are killed in the ablation process due to improper microwave ablation power and microwave ablation time is avoided.
The method for determining the microwave ablation parameters acquires the size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time. According to the mode, the microwave ablation power and the microwave ablation time matched with the area to be ablated can be obtained by inputting the size data of the area to be ablated into the ablation prediction model, and therefore the ablation efficiency can be improved.
The embodiment of the invention also provides another method for determining microwave ablation parameters, which is realized on the basis of the method, and the ablation prediction model comprises the following steps: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; in actual implementation, mathematical models respectively corresponding to the long-diameter submodel, the short-diameter submodel, and the volumetric submodel may be expressed as follows:
short diameter submodel: l (L) Short diameter ＝0.6870+0.0044t+0.0226P-0.0001P 2 ，R 2 ＝0.9612；
Long diameter submodel: l (L) Long diameter ＝1.3859+0.003t+0.0495P-0.0002P 2 ，R 2 ＝0.9830；
Volume submodel: v= -7.3842+0.0397t+0.1962p+0.0005pt-0.0005P 2 ，R 2 ＝0.9864；
Wherein L is Short diameter Short diameter data representing the region to be ablated; l (L) Long diameter Long diameter data representing the region to be ablated; v represents volumetric data of the region to be ablated; p represents microwave ablation power matched with the region to be ablated; t represents the microwave ablation time matched to the area to be ablated.
The method comprises the following steps:
step one, acquiring size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated.
And secondly, inputting the long-diameter data into the long-diameter submodel to obtain a first submodel.
And thirdly, inputting the short-diameter data into the short-diameter submodel to obtain a second submodel.
And step four, inputting the volume data into the volume submodel to obtain a third submodel.
In actual implementation, long diameter data of the area to be ablated can be substituted into the long diameter submodel, short diameter data of the area to be ablated can be substituted into the short diameter submodel, and volume data of the area to be ablated can be substituted into the volume submodel, so that a first submodel, a second submodel and a third submodel which are respectively corresponding to each other are obtained.
The ablation predictive model is determined by:
step A, obtaining a plurality of groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration.
The simulation power and the simulation ablation duration in each set of simulation parameters can be selected and set according to actual requirements, for example, the simulation power ranges from 30W to 80W, the step is 10W, the simulation ablation duration is 20min, and then the multiple sets of simulation parameters can be: 30W/20min;40W/20min;50W/20min;60W/20min;70W/20min;80W/20min.
Step B, inputting the group of simulation parameters into a pre-verified simulation model aiming at each group of simulation parameters, and reading simulation cross section sizes of designated cross sections of simulation ablation areas corresponding to a plurality of first time points in simulation ablation time under simulation power; wherein the simulated ablation area is ellipsoidal; the simulated cross-sectional dimensions include: a simulated long diameter dimension and a simulated short diameter dimension.
The first time point can be selected according to actual requirements, for example, a time point can be selected every other minute within the simulation ablation duration to obtain a plurality of first time points; the specified cross-section is typically the largest cross-sectional area of the ellipsoidal simulated ablation region; in actual implementation, after each group of simulation parameters is obtained, each group of simulation parameters can be respectively input into a pre-verified simulation model, simulation is performed, and in the simulation process, the simulation long diameter size and the simulation short diameter size of the designated section of the simulation ablation area corresponding to each first time point are read.
The simulation model verified in advance is determined by the following steps B1 to B4:
step B1, obtaining at least one group of initial parameters; wherein each set of initial parameters comprises: an initial power and an initial ablation duration.
A diagram of a microwave ablation and simulation data acquisition system as shown in fig. 2, wherein the microwave ablation and simulation data acquisition system comprises: a microwave ablation module and a simulation module. Taking the sample as an example of isolated pork liver, fig. 2 includes: the device comprises a microwave ablation device 20, a microwave ablation needle 21, an isolated pig liver sample 23, a PC 24 and an ablation simulation software interface 25. The microwave ablation device 20 is connected with the microwave ablation needle 21, the microwave ablation needle 21 is inserted into an isolated pig liver sample 23 for ablation in the experiment, the microwave ablation device 20 is connected with the PC 24 so as to run ablation simulation software, and an ablation simulation software interface 25 is displayed. Wherein the microwave ablation module comprises a microwave ablation device 20 (for example, the microwave ablation device can comprise a 2450MHz microwave source, a water pump, upper computer software), a microwave ablation needle 21 and the like; the simulation modules include PC (Personal Computer ) machine 24 and COMSOL (COMSOL Multiphysics, a general purpose multiple physical field simulation software) ablation simulation software, and the like. In simulation, parameter settings for the ablation model typically include: the dielectric constant, electrical conductivity and thermophysical parameters of liver include thermal conductivity, density and specific heat capacity. The initial power and the initial ablation duration in each set of initial parameters can be selected and set according to actual requirements, for example, the initial power is 50W, the initial ablation duration is 20min, and the like.
Step B2, inputting the initial parameters into a preset initial simulation model aiming at each initial parameter group, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation duration under initial power;
the initial simulation result comprises: an initial cross-sectional dimension and an initial volumetric dimension; the step B2 can be specifically realized by the following steps B21 and B22:
step B21, inputting the initial parameters into a preset initial simulation model for each group of initial parameters, and reading initial section sizes of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation time under initial power; wherein the initial ablation region is ellipsoidal; the initial cross-sectional dimensions include: an initial long diameter dimension and an initial short diameter dimension.
The second time point may be selected according to actual requirements, for example, a time point may be selected every one minute within the simulated ablation duration to obtain a plurality of second time points, or a plurality of time points such as 5min,8min,10min, etc. are selected as a plurality of second time points; the above specified cross-section is typically the largest cross-sectional area of the ellipsoidal initial ablation region; in actual implementation, after each set of initial parameters is obtained, each set of initial parameters can be respectively input into a preset initial simulation model, and simulation is performed, in the simulation process, under each set of initial parameters, the initial long diameter size and the initial short diameter size of the designated section of the initial ablation area corresponding to each second time point respectively are read, for example, a set of initial parameters is set, wherein initial power is 50W, initial ablation duration is 20min, and from the initial ablation duration, the initial long diameter size and the initial short diameter size of the designated section of the initial ablation area corresponding to each simulation 5min,8min,10min,12min and 15min are selected.
And step B22, determining the initial volume size corresponding to the initial ablation region according to the initial cross-sectional size.
Since the initial ablation region may be generally considered as an ellipsoid, after obtaining the initial long diameter size and the initial short diameter size of the designated section of the initial ablation region corresponding to each second time point, the initial volume size of the initial ablation region corresponding to each second time point may be calculated by using an ellipsoid volume calculation formula v=4pi abc/3, where a=initial long diameter size/2, b=initial short diameter size/2, and c=initial short diameter size/2; for example, continuing with the example in the previous step, according to the ellipsoidal volume calculation formula, the initial power is calculated to be 50W, and the initial volume sizes corresponding to the initial ablation regions at the second time points of 5min,8min,10min,12min and 15min respectively.
And B3, acquiring an actual ablation result of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power.
The actual ablation results include: actual cross-sectional dimensions and actual volumetric dimensions; this step B3 can be specifically realized by the following steps B31 and B32:
Step B31, acquiring the actual section size of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under initial power according to each group of initial parameters; wherein, actual ablation area is oval, and actual cross-section size includes: an actual long diameter dimension and an actual short diameter dimension;
in actual operation, in order to ensure the consistency of an isolated pig liver sample, a fresh pig liver on the same day is generally adopted, a microwave ablation needle is inserted into the isolated pig liver about 8cm before ablation, the ablation power is set to be 50W, 5min,8min,10min,12min and 15min are respectively ablated at each power time, four groups of power time data correspond to one isolated pig liver sample, the isolated pig liver sample is cut along a needle track after ablation, then the actual long diameter size and the actual short diameter size of a solidification area of a designated section are read, and the average value of the four groups of actual long diameter sizes and the average value of the actual short diameter size of the same power time data are obtained, so that relatively more accurate values can be obtained.
And step B32, determining the actual volume size corresponding to the actual ablation region according to the actual cross-sectional size.
Since the actual ablation region may be generally considered as an ellipsoid, after the actual long diameter size and the actual short diameter size of the designated section of the actual ablation region corresponding to each second time point are obtained, the actual volume size of the actual ablation region corresponding to each second time point may be calculated by using an ellipsoid volume calculation formula v=4pi abc/3, where a=the actual long diameter size/2, b=the actual short diameter size/2, and c=the actual short diameter size/2; for example, continuing with the example in the previous step, according to the ellipsoidal volume calculation formula, the actual volume size corresponding to the actual ablation region when the initial power is 50W and the second time point is 5min,8min,10min,12min, and 15min respectively can be calculated.
And step B4, comparing the initial simulation result with the actual ablation result to obtain a comparison result, so as to determine a pre-verified simulation model based on the comparison result.
This step B4 can be specifically realized by the following steps 41 to 45:
and step B41, calculating the difference between the initial long diameter size and the actual long diameter size corresponding to each second time point to obtain a plurality of first difference results.
For example, continuing to take the above example as an example, the difference between the simulated initial long diameter size and the actual long diameter size may be calculated when the initial power is 50W and the second time point is 5min, so as to obtain a corresponding first difference result, and so on, so as to obtain the difference between the initial long diameter size and the actual long diameter size when the initial power is 50W and the second time point is 5min,8min,10min,12min,15min, respectively, so as to obtain a plurality of first difference results.
And step B42, calculating the difference between the initial short diameter size and the actual short diameter size corresponding to each second time point to obtain a plurality of second difference results.
For example, continuing to take the above example as an example, the difference between the simulated initial short diameter size and the actual short diameter size may be calculated when the initial power is 50W and the second time point is 5min, so as to obtain a corresponding second difference result, and so on, so as to obtain the difference between the initial short diameter size and the actual short diameter size when the initial power is 50W and the second time point is 5min,8min,10min,12min,15min, respectively, so as to obtain a plurality of second difference results.
And step B43, comparing the difference between the initial volume size and the actual volume size corresponding to each second time point to obtain a plurality of third difference results.
For example, continuing to take the above example as an example, when the initial power is 50W and the second time point is 5min, the difference between the simulated initial volume size and the actual volume size can be calculated, so as to obtain a corresponding third difference result, and so on, when the initial power is 50W and the second time point is 5min,8min,10min,12min,15min, the difference between the corresponding initial volume size and the actual volume size can be respectively obtained, so as to obtain a plurality of third difference results.
Step B44, if at least one first difference result exceeds a preset first threshold range, and/or at least one second difference result exceeds a preset second threshold range, and/or at least one third difference result exceeds a preset third threshold range, receiving a model parameter adjustment instruction for the initial simulation model to adjust the initial simulation model to obtain a second simulation model.
The first threshold range, the second threshold range and the third threshold range may be set according to actual requirements, and in actual implementation, if at least one first difference result exceeds the first threshold range, and/or at least one second difference result exceeds the second threshold range, and/or at least one third difference result exceeds the third threshold range, parameters of the initial simulation model may be considered to be still required to be adjusted, and a user may adjust the parameters of the initial simulation model to obtain an adjusted second simulation model.
And B45, taking the second simulation model as a new initial simulation model, and repeatedly executing the step of acquiring at least one group of initial parameters until each first difference value result belongs to a first threshold range, each second difference value result belongs to a second threshold range, and each third difference value result belongs to a third threshold range, so as to obtain a pre-verified simulation model.
In actual implementation, after the parameters are adjusted to obtain the second simulation model, the step of obtaining at least one group of initial parameters can be repeatedly executed until each first difference result accords with the first threshold range, each second difference result accords with the second threshold range, and each third difference result accords with the third threshold range, so that the simulation model obtained after the parameters are adjusted last time can be considered to be used as a verified simulation model; as shown in fig. 3, the actual ablation section is compared with the simulated section, wherein the actual ablation section is on the left side, and the simulated section is on the right side, so that the two effects are relatively similar.
As shown in a comparison schematic diagram of actual ablation data and simulation data shown in fig. 4 (a), and a comparison schematic diagram of actual ablation data and simulation data shown in fig. 4 (b), comparing the effect of the actual ablation with a verified simulation model, and showing that the specified section of the actual ablation area is 0.18cm in average error, 0.2cm in maximum error and 0.1cm in minimum error compared with the specified section of the ablation simulation using the simulation model through experimental results; the average error of the long diameter is 0.12cm, the maximum error is 0.3cm, and the minimum error is 0; average error of volume is 2.93cm 3 Maximum error of 3.8cm 3 Minimum error of 1.58cm 3 . In order to ensure that the complete excision of the tumor can cut more than 1-2cm to the periphery of the tumor, therefore, the error of ablation simulation and actual ablation by adopting a simulation model is in a trusted range, and the accuracy of simulation is verified by comparing the simulation prediction ablation effect and the actual ablation effect.
And C, determining the simulation volume size corresponding to the simulation ablation area according to the simulation cross section size.
After the simulated cross section size is obtained, the simulated long diameter size and the simulated short diameter size are obtained, and then the simulated volume size corresponding to the ellipsoidal simulated ablation area can be calculated by using an ellipsoidal volume calculation formula.
And D, performing curve fitting on a plurality of simulated long diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a long diameter submodel.
E, curve fitting is carried out on a plurality of simulated short diameter sizes of a plurality of first time points corresponding to each group of simulation parameters, and a short diameter sub-model is obtained;
and F, performing curve fitting on a plurality of simulation volume sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a volume sub-model.
Curve fitting of different functions is performed on a plurality of simulated long diameter sizes, a plurality of simulated short diameter sizes and a plurality of simulated volume sizes of a plurality of first time points corresponding to each set of simulation parameters, wherein the curve fitting is performed on a short diameter fitting function comparison diagram shown in fig. 5 (a), a long diameter fitting function comparison diagram shown in fig. 5 (b), a volume fitting function comparison diagram shown in fig. 5 (c), and fitting curves of the simulated short diameter sizes, the simulated long diameter sizes and the plurality of simulated volume sizes are respectively shown in fig. 5 (a), 5 (b) and 5 (c). The short diameter and the long diameter begin to adopt binomial fitting for 0-20min, and as a result, the short diameter and the long diameter change from the first minute to the second minute too quickly, the fitting degree is not high, and the short diameter R is found 2 = 0.8959, long diameter R 2 = 0.7789, a trigonometric fit was tried, with a short diameter R improved compared to the binomial fit 2 = 0.9554, long diameter R 2 0.8676, the effect is not ideal. The fitting degree is not high because the first minute changes too fast, and the fitting effect is better in 1-20min by comparing the fitting effect in 0-20min through fitting a curve in 1-20min, so that the simulated long diameter size and the simulated short diameter size are fitted by adopting a binomial function in 1-20 min. The volume is compared with binomial and trinomial fitting from 0 to 20min, and the effect is not very different, so that a simpler binomial fitting is adopted, as shown in fig. 6, which is a schematic diagram of a final ablation prediction model provided by the embodiment of the present invention, and from left to right, the two are respectively the volumesProduct and power, time 3D model, short diameter and power, time 3D model, and long diameter and power, time 3D model.
In the embodiment, microwave ablation power of 30W-80W is adopted, data are ablated for 1-20min, data fitting of different functions is carried out through simulated short diameter size, simulated long diameter size and simulated volume size of an ablation area, the fitting effect is compared, and the highest fitting degree is selected, namely R 2 Closest to 1, the function closest to the actual ablation process is used as the 3D ablation prediction model. And comparing the actual ablation area size with the ablation area size predicted by the ablation prediction model, determining the maximum absolute error, the minimum absolute error and the average error, and verifying the reliability of the 3D ablation prediction model.
And fifthly, calculating microwave ablation parameters matched with the area to be ablated based on the first sub-model, the second sub-model and the third sub-model.
And performing joint calculation on the first sub-model, the second sub-model and the third sub-model to obtain microwave ablation power and microwave ablation time, and performing ablation on an area to be ablated according to the microwave ablation power and the microwave ablation time to obtain an ablation result matched with the area to be ablated.
A comparative schematic of ablation data is shown in fig. 7, which is a further verification of whether the above 3D ablation effect prediction model can truly provide assistance for clinical use. At present, commercial microwave ablation equipment manufacturers can provide different powers and time ablation effect comparison tables for doctors, and provide ablation scheme references for the doctors. According to the published material search, in the specification of the ablation needle of the head enterprise in the microwave ablation field of A and B, the provided 2.0mm needle type ablation effect data are compared with the 3D ablation prediction model data established by the invention by selecting the data of 60W, 5min,8min,10min,12min and 15min of ablation, and determining errors, so that the accuracy and the practicability of the 3D ablation prediction model are further verified. In the figure, the dark filled part is the ablation short diameter range provided by the manufacturer B, the light filled part is the ablation short diameter range provided by the manufacturer A, and the model prediction value can be found to be equivalent to the ablation range of the current commercial equipment, wherein the model prediction value is relatively more accurate within the range of 300s-600s, and the deviation above 600s is relatively larger, but is also within the practical allowable range.
The simulation data of the 3D ablation prediction model is compared with the data of the first manufacturer, and the error of the short diameter is as follows: the average error was 0.84cm, the maximum error was 1.8cm, and the minimum error was 0.2cm. The simulation data of the 3D ablation prediction model is compared with the data of manufacturer B, and the error of the short diameter is as follows: the average error is 0.46cm, the maximum error is 1.1cm, and the minimum error is 0.1cm; the errors were all below 2cm and therefore acceptable, while it can be seen that the deviation in effect of the ablation predictive model after 10min was relatively greater than 1cm. In the figure, the short diameter is given by an upper limit and a lower limit, and the average value of the upper limit and the lower limit can be taken and compared with the result of the 3D ablation prediction model.
The simulation data of the 3D ablation prediction model is compared with the data of the first manufacturer, and the error of the long diameter is as follows: the average error was 0.5cm, the maximum error was 1.0cm, and the minimum error was 0.3cm. The simulation data of the 3D ablation prediction model is compared with the data of manufacturer B, and the error of the long diameter is as follows: the average error is 0.38cm, the maximum error is 0.6cm, and the minimum error is 0.2cm; the error is below 2cm, so that the method is acceptable, and meanwhile, the effect of the long diameter of the model is more approximate to actual application data than that of the short diameter. The volume may be calculated from the long diameter and the short diameter. Therefore, the 3D ablation prediction model provided by the invention has clinical practical value, particularly has higher accuracy and higher credibility in the period of 3-10min of common microwave ablation, and ensures the effectiveness and practicability of the verification model.
The method for determining the microwave ablation parameters discloses a simulation-based microwave ablation 3D ablation prediction model, and a set of 3D ablation prediction model is established according to the microwave ablation simulation model. Firstly, establishing a set of microwave ablation and simulation data acquisition system, and acquiring actual ablation data and simulation ablation data of different ablation powers and times; and comparing actual data with simulation data errors, verifying the effectiveness of a simulation model, obtaining a 30W-80W 3D ablation prediction model through a data fitting method, and finally verifying the accuracy of the model. The invention solves the problem that in the field of microwave ablation, a doctor can only select the pain point of the ablation power and time through experience before operation, and has important significance for making an operation scheme before operation. Finding the corresponding relation between the size of the ablation area and the ablation power and time of microwave ablation has great significance to clinic.
The microwave ablation simulation is verified as a method for effectively simulating an actual ablation process, and an ablation effect is simulated by calculating a tissue electromagnetic field and a temperature field by using an electromagnetic wave transmission model and a Pennes equation. The COMSOL software is specialized multimode standing model finite element analysis software for realizing ablation simulation process simulation. The ablation simulation technology is subjected to a static-to-dynamic optimization process, and the optimized dynamic ablation simulation technology is adopted to simulate the ablation process by considering the dielectric characteristics and thermophysical parameters of tissues and perform experimental verification. And then, on the basis of the verified simulation model, exploring the size relation of different power time and ablation areas, and establishing a set of simulation-based ablation area 3D prediction model.
The mode provides a set of 3D prediction models of ablation areas based on simulation, establishes a prediction relation model of ablation power, time and final ablation effect, and particularly provides a relation model of the most direct ablation power (the most common power) of 30W-80W, wherein the relation model of the ablation power, the ablation time and the ablation area size can directly calculate the ablation area size at any time under the powers, and is not limited to individual power and time. Compared with the ablation effect reference table provided by the current commercial microwave ablation manufacturer, the size of the ablation area provided by the model is more accurate, and is not only a approximate range.
The size of the ablation region is predicted by building an ablation prediction model. In microwave ablation operation, certain ablation power and time are set, and the size of an ablated region is predicted, so that the problem that a doctor cannot know how much power and time to set before operation and can completely ablate tumors without excessively ablating, and normal tissues are damaged can be solved.
The method utilizes a simulation technology to obtain the size of an ablation area under each ablation dose (power and time), so that the size of the ablation area under all power time is obtained through data modeling. When the problem of setting ablation parameters by doctors is solved by the existing commercial products, a set of table is provided, ablation sizes of 55w,60w,70w,80w,3min,5min,8min,10min,12min and 15min corresponding to in-vitro pig liver experiments are displayed, and reference is made to the doctors, but the sizes of ablation areas under other parameters still cannot be confirmed. And the prior art can only provide a reference of approximate scope through actual isolated pig liver experiments. The experiment is restricted by various conditions such as environment, pig liver and the like, the given ablation area can only be a rough range, and the experiment cannot traverse the data of all ablation parameters. According to the scheme, based on a simulation technology, mathematical calculation is performed by establishing a 3D ablation prediction model, so that a more accurate simulation result can be obtained.
An embodiment of the present invention provides a device for determining microwave ablation parameters, as shown in fig. 8, where the device includes: an acquisition module 80, configured to acquire size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; the output module 81 is configured to input the size data into a preset ablation prediction model, so as to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
The device for determining the microwave ablation parameters acquires the size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data, and volume data of the region to be ablated; inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time. According to the device, the microwave ablation power and the microwave ablation time matched with the area to be ablated can be obtained by inputting the size data of the area to be ablated into the ablation prediction model, so that the ablation efficiency can be improved.
Further, the ablation prediction model includes: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; the output module 81 is further configured to: inputting the long-diameter data into a long-diameter submodel to obtain a first submodel; inputting the short-diameter data into a short-diameter submodel to obtain a second submodel; inputting the volume data into a volume sub-model to obtain a third sub-model; and calculating microwave ablation parameters matched with the region to be ablated based on the first sub-model, the second sub-model and the third sub-model.
Further, the apparatus further comprises an ablation prediction model determination module, the ablation prediction model is determined by the ablation prediction model determination module, and the ablation prediction model determination module is used for: obtaining a plurality of groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration; inputting the group of simulation parameters into a pre-verified simulation model aiming at each group of simulation parameters, and reading simulation cross section sizes of designated cross sections of simulation ablation areas corresponding to a plurality of first time points in simulation ablation time under simulation power; wherein the simulated ablation area is ellipsoidal; the simulated cross-sectional dimensions include: simulating a long diameter dimension and a short diameter dimension; determining the simulation volume size corresponding to the simulation ablation area according to the simulation cross section size; performing curve fitting on a plurality of simulated long diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a long diameter sub-model; performing curve fitting on a plurality of simulated short diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a short diameter sub-model; and performing curve fitting on a plurality of simulation volume sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain a volume sub-model.
Further, the device also comprises a simulation model determining module, wherein the simulation model which is verified in advance is determined by the simulation model determining module, and the simulation model determining module is used for: acquiring at least one set of initial parameters; wherein each set of initial parameters comprises: an initial power and an initial ablation duration; inputting the initial parameters into a preset initial simulation model aiming at each initial parameter group, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation duration under initial power; acquiring an actual ablation result of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power; and comparing the initial simulation result with the actual ablation result to obtain a comparison result so as to determine a pre-verified simulation model based on the comparison result.
Further, the initial simulation result includes: an initial cross-sectional dimension and an initial volumetric dimension; the simulation model determination module is further configured to: inputting the initial parameters into a preset initial simulation model aiming at each initial parameter, and reading initial section sizes of designated sections of initial ablation areas corresponding to a plurality of second time points in initial ablation time under initial power; wherein the initial ablation region is ellipsoidal; the initial cross-sectional dimensions include: an initial long diameter dimension and an initial short diameter dimension; and determining the initial volume size corresponding to the initial ablation region according to the initial cross-sectional size.
Further, the actual ablation results include: actual cross-sectional dimensions and actual volumetric dimensions; the simulation model determination module is further configured to: acquiring the actual section size of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power according to each group of initial parameters; wherein, actual ablation area is oval, and actual cross-section size includes: an actual long diameter dimension and an actual short diameter dimension; and determining the actual volume size corresponding to the actual ablation region according to the actual cross-sectional size.
Further, the simulation model determining module is further configured to: calculating the difference value between the initial long diameter size and the actual long diameter size corresponding to each second time point to obtain a plurality of first difference value results; calculating the difference value between the initial short diameter size and the actual short diameter size corresponding to each second time point to obtain a plurality of second difference value results; comparing the difference value between the initial volume size and the actual volume size corresponding to each second time point to obtain a plurality of third difference value results; if at least one first difference result exceeds a preset first threshold range, and/or at least one second difference result exceeds a preset second threshold range, and/or at least one third difference result exceeds a preset third threshold range, receiving a model parameter adjustment instruction for the initial simulation model to adjust the initial simulation model to obtain a second simulation model; and taking the second simulation model as a new initial simulation model, and repeatedly executing the step of acquiring at least one group of initial parameters until each first difference result belongs to a first threshold range, each second difference result belongs to a second threshold range, and each third difference result belongs to a third threshold range, so as to obtain a pre-verified simulation model.
The implementation principle and the generated technical effects of the device for determining the microwave ablation parameters provided by the embodiment of the invention are the same as those of the embodiment of the method for determining the microwave ablation parameters, and for the sake of brief description, reference may be made to corresponding contents in the embodiment of the method for determining the microwave ablation parameters where the embodiment of the device for determining the microwave ablation parameters is not mentioned.
The embodiment of the present invention further provides an electronic device, referring to fig. 9, where the electronic device includes a processor 130 and a memory 131, where the memory 131 stores machine executable instructions that can be executed by the processor 130, and the processor 130 executes the machine executable instructions to implement the method for determining microwave ablation parameters described above.
Further, the electronic device shown in fig. 9 further includes a bus 132 and a communication interface 133, and the processor 130, the communication interface 133, and the memory 131 are connected through the bus 132.
The memory 131 may include a high-speed random access memory (RAM, random Access Memory), and may further include a non-volatile memory (non-volatile memory), such as at least one magnetic disk memory. The communication connection between the system network element and at least one other network element is implemented via at least one communication interface 133 (which may be wired or wireless), and may use the internet, a wide area network, a local network, a metropolitan area network, etc. Bus 132 may be an ISA bus, a PCI bus, an EISA bus, or the like. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, only one bi-directional arrow is shown in fig. 9, but not only one bus or one type of bus.
The processor 130 may be an integrated circuit chip with signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuitry in hardware or instructions in software in processor 130. The processor 130 may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU for short), a network processor (Network Processor, NP for short), etc.; but also digital signal processors (Digital Signal Processor, DSP for short), application specific integrated circuits (Application Specific Integrated Circuit, ASIC for short), field-programmable gate arrays (Field-Programmable Gate Array, FPGA for short) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The disclosed methods, steps, and logic blocks in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be embodied directly in the execution of a hardware decoding processor, or in the execution of a combination of hardware and software modules in a decoding processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in the memory 131, and the processor 130 reads the information in the memory 131, and in combination with its hardware, performs the steps of the method of the foregoing embodiment.
The embodiment of the invention also provides a machine-readable storage medium, which stores machine-executable instructions that, when being called and executed by a processor, cause the processor to implement the method for determining microwave ablation parameters, and the specific implementation can be referred to the method embodiment and will not be repeated herein.
The method, the device and the computer program product of the electronic device for determining the microwave ablation parameters provided by the embodiment of the invention comprise a computer readable storage medium storing program codes, and the instructions included in the program codes can be used for executing the method described in the method embodiment, and specific implementation can be referred to the method embodiment and will not be repeated here.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (10)
1. A method of determining microwave ablation parameters, the method comprising:
acquiring size data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data and volume data of the region to be ablated;
inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
2. The method of claim 1, wherein the ablation predictive model comprises: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; the step of inputting the size data into a preset ablation prediction model to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model comprises the following steps of:
Inputting the long-diameter data into the long-diameter submodel to obtain a first submodel;
inputting the short-diameter data into the short-diameter submodel to obtain a second submodel;
inputting the volume data into the volume sub-model to obtain a third sub-model;
and calculating microwave ablation parameters matched with the region to be ablated based on the first sub-model, the second sub-model and the third sub-model.
3. The method of claim 2, wherein the ablation predictive model is determined by:
obtaining a plurality of groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration;
inputting the group of simulation parameters into a pre-verified simulation model aiming at each group of simulation parameters, and reading simulation cross section sizes of designated cross sections of simulation ablation areas corresponding to a plurality of first time points in the simulation ablation duration under the simulation power; wherein the simulated ablation area is ellipsoidal; the simulated cross-sectional dimensions include: simulating a long diameter dimension and a short diameter dimension;
determining a simulation volume size corresponding to the simulation ablation area according to the simulation cross section size;
Performing curve fitting on a plurality of simulated long diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain the long diameter submodel;
performing curve fitting on a plurality of simulated short diameter sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain the short diameter submodel;
and performing curve fitting on a plurality of simulation volume sizes of a plurality of first time points corresponding to each group of simulation parameters to obtain the volume sub-model.
4. A method according to claim 3, wherein the pre-validated simulation model is determined by:
acquiring at least one set of initial parameters; wherein each set of initial parameters comprises: an initial power and an initial ablation duration;
inputting the initial parameters into a preset initial simulation model aiming at each initial parameter group, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in the initial ablation duration under the initial power;
acquiring an actual ablation result of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power;
and comparing the initial simulation result with the actual ablation result to obtain a comparison result so as to determine a pre-verified simulation model based on the comparison result.
5. The method of claim 4, wherein the initial simulation results comprise: an initial cross-sectional dimension and an initial volumetric dimension; inputting the initial parameters into a preset initial simulation model aiming at each initial parameter, and acquiring initial simulation results of designated sections of initial ablation areas corresponding to a plurality of second time points in the initial ablation duration under the initial power, wherein the steps comprise:
inputting the initial parameters into a preset initial simulation model aiming at each initial parameter, and reading initial section sizes of designated sections of initial ablation areas corresponding to a plurality of second time points in the initial ablation duration under the initial power; wherein the initial ablation region is ellipsoidal; the initial cross-sectional dimensions include: an initial long diameter dimension and an initial short diameter dimension;
and determining the initial volume size corresponding to the initial ablation region according to the initial cross-sectional size.
6. The method of claim 5, wherein the actual ablation results comprise: actual cross-sectional dimensions and actual volumetric dimensions;
the step of obtaining the actual ablation result of the designated section of the actual ablation area of the preset sample to be ablated corresponding to each second time point under the initial power comprises the following steps:
Acquiring the actual section size of a designated section of an actual ablation area of a preset sample to be ablated corresponding to each second time point under the initial power according to each group of initial parameters; wherein the actual ablation region is elliptical, and the actual cross-sectional dimension comprises: an actual long diameter dimension and an actual short diameter dimension;
and determining the actual volume size corresponding to the actual ablation region according to the actual cross-sectional size.
7. The method of claim 6, wherein comparing the initial simulation result and the actual ablation result to obtain a comparison result to determine a pre-validated simulation model based on the comparison result comprises:
calculating the difference value between the initial long diameter size and the actual long diameter size corresponding to each second time point to obtain a plurality of first difference value results;
calculating the difference value between the initial short diameter size and the actual short diameter size corresponding to each second time point to obtain a plurality of second difference value results;
comparing the difference value between the initial volume size and the actual volume size corresponding to each second time point to obtain a plurality of third difference value results;
If at least one first difference result exceeds a preset first threshold range, and/or at least one second difference result exceeds a preset second threshold range, and/or at least one third difference result exceeds a preset third threshold range, receiving a model parameter adjustment instruction aiming at the initial simulation model to adjust the initial simulation model to obtain a second simulation model;
and taking the second simulation model as a new initial simulation model, and repeatedly executing the step of acquiring at least one group of initial parameters until each first difference result belongs to the first threshold range, each second difference result belongs to the second threshold range, and each third difference result belongs to the third threshold range, so as to obtain a pre-verified simulation model.
8. A device for determining microwave ablation parameters, the device comprising:
the acquisition module is used for acquiring the size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the size data includes: long diameter data, short diameter data and volume data of the region to be ablated;
The output module is used for inputting the size data into a preset ablation prediction model so as to output microwave ablation parameters matched with the area to be ablated through the ablation prediction model; wherein the microwave ablation parameters include: microwave ablation power and microwave ablation time.
9. An electronic device comprising a processor and a memory, the memory storing machine executable instructions executable by the processor, the processor executing the machine executable instructions to implement the method of determining microwave ablation parameters of any of claims 1-7.
10. A machine-readable storage medium storing machine-executable instructions that, when invoked and executed by a processor, cause the processor to implement the method of determining microwave ablation parameters according to any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311618626.4A CN117637173A (en) | 2023-11-29 | 2023-11-29 | Method and device for determining microwave ablation parameters and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311618626.4A CN117637173A (en) | 2023-11-29 | 2023-11-29 | Method and device for determining microwave ablation parameters and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117637173A true CN117637173A (en) | 2024-03-01 |
Family
ID=90037165
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311618626.4A Pending CN117637173A (en) | 2023-11-29 | 2023-11-29 | Method and device for determining microwave ablation parameters and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117637173A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111685870A (en) * | 2020-06-23 | 2020-09-22 | 南京诺源医疗器械有限公司 | Microwave ablation carbonization regulation and control method for liver |
| WO2021258639A1 (en) * | 2020-06-23 | 2021-12-30 | 南京诺源医疗器械有限公司 | Therapeutic effect evaluation method in microwave ablation based on simulated temperature field |
| EP4011308A1 (en) * | 2020-12-10 | 2022-06-15 | Koninklijke Philips N.V. | Heat distribution model databases for planning thermal ablation |
| CN116115328A (en) * | 2022-12-29 | 2023-05-16 | 融和医疗科技(浙江)有限公司 | Object ablation system, control method, device, medium and electronic equipment |
-
2023
- 2023-11-29 CN CN202311618626.4A patent/CN117637173A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111685870A (en) * | 2020-06-23 | 2020-09-22 | 南京诺源医疗器械有限公司 | Microwave ablation carbonization regulation and control method for liver |
| WO2021258639A1 (en) * | 2020-06-23 | 2021-12-30 | 南京诺源医疗器械有限公司 | Therapeutic effect evaluation method in microwave ablation based on simulated temperature field |
| EP4011308A1 (en) * | 2020-12-10 | 2022-06-15 | Koninklijke Philips N.V. | Heat distribution model databases for planning thermal ablation |
| CN116115328A (en) * | 2022-12-29 | 2023-05-16 | 融和医疗科技(浙江)有限公司 | Object ablation system, control method, device, medium and electronic equipment |
Non-Patent Citations (1)
| Title |
|---|
| 赵金哲 等: "微波消融有效消融体积模型实验研究", 中国生物医学工程学报, vol. 33, no. 1, 28 February 2014 (2014-02-28), pages 52 - 55 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109498155B (en) | Rapid planning system and method in ablation treatment operation | |
| Ye et al. | The accuracy of computer image-guided template for mandibular angle ostectomy | |
| JP6602299B2 (en) | Efficient treatment planning trade-off analysis | |
| CN105825073B (en) | A kind of online radiotherapy planning quality control system | |
| US10607738B2 (en) | System and method for minimally invasive thermal ablation treatment planning | |
| KR20130140067A (en) | Systems and methods for planning image-guided interventional procedures | |
| Singh et al. | Numerical study to establish relationship between coagulation volume and target tip temperature during temperature-controlled radiofrequency ablation | |
| Rijnen et al. | Quality and comfort in head and neck hyperthermia: A redesign according to clinical experience and simulation studies | |
| Deuflhard et al. | Mathematical cancer therapy planning in deep regional hyperthermia | |
| CN117297758B (en) | Method and device for determining microwave ablation parameters and electronic equipment | |
| Wang et al. | Numerical evaluation of ablation zone under different tip temperatures during radiofrequency ablation | |
| CN117637173A (en) | Method and device for determining microwave ablation parameters and electronic equipment | |
| Radosevic et al. | Short pulsed microwave ablation: computer modeling and ex vivo experiments | |
| Masnok et al. | Catheter contact area strongly correlates with lesion area in radiofrequency cardiac ablation: an ex vivo porcine heart study | |
| CN116936112B (en) | Regression model training and ablation planning method, device, equipment and medium | |
| CN113693708A (en) | Radio frequency output adjusting method and device of radio frequency ablation equipment and computer storage medium | |
| Altrogge et al. | Towards optimization of probe placement for radio-frequency ablation | |
| Singh et al. | Numerical investigation of convective cooling in minimizing skin burns during radiofrequency ablation of breast tumor | |
| CN113952023B (en) | Dynamic simulation method, device, storage medium and equipment for ablation needle | |
| Ziegle et al. | RF-ablation pattern shaping employing switching channels of dual bipolar needle electrodes: ex vivo results | |
| CN113974820A (en) | Simulated ablation method, device, storage medium and equipment based on residual fitting | |
| CN114983559A (en) | Data processing method, device, system, equipment and storage medium for radio frequency ablation | |
| EP4061265B1 (en) | Heat distribution model databases for planning thermal ablation | |
| Nour et al. | Thermal damage modeling analysis and validation during treatment of tissue tumors | |
| Zaridze et al. | Electro-thermal computational suit for investigation of rf power absorption and associated temperature change in human body |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |