CN113974820A

CN113974820A - Simulated ablation method, device, storage medium and equipment based on residual fitting

Info

Publication number: CN113974820A
Application number: CN202111064065.9A
Authority: CN
Inventors: 周翔; 焦晓康
Original assignee: Yidu Cloud Beijing Technology Co Ltd
Current assignee: Yidu Cloud Beijing Technology Co Ltd
Priority date: 2021-09-10
Filing date: 2021-09-10
Publication date: 2022-01-28
Anticipated expiration: 2041-09-10
Also published as: CN113974820B

Abstract

The invention discloses a simulated ablation method, a simulated ablation device, a computer-readable storage medium and equipment based on residual fitting, wherein the method comprises the following steps: simulating a corresponding target ablation region according to the target ablation tissue; carrying out temperature field simulation according to the target ablation region to obtain a temperature distribution value; performing ablation simulation on a first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region; and performing ablation simulation on the second ablation needle according to the first residual error region to obtain a second needle application parameter set, a second coverage region and a second residual error region.

Description

Simulated ablation method, device, storage medium and equipment based on residual fitting

Technical Field

The invention relates to the technical field of medical instruments, in particular to a simulated ablation method, a simulated ablation device, a storage medium and equipment based on residual error fitting.

Background

The ablation technology is divided into two main technologies of heat ablation and cold ablation, and the principle is that an ablation needle is inserted into target tissue, and local physical high temperature or low temperature is generated by the ablation needle to disintegrate a cell tissue structure to cause direct necrosis of tissue cells, so that the purpose of local ablation is achieved.

However, cold ablation of target tissue requires accurate prediction of the needle application parameter set of the ablation needle according to experience, so that accurate cryoablation of the target tissue is realized.

Disclosure of Invention

The embodiment of the invention provides a simulated ablation method, a simulated ablation device, a storage medium and equipment based on residual error fitting, which can be used for simulating ablation so as to realize accurate prediction of an injection parameter set.

According to a first aspect of the embodiments of the present invention, there is provided a simulated ablation method based on residual fitting, the method including: simulating a corresponding target ablation region according to the target ablation tissue; carrying out temperature field simulation according to the target ablation region to obtain a temperature distribution value; performing ablation simulation on a first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region; and performing ablation simulation on the second ablation needle according to the first residual error region to obtain a second needle application parameter set, a second coverage region and a second residual error region.

According to an embodiment of the invention, after obtaining the second coverage area, the method further comprises: determining a simulated ablation region from the first coverage region and the second coverage region; and determining a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set under the condition that the simulated ablation region meets a preset index.

According to an embodiment of the invention, after determining the simulated ablation zone, the method further comprises: when the simulated ablation region does not meet the preset index, performing ablation simulation on other ablation needles according to the second residual error region to obtain other needle application parameter sets, other coverage regions and other residual error regions; re-determining a simulated ablation zone from the first coverage zone, the second coverage zone, and the other coverage zones.

According to an embodiment of the present invention, performing a temperature field simulation according to the target ablation region to obtain a temperature distribution value includes: performing heat conduction simulation on the target ablation area according to a biological heat conduction equation, and determining area simulation heat capacity, area simulation heat conductivity and ablation heat value; and determining a temperature distribution value according to the region simulated heat capacity, the region simulated heat conductivity and the ablation heat value.

According to an embodiment of the present invention, the performing ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region to obtain a first set of needle application parameters, a first coverage area, and a first residual area includes: performing parameter simulation on the first ablation needle according to the temperature distribution value to obtain a simulation data set corresponding to the first ablation needle; wherein the simulated data set comprises a simulated needle application parameter set and a corresponding simulated ablation zone; screening the simulated ablation region according to the target ablation region, and determining a region coverage value; determining the simulated ablation region with the maximum corresponding region coverage value as the first ablation region; determining a simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set; determining the first residual region for integration of the target ablation region and the first coverage region.

According to an embodiment of the invention, after simulating the corresponding target ablation zone according to the target ablation tissue, the method further comprises: simulating a target reserved tissue according to the target ablation region to obtain a target reserved region; screening the simulated ablation region according to the target retention region, and determining a region penalty value; and correcting the area coverage value according to the area penalty value to obtain a corrected coverage value, wherein the corrected coverage value is used for determining the first coverage area.

According to an embodiment of the present invention, the set of needle parameters includes ablation needle power, ablation needle size, ablation time, and ablation needle angle.

According to an embodiment of the invention, the method further comprises: determining a first ablation region corresponding to a first ablation parameter, and simulating the first ablation region according to the target ablation region to obtain a coordinate point set corresponding to the first coverage region; and determining the ablation center coordinate corresponding to the first coverage area according to the coordinate point set.

According to an embodiment of the invention, the method further comprises: and generating a control instruction according to the first needle application parameter set, wherein the control instruction is used for instructing the first ablation needle to execute a specific operation.

According to an embodiment of the invention, after obtaining the temperature distribution value, the method further comprises: performing cold ablation simulation on a third ablation needle according to the temperature distribution value and the target ablation region to obtain a third needle application parameter set and a third ablation region; performing thermal ablation simulation on a fourth ablation needle according to the third ablation region and the target ablation region to obtain a fourth needle application parameter set and a fourth ablation region; determining a third coverage area from the third ablation area and the fourth ablation area.

According to an embodiment of the present invention, after obtaining the first residual region, the method further includes: performing mobile ablation simulation on the first ablation needle according to the first residual error region to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual error region; wherein the positions of the first ablation center corresponding to the first needle application parameter set and the fifth ablation center corresponding to the fifth needle application parameter set are different.

There is further provided, in accordance with a second aspect of the embodiments of the present invention, a simulated ablation apparatus based on residual fitting, the apparatus including: the region simulation module is used for simulating a corresponding target ablation region according to the target ablation tissue; the temperature field simulation module is used for carrying out temperature field simulation according to the target ablation region to obtain a temperature distribution value; the ablation simulation module is used for carrying out ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region; the ablation simulation module is further configured to perform ablation simulation on a second ablation needle according to the first residual region, so as to obtain a second needle application parameter set, a second coverage region and a second residual region.

According to an embodiment of the present invention, the apparatus further comprises a determining module for determining a simulated ablation region from the first coverage region and the second coverage region; the determining module is further configured to determine a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set when the simulated ablation region meets a preset index.

According to an embodiment of the present invention, the determining module is further configured to perform ablation simulation on other ablation needles according to the second residual error region to obtain other needle application parameter sets, other coverage regions, and other residual error regions when the simulated ablation region does not meet a preset index; the determining module is further configured to re-determine a simulated ablation region based on the first coverage area, the second coverage area, and the other coverage areas.

