CN117297758B - 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 PDF

Info

Publication number
CN117297758B
CN117297758B CN202311606412.5A CN202311606412A CN117297758B CN 117297758 B CN117297758 B CN 117297758B CN 202311606412 A CN202311606412 A CN 202311606412A CN 117297758 B CN117297758 B CN 117297758B
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.)
Active
Application number
CN202311606412.5A
Other languages
Chinese (zh)
Other versions
CN117297758A (en
Inventor
钱露
蔡惠明
钱志余
王银芳
张舒
杜凯
曹勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Nuoyuan Medical Devices Co Ltd
Original Assignee
Nanjing Nuoyuan Medical Devices Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Nuoyuan Medical Devices Co Ltd filed Critical Nanjing Nuoyuan Medical Devices Co Ltd
Priority to CN202311606412.5A priority Critical patent/CN117297758B/en
Publication of CN117297758A publication Critical patent/CN117297758A/en
Application granted granted Critical
Publication of CN117297758B publication Critical patent/CN117297758B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B18/18—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by applying electromagnetic radiation, e.g. microwaves
    • A61B18/1815—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by applying electromagnetic radiation, e.g. microwaves using microwaves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1072—Measuring physical dimensions, e.g. size of the entire body or parts thereof measuring distances on the body, e.g. measuring length, height or thickness
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1073—Measuring volume, e.g. of limbs
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00571—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for achieving a particular surgical effect
    • A61B2018/00577—Ablation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00636—Sensing and controlling the application of energy
    • A61B2018/00696—Controlled or regulated parameters
    • A61B2018/00702—Power or energy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00636—Sensing and controlling the application of energy
    • A61B2018/00696—Controlled or regulated parameters
    • A61B2018/00761—Duration

Abstract

The invention provides a method and a device for determining microwave ablation parameters and electronic equipment, wherein the method comprises the following steps: acquiring section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; calculating volume data of the area to be ablated according to the long diameter data and the short diameter data; inputting the section size data and the volume 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

Method and device for determining microwave ablation parameters and electronic equipment
Technical Field
The present invention relates to the field of microwave ablation technology, 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 section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the above cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; calculating volume data of the area to be ablated according to the long diameter data and the short diameter data; inputting the section size data and the volume data into a preset ablation prediction model to output microwave ablation parameters matched with the region to be ablated through the ablation prediction model; wherein, the microwave ablation parameters include: microwave ablation power and microwave ablation time; the ablation prediction model is constructed based on preset multiple groups of simulation parameters and a preset verified simulation model.
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 section size data and the volume data into a preset ablation prediction model to output microwave ablation parameters matched with the region to be ablated through the ablation prediction model comprises the following steps: 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.
Further, the above ablation prediction model is determined by: acquiring the multiple groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration; inputting the group of simulation parameters into the verified simulation model aiming at each group of simulation parameters, and reading simulation section sizes of designated 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.
Further, the verified simulation model is determined by the following method: 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 the 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 for each initial parameter, and obtaining initial simulation results of designated sections of the initial ablation region 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 described above 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 result includes: 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 area is elliptical, and the actual cross-sectional dimension includes: an actual long diameter dimension and an actual short diameter dimension; and determining the actual volume size corresponding to the actual ablation area 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 the verified simulation model based on the comparison result includes: 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; 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; 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; if at least one of the first difference results exceeds a preset first threshold range, and/or at least one of the second difference results exceeds a preset second threshold range, and/or at least one of the third difference results 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 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, thereby obtaining the verified simulation model.
The invention provides a device for determining microwave ablation parameters, which comprises: the acquisition module is used for acquiring the section size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the above cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; the volume calculation module is used for calculating the volume data of the area to be ablated according to the long diameter data and the short diameter data; the output module is used for inputting the section size data and the volume 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 ablation prediction model is constructed based on preset multiple groups of simulation parameters and a preset verified simulation model.
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 method and device for determining the microwave ablation parameters and the electronic equipment provided by the embodiment of the invention comprise the following steps: acquiring section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; calculating to obtain volume data of the area to be ablated according to the long diameter data and the short diameter data; inputting the section size data and the volume 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 (a) is a schematic diagram of model verification of a short diameter submodel and a long diameter submodel provided by an embodiment of the present invention;
FIG. 6 (b) is a schematic diagram of model verification of a volumetric model according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of comparison of ablation data according to an embodiment of the present invention;
FIG. 8 (a) is a schematic diagram of a 2D short diameter submodel relationship between 30W and 80W for microwave ablation power provided by an embodiment of the present invention;
FIG. 8 (b) is a schematic diagram of a 2D long diameter submodel relationship with microwave ablation power between 30W and 80W according to an embodiment of the present invention;
FIG. 8 (c) is a schematic diagram of a 2D volumetric sub-model relationship with microwave ablation power of 30W-80W according to an embodiment of the present invention
Fig. 9 is a schematic structural diagram of a device for determining microwave ablation parameters according to an embodiment of the present invention;
fig. 10 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Reference numerals illustrate: 20-a microwave ablation device; 21-a microwave ablation needle; 23-isolated pig liver samples; 24-PC; 25-ablation simulation software interface; 91-an acquisition module; 92-a volume calculation module; 93-an output module; 130-a processor; 131-memory; 132-bus; 133-communication interface.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
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.
Example 1
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 S101: acquiring section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the above cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated.
In this embodiment, the area to be ablated may be an area including a tumor or the like that 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 and short diameter 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 S102: and calculating the volume data of the area to be ablated according to the long diameter data and the short diameter data.
Here, according to a preset ellipsoidal volume calculation formula, the volume data of the region to be ablated can be calculated according to the long diameter data and the short diameter data.
Step S103: inputting the section size data and the volume data into a preset ablation prediction model to output microwave ablation parameters matched with the region to be ablated through the ablation prediction model; wherein, the microwave ablation parameters include: microwave ablation power and microwave ablation time; the ablation prediction model is constructed based on preset multiple groups of simulation parameters and a preset verified simulation model.
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 provided by the embodiment of the invention comprises the following steps: acquiring section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; calculating to obtain volume data of the area to be ablated according to the long diameter data and the short diameter data; inputting the section size data and the volume 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.
Example 2
On the basis of the method for determining the microwave ablation parameters shown in fig. 1, the embodiment of the invention also provides another method for determining the microwave ablation parameters, which is realized on the basis of the method.
In this embodiment, the ablation prediction model includes: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; in actual implementation, if the microwave ablation power is 50W, the mathematical models corresponding to the long-diameter submodel, the short-diameter submodel and the volume submodel respectively can be expressed as follows:
Short diameter submodel: l (L) Short diameter =0.716lnt+0.2179，R 2 =0.9747；
Long diameter submodel: l (L) Long diameter =0.8974lnt−2.336，R 2 =0.995；
Volume submodel:；
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; 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 and short diameter data of the region to be ablated.
And step two, calculating to obtain volume data of the area to be ablated according to the long-diameter data and the short-diameter data.
And thirdly, inputting the long-diameter data into the long-diameter submodel to obtain a first submodel.
And step four, inputting the short-diameter data into the short-diameter submodel to obtain a second submodel.
And fifthly, inputting the volume data into the volume sub-model to obtain a third sub-model.
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.
And B2, inputting the initial parameters into a preset initial simulation model according to each initial parameter, and acquiring initial simulation results of the designated section of the initial ablation region corresponding to a plurality of second time points in the 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 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 are obtained, 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, the error of ablation simulation and actual ablation by adopting a simulation model is in a trusted rangeAnd in addition, the simulation accuracy 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 are found to change too quickly from the first minute to the second minute, the fitting degree is not high, the short diameter r2= 0.9296, the long diameter r2=0.74, the trend of the logarithmic curve is found to be closer to the ablation effect by fitting a curve for 1-20min to compare binomial fitting with logarithmic fitting, so that the short diameter and the long diameter are fitted by adopting logarithmic function, the volume is compared with binomial fitting and binomial fitting from 0-20min, the effect is not very different, and simpler binomial fitting is adopted.
Further, fig. 6 (a) is a schematic diagram of model verification of a short diameter submodel and a long diameter submodel according to an embodiment of the present invention. Comparing actual experimental data with predicted data of the 2D short diameter submodel, finding that the average error of the short diameter is 0.08cm, the maximum error is 0.2cm, the minimum error is 0, and finding that the average error of the long diameter is 0.12cm, the maximum error is 0.3cm, and the minimum error is 0;
further, fig. 6 (b) is a schematic diagram of model verification of a volumetric sub-model according to an embodiment of the present invention. As seen in FIG. 6 (b), the average error in volume was 1.88cm 3 The maximum error is 3.88cm3, and the minimum error is 0.26cm 3 . Errors are all within 1cm, and are acceptable, so that the validity of the 2D model is verified.
According to the embodiment, microwave ablation power is 30W-80W, data are ablated for 1-20min, data fitting of different functions is conducted through simulated short diameter size, simulated long diameter size and simulated volume size of an ablation area, the fitting effect is compared, and the function closest to the actual ablation process is selected as a 2D ablation prediction model, wherein the fitting effect is highest. 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 2D ablation prediction model.
In this example, 2D ablation effect prediction models of 30w,40w,50w,60w,70w,80w were constructed using the method described above, with each submodel referring to table 1 below, and table 1 above is a submodel list:
table 1:
and step six, calculating to obtain 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.
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 2D 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 2D 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 2D ablation prediction model are further verified. In fig. 7, the dark filled portion is the short diameter range of ablation provided by manufacturer b, the light filled portion is the short diameter range provided by manufacturer a, and the model predictive value can be found to be comparable to the ablation range of current commercial equipment.
The simulation data of the 2D 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 2D 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 fig. 7, the short path is given by the upper and lower limits, and the average value of the upper and lower limits can be taken and compared with the result of the 2D ablation prediction model.
The simulation data of the 2D 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 2D 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 2D 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 2D ablation prediction model, and a set of 2D 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 2D 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.
Further, fig. 8 (a) is a schematic diagram of a 2D short diameter submodel relationship of microwave ablation power between 30W and 80W according to an embodiment of the present invention; FIG. 8 (b) is a schematic diagram of a 2D long diameter submodel relationship with microwave ablation power between 30W and 80W according to an embodiment of the present invention; fig. 8 (c) is a schematic diagram of a 2D volumetric sub-model relationship of microwave ablation power between 30W and 80W according to an embodiment of the present invention.
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 2D prediction model.
The mode provides a set of simulation-based 2D prediction model of the ablation area, establishes a prediction relation model of the 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 2D ablation prediction model, so that a more accurate simulation result can be obtained.
Example 3
In this embodiment, fig. 9 is a schematic structural diagram of a device for determining microwave ablation parameters according to an embodiment of the present invention.
As can be seen in fig. 9, the device comprises:
an acquiring module 91, configured to acquire cross-sectional dimension data of an area to be ablated; wherein the area to be ablated is ellipsoidal; the above cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated.
And a volume calculating module 92, configured to calculate volume data of the region to be ablated according to the long diameter data and the short diameter data.
The output module 93 is configured to input the section size data and the volume data into a preset ablation prediction model, so as to output microwave ablation parameters matched with the region to be ablated through the ablation prediction model; wherein, the microwave ablation parameters include: microwave ablation power and microwave ablation time; the ablation prediction model is constructed based on preset multiple groups of simulation parameters and a preset verified simulation model.
The acquisition module 91, the volume calculation module 92, and the output module 93 are sequentially connected.
In one embodiment, the ablation prediction model includes: a long diameter sub-model, a short diameter sub-model, and a volume sub-model; the output module 93 is further configured to input 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.
In one embodiment, the apparatus further comprises an ablation prediction model building module; the ablation prediction model construction module is used for acquiring the multiple groups of simulation parameters; wherein each set of simulation parameters includes: simulation power and simulation ablation duration; inputting the group of simulation parameters into the verified simulation model aiming at each group of simulation parameters, and reading simulation section sizes of designated 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.
In one embodiment, the apparatus further includes a simulation model determination module; the simulation model determining module is used for acquiring at least one group 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 the verified simulation model based on the comparison result.
In one embodiment, the initial simulation result includes: an initial cross-sectional dimension and an initial volumetric dimension; the simulation model determining module is further configured to input, for each set of initial parameters, the set of initial parameters into a preset initial simulation model, and read initial cross-sectional dimensions of designated cross-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 described above 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.
In one embodiment, the actual ablation result includes: actual cross-sectional dimensions and actual volumetric dimensions; the simulation model determining module is further configured to obtain, for each set of initial parameters, an actual cross-section size of a designated cross-section of an actual ablation area of a preset sample to be ablated, where the designated cross-section corresponds to each second time point under the initial power; wherein the actual ablation area is elliptical, and the actual cross-sectional dimension includes: an actual long diameter dimension and an actual short diameter dimension; and determining the actual volume size corresponding to the actual ablation area according to the actual cross-sectional size.
In one embodiment, the simulation model determining module is further configured to calculate a difference between the initial long diameter size and the actual long diameter size corresponding to each of the second time points, to obtain a plurality of first difference results; 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; 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; if at least one of the first difference results exceeds a preset first threshold range, and/or at least one of the second difference results exceeds a preset second threshold range, and/or at least one of the third difference results 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 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, thereby obtaining the 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.
Example 3
The embodiment of the present invention further provides an electronic device, referring to fig. 10, 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. 10 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. 10, but not only one bus or 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 (7)

1. A method of determining microwave ablation parameters, the method comprising:
acquiring section size data of a region to be ablated; wherein the area to be ablated is ellipsoidal; the cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; the long diameter data is used for indicating the longest diameter of the area to be ablated; the short diameter data is used for indicating the shortest diameter of the area to be ablated;
calculating volume data of the area to be ablated according to the long diameter data and the short diameter data;
acquiring a plurality of preset 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 verified simulation model aiming at each group of simulation parameters, and reading simulation section sizes of designated 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; the verified simulation model is determined by the following method:
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;
comparing the initial simulation result with the actual ablation result to obtain a comparison result so as to determine the verified simulation model based on the comparison result;
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 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;
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;
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;
based on the first sub-model, the second sub-model and the third sub-model, calculating to obtain microwave ablation parameters matched with the region to be ablated, wherein the microwave ablation parameters comprise: microwave ablation power and microwave ablation time.
2. The method of claim 1, 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.
3. The method of determining microwave ablation parameters according to claim 2, wherein 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 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.
4. The method of determining microwave ablation parameters according to claim 3, wherein comparing the initial simulation result and the actual ablation result to obtain a comparison result, and determining the 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 the verified simulation model.
5. A device for determining microwave ablation parameters, the device comprising:
the acquisition module is used for acquiring the section size data of the area to be ablated; wherein the area to be ablated is ellipsoidal; the cross-sectional dimension data includes: long diameter data and short diameter data of the region to be ablated; the long diameter data is used for indicating the longest diameter of the area to be ablated; the short diameter data is used for indicating the shortest diameter of the area to be ablated;
the volume calculation module is used for calculating the volume data of the area to be ablated according to the long diameter data and the short diameter data;
the output module is used for acquiring a plurality of preset 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 verified simulation model aiming at each group of simulation parameters, and reading simulation section sizes of designated 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; the verified simulation model is determined by the following method: 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; comparing the initial simulation result with the actual ablation result to obtain a comparison result so as to determine the verified simulation model based on the comparison result; 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; 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; 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; based on the first sub-model, the second sub-model and the third sub-model, calculating to obtain microwave ablation parameters matched with the region to be ablated, wherein the microwave ablation parameters comprise: microwave ablation power and microwave ablation time.
6. 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-4.
7. 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-4.
CN202311606412.5A 2023-11-29 2023-11-29 Method and device for determining microwave ablation parameters and electronic equipment Active CN117297758B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311606412.5A CN117297758B (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
CN202311606412.5A CN117297758B (en) 2023-11-29 2023-11-29 Method and device for determining microwave ablation parameters and electronic equipment

Publications (2)

Publication Number Publication Date
CN117297758A CN117297758A (en) 2023-12-29
CN117297758B true CN117297758B (en) 2024-01-26

Family

ID=89250286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311606412.5A Active CN117297758B (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) CN117297758B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003099116A1 (en) * 2002-05-17 2003-12-04 Case Western Reserve University A method for delineating an ablation lesion
CN111027225A (en) * 2019-12-25 2020-04-17 南京亿高微波系统工程有限公司 Ablation simulation method based on data fitting
WO2020247885A1 (en) * 2019-06-06 2020-12-10 TriAgenics, Inc. Ablation probe systems
CN115546213A (en) * 2022-11-30 2022-12-30 海杰亚（北京）医疗器械有限公司 Cryoablation range prediction method and device, storage medium and electronic equipment
CN116196096A (en) * 2023-01-05 2023-06-02 上海介航机器人有限公司 Ablation needle deployment method, system, computer device and storage medium
CN116983085A (en) * 2023-09-25 2023-11-03 上海微创电生理医疗科技股份有限公司 Ablation stove size prediction system and ablation equipment

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220151688A1 (en) * 2008-04-29 2022-05-19 Virginia Tech Intellectual Properties Inc. System and method for estimating tissue heating of a target ablation zone for electrical-energy based therapies
US20230355301A1 (en) * 2022-05-03 2023-11-09 Medtronic, Inc. Ablation assistance using numeric modeling

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003099116A1 (en) * 2002-05-17 2003-12-04 Case Western Reserve University A method for delineating an ablation lesion
WO2020247885A1 (en) * 2019-06-06 2020-12-10 TriAgenics, Inc. Ablation probe systems
CN111027225A (en) * 2019-12-25 2020-04-17 南京亿高微波系统工程有限公司 Ablation simulation method based on data fitting
CN115546213A (en) * 2022-11-30 2022-12-30 海杰亚（北京）医疗器械有限公司 Cryoablation range prediction method and device, storage medium and electronic equipment
CN116196096A (en) * 2023-01-05 2023-06-02 上海介航机器人有限公司 Ablation needle deployment method, system, computer device and storage medium
CN116983085A (en) * 2023-09-25 2023-11-03 上海微创电生理医疗科技股份有限公司 Ablation stove size prediction system and ablation equipment

Also Published As

Publication number Publication date
CN117297758A (en) 2023-12-29

Similar Documents

Publication Publication Date Title
CN109498155B (en) Rapid planning system and method in ablation treatment operation
CN107278303B (en) Therapy planning system and method for generating optimal treatment plan
CN105825073B (en) A kind of online radiotherapy planning quality control system
Amabile et al. Tissue shrinkage in microwave ablation of liver: an ex vivo predictive model
Prakash et al. Considerations for theoretical modelling of thermal ablation with catheter-based ultrasonic sources: Implications for treatment planning, monitoring and control
Ye et al. The accuracy of computer image-guided template for mandibular angle ostectomy
US20160335413A1 (en) System and method for minimally invasive thermal ablation treatment planning
KR20130140067A (en) Systems and methods for planning image-guided interventional procedures
JP2016538944A (en) Efficient treatment planning trade-off analysis
Farina et al. Tissue shrinkage in microwave thermal ablation: comparison of three commercial devices
Rijnen et al. Quality and comfort in head and neck hyperthermia: A redesign according to clinical experience and simulation studies
Singh et al. Numerical study to establish relationship between coagulation volume and target tip temperature during temperature-controlled radiofrequency ablation
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
Sebaratnam et al. Direct cost-analysis of Mohs micrographic surgery and traditional excision for basal cell carcinoma at initial margin clearance
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
CN116936112B (en) Regression model training and ablation planning method, device, equipment and medium
CN113974820A (en) Simulated ablation method, device, storage medium and equipment based on residual fitting
JP2022506906A (en) A system that assists in the provision of template treatment parameters for ablation treatment
Ziegle et al. RF-ablation pattern shaping employing switching channels of dual bipolar needle electrodes: ex vivo results
Poni et al. Rapid SAR optimization for hyperthermic oncology: combining multi-goal optimization and time-multiplexed steering for hotspot suppression
CN114983559A (en) Data processing method, device, system, equipment and storage medium for radio frequency ablation
Kang et al. Contour and osteotomy of free fibula transplant using a ruler template
Kusters et al. Evaluation of two independent dose prediction methods to personalize the automated radiotherapy planning process for prostate cancer

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
GR01 Patent grant
GR01 Patent grant