According to an embodiment of the present invention, the temperature field simulation module includes: performing heat conduction simulation on the target ablation area according to a biological heat conduction equation, and determining area simulation heat capacity, area simulation heat conductivity and ablation heat value; and determining a temperature distribution value according to the region simulated heat capacity, the region simulated heat conductivity and the ablation heat value.

According to an embodiment of the invention, the ablation simulation module comprises: the simulation submodule is used for carrying out parameter simulation on the first ablation needle according to the temperature distribution value to obtain a simulation data set corresponding to the first ablation needle; wherein the simulated data set comprises a simulated needle application parameter set and a corresponding simulated ablation zone; the screening submodule is used for screening the simulated ablation region according to the target ablation region and determining a region coverage value; the determining submodule is used for determining the simulated ablation region with the maximum corresponding region coverage value as the first ablation region; the determining submodule is further configured to determine an analog injection parameter set corresponding to the first coverage area as a first injection parameter set; an integration sub-module for integrating the target ablation region and the first coverage region to determine the first residual region.

According to an embodiment of the present invention, the region simulation module is further configured to simulate a target retaining tissue according to the target ablation region to obtain a target retaining region; the device further comprises: the screening module is used for screening the simulated ablation region according to the target reserved region and determining a region penalty value; and the correction module is used for correcting the area coverage value according to the area penalty value to obtain a corrected coverage value, and the corrected coverage value is used for determining the first coverage area.

According to an embodiment of the present invention, the region simulation module is further configured to determine a first ablation region corresponding to a first ablation parameter, and simulate the first ablation region according to the target ablation region to obtain a coordinate point set corresponding to the first coverage region; the determining module is further configured to determine ablation center coordinates corresponding to the first coverage area according to the set of coordinate points.

According to an embodiment of the invention, the apparatus further comprises: and the generating module is used for generating a control instruction according to the first needle application parameter set, wherein the control instruction is used for instructing the first ablation needle to execute a specific operation.

According to an embodiment of the present invention, the ablation simulation module is further configured to perform cold ablation simulation on a third ablation needle according to the temperature distribution value and the target ablation region, so as to obtain a third needle application parameter set and a third ablation region; performing thermal ablation simulation on a fourth ablation needle according to the third ablation region and the target ablation region to obtain a fourth needle application parameter set and a fourth ablation region; determining a third coverage area from the third ablation area and the fourth ablation area.

According to an embodiment of the present invention, the ablation simulation module is further configured to perform a moving ablation simulation on the first ablation needle according to the first residual region, so as to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual region; wherein the positions of the first ablation center corresponding to the first needle application parameter set and the fifth ablation center corresponding to the fifth needle application parameter set are different.

According to a third aspect of embodiments of the present invention, there is further provided an apparatus, including: one or more processors; a storage device for storing one or more programs that, when executed by the one or more processors, cause the one or more processors to implement any of the above-described methods for simulated ablation based on residual fit.

There is yet further provided, in accordance with the fourth aspect of the embodiments of the present invention, a computer-readable storage medium, comprising a set of computer-executable instructions, which when executed, perform any of the above-described methods for simulated ablation based on residual fit.

The simulated ablation method, the simulated ablation device, the storage medium and the equipment based on residual fitting provided by the embodiment of the invention can obtain a temperature distribution value by simulating a target ablation region corresponding to a target ablation tissue and simulating a temperature field of the target ablation region; and performing ablation simulation on the ablation needle through the temperature distribution value, and determining a corresponding needle application parameter set, a coverage area and a residual error area obtained through simulation. By analogy, the ablation is simulated so as to realize the purpose of accurately predicting the injection parameter set. By means of the method for simulating the target ablation area based on residual fitting, the method can achieve the needle application simulation scheme of the target ablation area through residual simulation prediction, can obtain the needle application mode of the ablation needle in advance, assists in reasonable arrangement and time setting of the ablation needle, and achieves the purpose of accurately ablating the target ablation area.

It is to be understood that the teachings of the present invention need not achieve all of the above-described benefits, but rather that specific embodiments may achieve specific technical results, and that other embodiments of the present invention may achieve benefits not mentioned above.

Drawings

The above and other objects, features and advantages of exemplary embodiments of the present invention will become readily apparent from the following detailed description read in conjunction with the accompanying drawings. Several embodiments of the invention are illustrated by way of example, and not by way of limitation, in the figures of the accompanying drawings and in which:

in the drawings, the same or corresponding reference numerals indicate the same or corresponding parts.

FIG. 1 is a schematic diagram showing a first implementation flow of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

FIG. 2 is a schematic diagram showing a second implementation flow of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

FIG. 3 is a schematic diagram showing a third implementation flow of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

FIG. 4 is a schematic diagram illustrating a fourth implementation flow of the simulated ablation method based on residual fitting according to the embodiment of the present invention;

FIG. 5 is a schematic diagram of an implementation module of a simulated ablation device based on residual fit according to an embodiment of the invention;

fig. 6 shows an implementation structure diagram of a device according to an embodiment of the present invention.

Detailed Description

The principles and spirit of the present invention will be described with reference to a number of exemplary embodiments. It is understood that these embodiments are given only to enable those skilled in the art to better understand and to implement the present invention, and do not limit the scope of the present invention in any way. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

The technical solution of the present invention is further elaborated below with reference to the drawings and the specific embodiments.

Fig. 1 is a schematic flow chart showing an implementation process of a simulated ablation method based on residual error fitting according to an embodiment of the present invention.

Referring to fig. 1, according to a first aspect of the embodiments of the present invention, there is provided a method of simulated ablation based on residual fitting, the method including: an operation 101 of simulating a corresponding target ablation zone according to a target ablation tissue; operation 102, performing temperature field simulation according to the target ablation region to obtain a temperature distribution value; operation 103, performing ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region; and operation 104, performing ablation simulation on the second ablation needle according to the first residual error region, and obtaining a second needle application parameter set, a second coverage region and a second residual error region.

According to the simulated ablation method based on residual fitting, provided by the embodiment of the invention, the temperature distribution value can be obtained by simulating the target ablation region corresponding to the target ablation tissue and simulating the temperature field of the target ablation region; the ablation simulation is carried out on the ablation needle through the temperature distribution value, the corresponding needle application parameter set, the corresponding coverage area and the corresponding residual error area can be obtained through the simulation, and the like, so that the purpose of simulating the ablation to accurately predict the needle application parameter set is realized. The method is suitable for medical equipment with an ablation function, and can also be suitable for a control device or an auxiliary device for controlling the medical equipment.

In the present method operation 101, the target ablation tissue is tissue of an organism that needs to be ablated, such as tumor tissue. The method can simulate the target ablation tissue according to at least one of the shape and the size of the target ablation tissue to obtain the target ablation region. Further, the method can simulate the shape and the size of the target ablation tissue in a Cartesian coordinate system to obtain a target ablation area. It is added that the method can acquire a three-dimensional image of the target ablation tissue through the image acquisition device to obtain a three-dimensional image of the target ablation tissue, and acquire the target ablation region through the simulation of the three-dimensional image. In one embodiment, a three-dimensional image of the target ablated tissue can be obtained by CT sectioning and three-dimensional reconstruction of the target ablated tissue.

In the method operation 102, a temperature field simulation is performed on the target ablation region, so that a temperature change rule of the target ablation region and a temperature distribution value corresponding to the target ablation region can be known. Specifically, the method performs temperature field simulation based on the thermal transfer process of the ablation needle in the target ablation tissue along with the time change, so as to know the corresponding temperature distribution value under the ablation condition of the target ablation region. Further, as the ablation process of the present application is conducted in a thermal conduction manner, the temperature distribution values can be modeled by a thermal conduction equation, and a numerical solution method can be generally used. When the ablation process of the present application is delivered in convective heat transfer, the temperature profile can be modeled by an energy equation if the velocity profile is known. According to the specific implementation scene of the target ablation region, the corresponding temperature field simulation equation can be selected to simulate the temperature distribution value.

In the method operation 103, with the known temperature distribution values and target ablation regions, a first ablation region corresponding to the first ablation needle for the first time may be determined by performing a simulation of a set of needle application parameters on the first ablation needle, wherein the set of needle application parameters is used to characterize parameter information related to the ablation needle and capable of affecting the ablation region, including but not limited to ablation needle power, ablation needle size, ablation time, ablation needle angle, and the like. The first needle application parameter set is used for characterizing the parameter information of the first ablation needle application. The first ablation region is used to characterize a theoretical ablation region corresponding to a first set of needle parameters. The first coverage area is used to characterize the area covered by the first ablation region when the first ablation region is overlaid on the target simulation area. That is, the first ablation zone is the same size as the first coverage zone when the first ablation zone completely overlies the target simulation zone, and the first ablation zone is a different size than the first coverage zone when the first ablation zone partially overlies the target simulation zone. Depending on the first set of needle parameters, the first coverage area may be the same size as the target ablation area, or the first coverage area may be a different size than the target ablation area. In particular, the first coverage area may be sized to exceed the target ablation area, and the first coverage area may be sized to be smaller than the target ablation area. Typically, since the shape of the ablation region to which the ablation needle is directed is fixed and the shape of the target ablation region is not fixed, the method typically sets the first coverage area to a size that does not exceed the target ablation region in order to avoid ablation to other remaining areas. The first residual region is used to characterize a region of the target ablation region not covered by the first ablation region.

In operation 104 of the method, in order to sufficiently ablate the target ablation region, the method needs to perform a second ablation simulation on the first residual region to perform ablation on the first residual region. In one possible implementation, the method can determine the first residual region as the target ablation region again, perform temperature field simulation on the target ablation region and ablation simulation on the second ablation needle, and realize ablation on the first residual region. In another possible implementation, since the target ablation region corresponds to the same target ablation tissue, and the temperature distribution values obtained by the temperature field simulation are already known in operation 102, the method can perform ablation simulation on the first residual region directly through the temperature distribution values obtained in operation 102, so as to achieve ablation of the first residual region.

It should be added that, according to the actual situation of the target ablation region, when the method can complete the complete ablation of the target ablation region by one ablation, that is, under the condition that the size of the first coverage region is equal to the target ablation region and the size of the first residual region is 0, operation 104 is not required. Similarly, when the sizes of the first coverage area and the second coverage area are equal to the size of the target ablation area and the size of the second residual error area is 0, ablation is not required to be continued, otherwise, ablation simulation … is performed on other ablation needles according to other residual error areas, and so on, and ablation simulation of the fourth ablation needle and ablation simulation of the fifth ablation needle can also be performed, which is not described in detail below until ablation of the target ablation area is completed. It should be added that the expressions "first" and "second" in the first ablation needle and the second ablation needle in the method are only used as a distinction between different ablation needles for convenience of understanding, and there is no substantial difference.

To facilitate further understanding of the above embodiments, a specific implementation scenario is provided below for description. In this scenario, the method is applied to an ablation medical device that receives CT slices from a CT machine, the CT slices containing images corresponding to the target ablated tissue. Firstly, the CT slices are three-dimensionally recombined to obtain a three-dimensional image corresponding to the target ablation tissue. Then, simulation is carried out in a Cartesian coordinate system according to the three-dimensional image, and a target ablation area is obtained. And then, carrying out temperature field simulation according to the tissue parameters corresponding to the target ablation tissue to obtain a temperature distribution value. And then, performing ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region, and performing simulation adjustment on the ablation region corresponding to the first ablation needle by simulating and adjusting the needle application parameter set of the first ablation needle so as to enable the size of the ablation region to be as close to or consistent with the size of the target ablation region as possible. Then, after completing the first ablation simulation, obtaining the ablation region left after the ablation of the first ablation needle. And simulating the next ablation needle according to the rest ablation regions and temperature distribution values by the same method, and repeating the steps until the complete ablation of the target ablation region is completed.

Fig. 2 is a schematic diagram illustrating an implementation flow of a simulated ablation method based on residual error fitting according to an embodiment of the present invention.

Referring to fig. 2, after obtaining the second coverage area, the method further includes, in accordance with an embodiment of the present invention: an operation 201, determining a simulated ablation region according to a first coverage region and a second coverage region; in operation 202, a needle application simulation plan is determined according to the first set of needle application parameters and the second set of needle application parameters when the simulated ablation region meets the preset index.

In method operation 201, in one implementation scenario, the method may determine the simulated ablation region by direct summation of the first coverage region and the second coverage region. In another implementation scenario, after the first coverage area and the second coverage area are simulated, the simulated ablation area is determined according to the actual overlapping range of the first coverage area and the second coverage area in the coordinate system, and specifically, the simulated ablation area is obtained by subtracting the actual overlapping range after the first coverage area and the second coverage area are accumulated. The specific determination mode of the simulated ablation region of the method can be adaptively adjusted according to the actual scene. After the simulated ablation zone is determined, the method determines the follow-up operation by judging whether the simulated ablation zone meets the preset index.

In one implementation scenario, the preset indicator may be that the simulated ablation region completely covers the target ablation region. In another implementation scenario, the preset index may be that the area of the residual region is 0. In another implementation scenario, the preset indicator may be the number of needle applications. In other implementation scenarios, the preset index of the method may also be a combination of the above multiple preset indexes. The specific index content of the preset index can be adaptively adjusted according to the actual scene.

In the method operation 202, when the simulated ablation region meets the preset criterion, that is, the simulated ablation region completely covers the target ablation region, the needle application simulation scheme is determined according to the needle application parameter set. Specifically, the required needle application parameter set is determined according to the simulation times, and the corresponding needle application simulation scheme is determined according to the first needle application parameter set and the second needle application parameter set. If only the first needling parameter set exists in the practical situation, determining a corresponding needling simulation scheme according to the first needling parameter set. If the third needling parameter set exists, the corresponding needling simulation scheme is determined according to the first needling parameter set, the second needling parameter set and the third needling parameter set. And will not be described in detail below. In subsequent operations, the medical device may output a corresponding protocol report to be provided to the staff in accordance with the needle application simulation protocol.

After operation 201, the method further includes, according to an embodiment of the invention: operation 203, when the simulated ablation region does not meet the preset index, performing ablation simulation on other ablation needles according to the second residual error region to obtain other needle application parameter sets, other coverage regions and other residual error regions; at operation 204, the simulated ablation zone is re-determined based on the first coverage zone, the second coverage zone, and the other coverage zones.

In operation 203 of the method, when the simulated ablation region does not meet the preset index, that is, the simulated ablation region does not completely cover the target ablation region, the method needs to continue to simulate other ablation needle parameter sets, so as to achieve ablation of the second residual region. Furthermore, the method can simulate at least two ablation needles at one time to determine corresponding needle application parameter sets, coverage areas and residual error areas; it is also possible to simulate only one ablation needle at a time to determine the corresponding set of needle parameters, coverage area and residual area.

In the method operation 204, the simulated ablation region is re-determined by integrating the coverage region corresponding to the ablation needle, so as to determine the preset index. It is to be added that the method of determining the simulated ablation zone in operation 204 is the same as in operation 201. It should be further added that the setting of operation 202 and operation 203 of the method is to distinguish each operation step, and there is no sequential relationship between operation 202 and operation 203.

According to an embodiment of the present invention, in operation 102, a temperature field simulation is performed according to a target ablation region to obtain a temperature distribution value, including: firstly, carrying out heat conduction simulation on a target ablation area according to a biological heat conduction equation, and determining area simulation heat capacity, area simulation heat conductivity and ablation heat value; then, a temperature distribution value is determined according to the region simulated heat capacity, the region simulated heat conductivity and the ablation heat value.

The method simulates a biological heat transfer process based on the change of an ablation needle along with time based on a biological heat conduction equation. Thereby obtaining the numerical simulation of the temperature field of each ablation needle in the heat transfer between biological tissues at different times.

Specifically, the method adopts a Pennes-based biological heat conduction equation to carry out temperature field simulation, and the biological heat conduction equation has the following specific formula:

wherein C is used for representing the tissue heat capacity corresponding to the target ablation tissue; t is used to characterize the tissue temperature corresponding to the target ablated tissue; t is used to characterize the time corresponding to the ablation needle; k is used to characterize the thermal conductivity corresponding to the target ablated tissue; x is used to characterize each point on the target ablated tissue, which can be characterized in coordinates (X, y, z) by modeling the target ablated tissue in a cartesian coordinate system; t (X, T) may be used to characterize the temperature at each point and at each instant in time.

C_bFor characterizing a blood heat capacity corresponding to the target ablated tissue; omega_bFor characterizing a perfusion volume of blood flow corresponding to the target ablated tissue; q_mHeat generation, T, for characterizing effective metabolism corresponding to target ablated tissue_aFor characterizing arterial temperature corresponding to the target ablated tissue.

Wherein the blood-related parameter C_d、ω_b、T_aAnd a parameter Q related to metabolism_mIt is difficult to obtain immediately, so in the simulation of the method, the biological heat conduction equation is optimized, specifically, the optimization mode may be omitted or replaced by a constant, and then the optimized heat conduction equation is:

wherein,

the heat required for ablation of each point in the target ablation region can be represented;

the self-heat corresponding to each point of the target ablation region can be characterized, and the Qr can be used for representing the heat required to be provided by the ablation needle corresponding to each point of the target ablation region.

When the tissue heat capacity C, the heat conductivity k and the heat Q of the ablation needle_rUnder the known condition, the temperature distribution value of the target ablation area, namely the distribution condition of the temperature value corresponding to each coordinate point of the target ablation area, can be obtained through analysis by a heat conduction equation.

Fig. 3 is a schematic flow chart showing an implementation process of a simulated ablation method based on residual fitting according to an embodiment of the present invention.

Referring to fig. 3, according to an embodiment of the present invention, in operation 103, performing ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region, and obtaining a first set of needle application parameters, a first coverage area and a first residual area, includes: operation 1031, performing parameter simulation on the first ablation needle according to the temperature distribution value to obtain a simulation data set corresponding to the first ablation needle; wherein the simulation data set comprises a simulation needle application parameter set and a corresponding simulation ablation region; operation 1032, screening the simulated ablation region according to the target ablation region, and determining a region coverage value; at operation 1033, determining the simulated ablation zone with the largest corresponding zone coverage value as the first coverage zone; operation 1034, determine a simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set; at operation 1035, a first residual region is determined for the target ablation region and the first coverage region integration.

After determining the temperature distribution value, the method can simulate the simulated ablation region corresponding to the ablation needle according to the temperature distribution value.

Specifically, the parameter set of the ablation needle and the heat Q of the ablation needle_rAn association may be made. In particular, the heat Q of the ablation needle_rAssociated with the power W of the ablation needle, the size d of the ablation needle and the ablation time t, since the shape of the temperature field corresponding to each ablation of the ablation needle is known, specifically, the shape of the temperature field corresponding to the ablation needle is an ellipsoid, a standard volume formula of the ellipsoid in a cartesian coordinate system can be obtained:

wherein, (x, y, z) is used for representing the ellipsoid boundary coordinate point of the ellipsoid geometric center at the origin, and abc is used for representing the radius of the ellipsoid.

Correspondingly, the simulated ablation region corresponding to the ablation needle is a set of all points in an ellipsoid, and the ablation center of the ablation needle, i.e. the vertex of the ablation needle, is the geometric center of the ellipsoid, so that it can be obtained that the set of all points of the simulated ablation region can be characterized as follows:

as above, where (x, y, z) is used to characterize all coordinate points of the simulated ablation region at the origin of the geometric center of the simulated ablation region, abc is used to characterize the radius of the simulated ablation region.

Mapping to by administering a set of parameters

The expansion is carried out by the ablation needle power W, the ablation needle size d and the ablation time t

A remapping is performed.

It is further contemplated that when the ablation needle is also simulated in a cartesian coordinate system, the ablation needle angle may also affect the shape of the simulated ablation zone and may be characterized by an angle θ 1 of the ablation needle away from the positive direction of the x-axis and an angle θ 2 of the ablation needle away from the positive direction of the z-axis.

Combined with the heat Q of the ablation needle_rThe corresponding fifth-order polynomial f can be constructed according to the constraints of the power W of the ablation needle, the size d of the ablation needle and the ablation time t₁，f₂，f₃Thereby to pair

A remapping is performed.

The specific fifth order polynomial is as follows:

x²×f₁(W，d，t，θ₁)+y²×f₂(W，d，t，θ₁)+z²×f₃(W，d，t，θ₂)≤1

according to the actual situation, the above-mentioned fifth order polynomial f₁，f₂，f₃Fitting is performed, such as fitting the fifth-order polynomials f1, f2, f3 according to the coordinate point sets of the actual ablation regions of the ablation needles under different parameter settings in actual cases, so as to obtain the corresponding parameters of the polynomials, and determine the target mapping function, specifically, each ablation needle is inThe simulated ablation zone in the resting state can be characterized as:

S＝f_a(W，d，t，θ₁，θ₂)

wherein S is used to characterize the simulated ablation zone corresponding to each ablation needle.

Based on the above formula, by setting the ablation needle to reach the region of the ablation temperature at different ablation needle powers W, ablation needle sizes d and ablation times t, specifically, the first ablation region corresponding to the first ablation needle can be characterized as S₁The second ablation zone for the second ablation needle may be characterized as S₂So on, it is not described in detail below.

Multiple groups of simulated needle application parameter sets corresponding to the ablation needles and corresponding simulated ablation regions can be obtained by presetting step traversal ablation needle power W, ablation needle size d, ablation time tposition angle theta 1 and declination angle theta 2, and the integration of the multiple groups of simulated needle application parameter sets and the corresponding simulated ablation regions is a simulated data set.

The simulated ablation zone in the simulated data set is screened through the target ablation zone to determine a simulated ablation zone, and in one embodiment, the screening criteria may be to determine a simulated ablation zone that maximizes coverage of the target ablation zone. Wherein, whether the coverage of the target ablation area is maximized can be reflected by the area coverage value.

The method can be characterized in particular as follows: s₁＝f_a(W₁，d₁，t₁，θ₁₁，θ₂₁)＝argmax(S_A∩S₁)

Wherein S is_AFor characterizing a target ablation zone, S₁For characterizing a simulated ablation zone corresponding to the ablation needle. S_A∩S₁For characterizing the area coverage value, argmax (S)_A∩S₁) For characterizing the area coverage value with the largest value.

The method also comprises the following steps: firstly, simulating a first coverage area according to a target ablation area to obtain a coordinate point set corresponding to the first coverage area; ablation center coordinates corresponding to the first coverage area are then determined from the set of coordinate points.

Specifically, by argmax (S)_A∩S₁) Under the condition of determining the maximum area coverage value, the actual parameters corresponding to the ablation needle and the simulated ablation area S can be obtained₁Further, the ablation region S is simulated by calculation₁The geometrical center of (a) corresponds to a coordinate point in a cartesian coordinate system, i.e. the ablation center coordinate can be determined, i.e. the point where the ablation needle generates the temperature source, typically the apex of the ablation needle.

Specifically, the geometric center (x, y, z) is found in the formed simulated ablation region S₁In the point set (X, Y, Z), it is determined by the following formula:

x＝(Xmax-Xmin)

y＝(Ymax-Ymin)

z＝(Zmax-Zmin)

wherein Xmax is the maximum X-axis coordinate, and Xmin is the minimum X-axis coordinate; ymax-is the maximum Y-axis coordinate, and Ymin is the minimum Y-axis coordinate; zmax is the maximum Z-axis coordinate and Zmin is the minimum Z-axis coordinate.

The first residual region is a difference between the target ablation region and the first coverage region, and can be characterized by the following formula:

S_r1＝S_A-f_a(W₁，d₁，t₁，θ₁₁，θ₁₂)，

wherein S is_r1For characterizing the first residual region, it can be understood that if S is_r1When the first ablation needle reaches the target tissue ablation, the first needle application parameter set corresponding to the first ablation needle can be directly output. If S_r1> 0, there is a remaining target ablation zone, S can be adjusted_r1Re-determining the target ablation region, obtaining a second ablation parameter corresponding to the second ablation needle through the ablation simulation provided in operation 103, and so on until S_r10. It is necessary to supplement if S_r1< 0, the scheme may damage the target reserved area, and the scheme is adopted or rejected according to the actual situation of the target reserved areaAnd (5) abandoning.

Fig. 4 is a schematic flow chart illustrating an implementation process of a simulated ablation method based on residual error fitting according to an embodiment of the present invention.

Referring to fig. 4, after simulating a corresponding target ablation zone according to target ablation tissue in operation 101, according to an embodiment of the present invention, the method further includes: operation 401, simulating a target reserved tissue according to a target ablation region to obtain a target reserved region; operation 402, screening the simulated ablation region according to the target retention region, and determining a region penalty value; in operation 403, the area coverage value is modified according to the area penalty value, so as to obtain a modified coverage value, where the modified coverage value is used to determine the first coverage area.

In the operation of operation 101, in order to further improve the accuracy of the ablation procedure, it is necessary to ensure that the target remaining area is not ablated in addition to the ablation target ablation area. During simulation, the method firstly simulates the target ablation tissue and then simulates the target retaining tissue to determine the relative position between the target ablation region and the target retaining region, thereby providing a basis for reducing the target retaining region as much as possible while ablating the target ablation region.

In one implementation scenario, assuming the target ablation tissue is tumor tissue, the positive direction of the X-axis can be kept consistent with the front of the living body when the cartesian coordinate system is established. The relative position between the target ablation zone and the target retention zone can be determined by the coordinate system.

In the method operation 402, in the case of obtaining the simulated ablation region, the simulated ablation region may be overlapped with the target retention region to determine a target retention region range to be ablated by the simulated ablation region, that is, a region penalty value, and specifically, the region penalty value may be characterized by an intersection of the target retention region and the simulated ablation region.

In the method operation 403, a corresponding first coverage area is determined by comprehensively considering the area penalty value and the area coverage value. Specifically, the area penalty value may be corrected by a difference calculation manner to determine a simulated ablation area that maximally covers the target ablation area, and specifically, the following formula may be used for characterization:

S₁＝f_a(W₁，d₁，t₁，θ₁₁，θ₂₁)＝argmax((S_A∩S₁)-(S_H∩S₁))

wherein S is_HFor characterizing the target reserved area. (S)_H∩S₁) For characterizing the regional penalty values.

Furthermore, the method can determine a tolerance value corresponding to the target retention area by performing tolerance simulation on the target retention area, determine a penalty weight corresponding to the area penalty value according to the tolerance value, and weight the area penalty value through the penalty weight, so that the area penalty value is more suitable for an actual scene. The value of the penalty weight may be set to be a larger value, such as any positive number greater than 1, if the tolerance value is higher, or to be a smaller value, such as any positive number less than 1, if the tolerance value is lower.

Specifically, after adding the penalty weight to the formula, the method can be characterized as follows:

S₁＝f_a(W₁，d₁，t₁，θ₁₁，θ₂₁)＝argmax((S_A∩S₁)-λ×(S_H∩S₁))

wherein λ is used to characterize the penalty weight.

Further, in the case of simulating the target ablation region and the target retention region in the cartesian coordinate system, if the relative position between the target ablation region and the target retention region is relatively far, so that the target retention region is completely impossible to intersect with the simulated ablation region, the region penalty value is an empty set, that is, the region penalty value is an empty set

After operation 103, the method further includes, according to an embodiment of the invention: and generating a control instruction according to the first needle application parameter set, wherein the control instruction is used for instructing the first ablation needle to perform a specific operation.

When the method is applied to a medical device with an ablation function, after the first set of needle application parameters is obtained according to the above embodiment, a corresponding control instruction may be generated according to the first set of needle application parameters to instruct the first ablation needle to perform a specific operation, specifically, the specific operation may be to control the first ablation needle to move to a specified geometric center, and perform an ablation operation according to the first set of needle application parameters.

After obtaining the temperature distribution values at operation 102, the method further includes: firstly, performing cold ablation simulation on a third ablation needle according to the temperature distribution value and the target ablation region to obtain a third needle application parameter set and a third ablation region; then, performing thermal ablation simulation on a fourth ablation needle according to the third ablation region and the target ablation region to obtain a fourth needle application parameter set and a fourth ablation region; a third coverage area is determined from the third ablation region and the fourth ablation region.

In an implementation scenario, as various shapes of the target ablation region may exist, through the cooperation of the cold ablation simulation and the heat ablation simulation of the third ablation needle and the fourth ablation needle, the overlapped parts of the third ablation needle and the fourth ablation needle can be mutually offset, so that a concave surface is formed on the third coverage region, and the fourth coverage region is not limited to be an ellipsoid, but can be formed into various shapes with concave surfaces, so as to be more suitable for the purpose of ablation of the target ablation region, and by regulating and controlling the position relationship between the first ablation region and the second ablation region, the ablation effect of completely conformal ablation of the target ablation region can be realized.

It is to be understood that the third ablation needle and the fourth ablation needle in this embodiment are only used for convenience of description, and in practical cases, the third ablation needle and the fourth ablation needle may be the same as or different from the first ablation needle and the second ablation needle.

Through the intersection part of the third ablation region and the fourth ablation region, the high temperature generated by the fourth ablation needle can be neutralized with the low temperature generated by the third ablation needle to form temperature field neutralization, so that the intersection ablation region corresponding to the third ablation needle cannot form ice balls, namely, the ablation region neutralized by the temperatures of the third ablation needle and the fourth ablation needle cannot ablate the target ablation region. The third coverage area is obtained by removing the temperature neutralization area of the fourth ablation area from the third ablation area, and the overlapping area of the third coverage area and the target ablation area is the third coverage area. The third coverage area is used to characterize the actual ablation zone of the third and fourth ablation needles in the target ablation zone.

It should be added that the third coverage area may be sized to exceed the target ablation area, and the third coverage area may be sized to be smaller than the target ablation area. The method generally sets the third coverage area to a size not exceeding the target ablation area.

By simulating the temperature neutralization between the third ablation region and the fourth ablation region, as the ablation regions corresponding to the ablation needle are all in an ellipsoidal shape, under the condition that the two ellipsoidal ablation regions are intersected, the overlapped part can be neutralized, so that the rest ablation regions can form a concave shape, namely a third covering region with a concave surface, and the concave surface can be better attached to the boundary of the target ablation region, thereby achieving the purpose of conforming to the boundary of the target ablation region.

According to an embodiment of the present invention, after obtaining the first residual error region in operation 103, the method further includes: performing mobile ablation simulation on the first ablation needle according to the first residual error region to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual error region; wherein the positions of the first ablation center corresponding to the first needle application parameter set and the fifth ablation center corresponding to the fifth needle application parameter set are different.

In one implementation scenario, the method can control the first ablation needle to perform moving ablation simulation, and enable the first ablation needle to dynamically move in the target ablation region according to the first ablation center and the fifth ablation center, so as to ablate the target ablation region at different positions, so that the region ablated by the first ablation needle is not limited to an ellipsoid shape, but can be formed into various shapes, thereby being suitable for the purpose of ablating the target ablation region. The method can realize multiple ablation through one ablation needle, and reduces the use amount of the ablation needle.

Specifically, in order to sufficiently ablate the target ablation region, the method needs to perform another ablation simulation on the first residual region to perform ablation on the first residual region. The method can determine the first residual area as the target ablation area again, and carries out temperature field simulation and fifth ablation simulation aiming at the first ablation needle on the target ablation area to realize ablation of the first residual area. In another possible implementation, since the target ablation region corresponds to the same target ablation tissue, and the temperature distribution values obtained by the temperature field simulation are already known in operation 102, the method can perform ablation simulation on the first residual region directly through the temperature distribution values obtained in operation 102, so as to achieve ablation of the first residual region. Specifically, the method can control the first ablation needle to move from the first ablation center to the fifth ablation center, and the ablation is performed at the fifth ablation center to realize the ablation of the first residual region, namely the fifth ablation center is determined according to the position of the first residual region. It will be appreciated that depending on the relative positions of the first and fifth ablation centers, the first ablation region corresponding to the first ablation center and the fifth ablation region corresponding to the fifth ablation center may include, but are not limited to, a partially overlapping state, a tangent state, or a separated state. It is to be understood that the present method may be applied to any ablation needle, not limited to the first ablation needle, depending on the actual situation.

Fig. 5 is a schematic diagram illustrating an implementation module of a simulated ablation device based on residual fitting according to an embodiment of the invention.

Referring to fig. 5, there is further provided, in accordance with a second aspect of the embodiments of the present invention, a simulated ablation apparatus based on residual fitting, the apparatus including: a region simulation module 501, configured to simulate a corresponding target ablation region according to a target ablation tissue; a temperature field simulation module 502, configured to perform temperature field simulation according to the target ablation region to obtain a temperature distribution value; an ablation simulation module 503, configured to perform ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region, so as to obtain a first needle application parameter set, a first coverage area, and a first residual error area; the ablation simulation module 503 is further configured to perform ablation simulation on a second ablation needle according to the first residual region, so as to obtain a second needle application parameter set, a second coverage region, and a second residual region.

According to an embodiment of the present invention, the apparatus further comprises a determining module 504 for determining a simulated ablation region according to the first coverage region and the second coverage region; the determining module 504 is further configured to determine a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set when the simulated ablation region meets a preset index.

According to an embodiment of the present invention, the determining module 504 is further configured to perform ablation simulation on other ablation needles according to the second residual error region to obtain other needle application parameter sets, other coverage regions, and other residual error regions when the simulated ablation region does not meet a preset index; the determining module 504 is further configured to re-determine a simulated ablation region according to the first coverage area, the second coverage area, and the other coverage areas.

According to an embodiment of the present invention, the temperature field simulation module 502: the device is used for carrying out heat conduction simulation on the target ablation area according to a biological heat conduction equation, and determining area simulation heat capacity, area simulation heat conductivity and ablation heat value; and determining a temperature distribution value according to the region simulated heat capacity, the region simulated heat conductivity and the ablation heat value.

According to an embodiment of the present invention, the ablation simulation module 503 includes: the simulation submodule 5031 is configured to perform parameter simulation on the first ablation needle according to the temperature distribution value to obtain a simulation data set corresponding to the first ablation needle; wherein the simulated data set comprises a simulated needle application parameter set and a corresponding simulated ablation zone; a screening submodule 5032, configured to screen the simulated ablation region according to the target ablation region, and determine a region coverage value; a determining submodule 5033, configured to determine the simulated ablation region with the largest corresponding region coverage value as the first ablation region; the determining sub-module 5033 is further configured to determine an analog administering parameter set corresponding to the first coverage area as a first administering parameter set; an integration sub-module 5034 configured to integrate the target ablation region and the first coverage region and determine the first residual region.

According to an embodiment of the present invention, the region simulation module 501 is further configured to simulate a target reserved tissue according to the target ablation region, so as to obtain a target reserved region; the device further comprises: a screening module 505, configured to screen the simulated ablation region according to the target retention region, and determine a region penalty value; a correcting module 506, configured to correct the area coverage value according to the area penalty value, to obtain a corrected coverage value, where the corrected coverage value is used to determine the first coverage area.

According to an embodiment of the present invention, the region simulation module 501 is further configured to determine a first ablation region corresponding to a first ablation parameter, and simulate the first ablation region according to the target ablation region to obtain a coordinate point set corresponding to the first coverage region; the determining module 504 is further configured to determine ablation center coordinates corresponding to the first coverage area according to the set of coordinate points.

According to an embodiment of the invention, the apparatus further comprises: a generating module 507, configured to generate a control instruction according to the first needle application parameter set, where the control instruction is used to instruct the first ablation needle to perform a specific operation.

According to an embodiment of the present invention, the ablation simulation module 503 is further configured to perform cold ablation simulation on a third ablation needle according to the temperature distribution value and the target ablation region, so as to obtain a third needle application parameter set and a third ablation region; performing thermal ablation simulation on a fourth ablation needle according to the third ablation region and the target ablation region to obtain a fourth needle application parameter set and a fourth ablation region; determining a third coverage area from the third ablation area and the fourth ablation area.

According to an embodiment of the present invention, the ablation simulation module 503 is further configured to perform a moving ablation simulation on the first ablation needle according to the first residual region, so as to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual region; wherein the positions of the first ablation center corresponding to the first needle application parameter set and the fifth ablation center corresponding to the fifth needle application parameter set are different.

Here, it should be noted that: the above description of an embodiment of a simulated ablation apparatus based on residual fitting is similar to the description of the embodiment of the method shown in fig. 1 to 5, and has similar beneficial effects to the embodiment of the method shown in fig. 1 to 5, and therefore, the description thereof is omitted. For technical details not disclosed in the embodiment of the simulated ablation device based on residual error fitting, please refer to the description of the embodiment of the method shown in fig. 1 to 5 of the present invention, which will be omitted for brevity.

Referring to fig. 6, according to a third aspect of the present invention, there is also provided an apparatus comprising: one or more processors; a storage device for storing one or more programs which, when executed by the one or more processors, cause the one or more processors to implement any of the above-described methods of simulated ablation based on residual fit.

On the hardware level, the apparatus comprises a processor 601, optionally an internal bus 603, a network interface 604, a memory 602. The Memory 602 may include a Memory, such as a Random-Access Memory (RAM), and may further include a non-volatile Memory, such as at least 1 disk Memory. Of course, the device may also include hardware required for other services.

The processor 601, the network interface 604, and the memory 602 may be connected to each other by an internal bus 603, and the internal bus 603 may be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, an EISA (Extended Industry Standard Architecture) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown, but this does not indicate only one bus or one type of bus.

A memory 602 for storing instructions for execution. In particular, a computer program that can be executed by executing instructions. The memory 602 may include both memory and non-volatile storage and provides execution instructions and data to the processor.

In a possible implementation manner, the processor 601 reads the corresponding execution instruction from the non-volatile memory into the memory and then runs the corresponding execution instruction, and may also obtain the corresponding execution instruction from other devices, so as to form the simulated ablation method based on residual error fitting on a logic level. The processor executes the execution instructions stored in the memory to implement the simulated ablation method based on residual error fitting provided in any embodiment of the invention through the executed execution instructions.

The method performed by the simulated ablation method based on residual fitting provided in the embodiment of fig. 6 according to the present invention can be applied to the processor 601 or implemented by the processor 601. The processor 601 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 601. The Processor may be a general-purpose Processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but also Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps and logic blocks disclosed 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 directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory, and the processor 601 reads the information in the memory 602, and performs the steps of the above method in combination with hardware thereof.

According to a fourth aspect of the present invention, there is also provided a computer-readable storage medium comprising a set of computer-executable instructions that, when executed, perform any of the above-described methods of simulated ablation based on residual fit.

It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.

In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above-described device embodiments are merely illustrative, for example, the division of the unit is only a logical functional division, and there may be other division ways in actual implementation, such as: multiple units or components may be combined, or may be integrated into another system, or some features may be omitted, or not implemented. In addition, the coupling, direct coupling or communication connection between the components shown or discussed may be through some interfaces, and the indirect coupling or communication connection between the devices or units may be electrical, mechanical or other forms.

The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units; can be located in one place or distributed on a plurality of network units; some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.

In addition, all the functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may be separately regarded as one unit, or two or more units may be integrated into one unit; the integrated unit can be realized in a form of hardware, or in a form of hardware plus a software functional unit.

Those of ordinary skill in the art will understand that: all or part of the steps for realizing the method embodiments can be completed by hardware related to program instructions, the program can be stored in a computer readable storage medium, and the program executes the steps comprising the method embodiments when executed; and the aforementioned storage medium includes: various media that can store program codes, such as a removable Memory device, a Read Only Memory (ROM), a magnetic disk, or an optical disk.

Alternatively, the integrated unit of the present invention may be stored in a computer-readable storage medium if it is implemented in the form of a software functional module and sold or used as a separate product. Based on such understanding, the technical solutions of the embodiments of the present invention may be essentially implemented or a part contributing to the prior art may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the methods described in the embodiments of the present invention. And the aforementioned storage medium includes: a removable storage device, a ROM, a magnetic or optical disk, or other various media that can store program code.

The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims

1. A method of simulated ablation based on residual fit, the method comprising:

simulating a corresponding target ablation region according to the target ablation tissue;

carrying out temperature field simulation according to the target ablation region to obtain a temperature distribution value;

performing ablation simulation on a first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region;

and performing ablation simulation on the second ablation needle according to the first residual error region to obtain a second needle application parameter set, a second coverage region and a second residual error region.

2. The method of claim 1, wherein after obtaining the second coverage area, the method further comprises:

determining a simulated ablation region from the first coverage region and the second coverage region;

and determining a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set under the condition that the simulated ablation region meets a preset index.

3. The method of claim 2, wherein after determining the simulated ablation zone, the method further comprises:

when the simulated ablation region does not meet the preset index, performing ablation simulation on other ablation needles according to the second residual error region to obtain other needle application parameter sets, other coverage regions and other residual error regions;

re-determining a simulated ablation zone from the first coverage zone, the second coverage zone, and the other coverage zones.

4. The method of claim 1, wherein performing a temperature field simulation based on the target ablation region to obtain temperature distribution values comprises:

performing heat conduction simulation on the target ablation area according to a biological heat conduction equation, and determining area simulation heat capacity, area simulation heat conductivity and ablation heat value;

and determining a temperature distribution value according to the region simulated heat capacity, the region simulated heat conductivity and the ablation heat value.

5. The method of claim 1, wherein performing ablation simulation on the first ablation needle according to the temperature distribution values and the target ablation region to obtain a first set of needle application parameters, a first coverage region and a first residual region comprises:

performing parameter simulation on the first ablation needle according to the temperature distribution value to obtain a simulation data set corresponding to the first ablation needle; wherein the simulated data set comprises a simulated needle application parameter set and a corresponding simulated ablation zone;

screening the simulated ablation region according to the target ablation region, and determining a region coverage value;

determining the simulated ablation region with the maximum corresponding region coverage value as the first ablation region;

determining a simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set;

determining the first residual region for integration of the target ablation region and the first coverage region.

6. The method of claim 5, wherein after simulating the corresponding target ablation zone from the target ablation tissue, the method further comprises:

simulating a target reserved tissue according to the target ablation region to obtain a target reserved region;

screening the simulated ablation region according to the target retention region, and determining a region penalty value;

and correcting the area coverage value according to the area penalty value to obtain a corrected coverage value, wherein the corrected coverage value is used for determining the first coverage area.

7. The method of claim 1, wherein the set of needle parameters comprises ablation needle power, ablation needle size, ablation time, and ablation needle angle.

8. The method of claim 1, further comprising:

determining a first ablation region corresponding to a first ablation parameter, and simulating the first ablation region according to the target ablation region to obtain a coordinate point set corresponding to the first coverage region;

and determining the ablation center coordinate corresponding to the first coverage area according to the coordinate point set.

9. The method of claim 1, further comprising:

and generating a control instruction according to the first needle application parameter set, wherein the control instruction is used for instructing the first ablation needle to execute a specific operation.

10. The method of claim 1, wherein after obtaining the temperature distribution values, the method further comprises:

performing cold ablation simulation on a third ablation needle according to the temperature distribution value and the target ablation region to obtain a third needle application parameter set and a third ablation region;

performing thermal ablation simulation on a fourth ablation needle according to the third ablation region and the target ablation region to obtain a fourth needle application parameter set and a fourth ablation region;

determining a third coverage area from the third ablation area and the fourth ablation area.

11. The method of claim 1, wherein after obtaining the first residual region, the method further comprises:

performing mobile ablation simulation on the first ablation needle according to the first residual error region to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual error region;

wherein the positions of the first ablation center corresponding to the first needle application parameter set and the fifth ablation center corresponding to the fifth needle application parameter set are different.

12. A simulated ablation device based on residual fit, the device comprising:

the region simulation module is used for simulating a corresponding target ablation region according to the target ablation tissue;

the temperature field simulation module is used for carrying out temperature field simulation according to the target ablation region to obtain a temperature distribution value;

the ablation simulation module is used for carrying out ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation region to obtain a first needle application parameter set, a first coverage region and a first residual error region;

the ablation simulation module is further configured to perform ablation simulation on a second ablation needle according to the first residual region, so as to obtain a second needle application parameter set, a second coverage region and a second residual region.

13. An apparatus, characterized in that the apparatus comprises:

one or more processors;

a storage device for storing one or more programs,

when executed by the one or more processors, cause the one or more processors to implement the method of any one of claims 1-11.

14. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the method according to any one of claims 1-11.