CN113974820B

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

Info

Publication number: CN113974820B
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: 2024-09-13
Anticipated expiration: 2041-09-10
Also published as: CN113974820A

Abstract

The invention discloses a simulated ablation method, a simulated ablation device, a computer-readable storage medium and computer-readable equipment based on residual fitting, wherein the method comprises the following steps: simulating a corresponding target ablation region according to the target ablation tissue; performing temperature field simulation according to the target ablation region to obtain a temperature distribution value; 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 region; according to the first residual region, the second ablation needle is subjected to ablation simulation to obtain a second needle application parameter set, a second coverage region and a second residual region, and by applying the method, cancellation fusion can be simulated to realize accurate prediction of the needle application parameter set.

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, device, storage medium and equipment based on residual fitting.

Background

The ablation technology is generally divided into two major technologies, namely, thermal ablation and cold ablation, and the principle is that an ablation needle is inserted into target tissues, and the ablation needle is utilized to generate local physical high temperature or low temperature to disintegrate cell tissue structures, so that tissue cells are directly necrotized, and the aim of local ablation is achieved.

However, the cold ablation of the target tissue needs to accurately predict the needle application parameter set of the ablation needle according to experience, so that the 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 simulated ablation storage medium and simulated ablation equipment based on residual fitting, which can simulate ablation so as to realize accurate prediction of a needle application parameter set.

According to a first aspect of an embodiment of the present invention, there is provided a simulated ablation method based on residual fitting, the method comprising: simulating a corresponding target ablation region according to the target ablation tissue; performing temperature field simulation according to the target ablation region to obtain a temperature distribution value; 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 region; and performing ablation simulation on the second ablation needle according to the first residual region to obtain a second needle application parameter set, a second coverage area and a second residual 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 under the condition that the simulated ablation area meets the preset index, determining a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set.

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

According to an embodiment of the present invention, performing temperature field simulation according to the target ablation region to obtain a temperature distribution value includes: performing heat conduction simulation on the target ablation region according to a biological heat conduction equation, and determining a region simulation heat capacity, a region simulation heat conductivity and an ablation heat quantity value; and determining a temperature distribution value according to the area simulation heat capacity, the area simulation heat conductivity and the ablation heat quantity 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 needle application parameter set, 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 region; screening the simulated ablation region according to the target ablation region, and determining a region coverage value; determining a simulated ablation region with the largest 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; and integrating the target ablation region and the first coverage region, and determining the first residual region.

According to an embodiment of the present invention, after simulating the corresponding target ablation zone according to the target ablation tissue, the method further comprises: simulating target reserved tissues according to the target ablation region to obtain a target reserved region; screening the simulated ablation region according to the target reserved region, and determining a region punishment 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 Shi Zhen parameter sets include 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 coordinates 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 control instructions according to the first needle application parameter set, wherein the control instructions are used for instructing the first ablation needle to execute specific operations.

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; and determining a third coverage area according to the third ablation area and the fourth ablation area.

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

According to a second aspect of embodiments of the present invention, there is also provided a simulated ablation device based on residual fitting, 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 performing temperature field simulation according to the target ablation area to obtain a temperature distribution value; the ablation simulation module is used for performing ablation simulation on the first ablation needle according to the temperature distribution value and the target ablation area to obtain a first needle application parameter set, a first coverage area and a first residual area; the ablation simulation module is further used for performing ablation simulation on the second ablation needle according to the first residual region 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 according to 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 area 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 region to obtain other needle application parameter sets, other coverage areas, and other residual regions when the simulated ablation region does not meet a preset index; the determining module is further configured to redetermine a simulated ablation area 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 includes: performing heat conduction simulation on the target ablation region according to a biological heat conduction equation, and determining a region simulation heat capacity, a region simulation heat conductivity and an ablation heat quantity value; and determining a temperature distribution value according to the area simulation heat capacity, the area simulation heat conductivity and the ablation heat quantity value.

According to an embodiment of the present invention, the ablation simulation module includes: the simulation sub-module is used for 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 region; the screening submodule is used for screening the simulated ablation area according to the target ablation area and determining an area coverage value; the determining submodule is used for determining a simulated ablation area with the largest corresponding area coverage value as the first ablation area; the determining submodule is further used for determining the simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set; and the integration sub-module is used for integrating the target ablation area and the first coverage area and determining the first residual area.

According to an embodiment of the present invention, the area simulation module is further configured to simulate the target reserved tissue according to the target ablation area, so as to obtain a target reserved area; the apparatus further comprises: the screening module is used for screening the simulated ablation area according to the target reserved area and determining an area punishment 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 area simulation module is further configured to determine a first ablation area corresponding to a first ablation parameter, simulate the first ablation area according to the target ablation area, and obtain a coordinate point set corresponding to the first coverage area; the determining module is further configured to determine an ablation center coordinate corresponding to the first coverage area according to the coordinate point set.

According to an embodiment of the invention, the device further comprises: the generation module is used for generating control instructions according to the first needle application parameter set, and the control instructions are used for indicating the first ablation needle to execute specific operations.

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 area, so as to obtain a third needle application parameter set and a third ablation area; 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; and determining a third coverage area according to 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 mobile 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 area, and a fifth residual region; the positions of a first ablation center corresponding to the first needle application parameter set and a 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 also provided an apparatus comprising: one or more processors; and 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 residual fitting-based simulated ablation methods described above.

According to a fourth aspect of embodiments of the present invention, there is further provided a computer readable storage medium comprising a set of computer executable instructions for performing any of the above-described residual fitting based simulated ablation methods when executed.

According to the simulated ablation method, the simulated ablation device, the storage medium and the simulated ablation equipment based on residual fitting, provided by the embodiment of the invention, the temperature distribution value can be obtained by simulating the target ablation area corresponding to the target ablation tissue and performing temperature field simulation on the target ablation area; and performing ablation simulation on the ablation needle through the temperature distribution value, and determining simulation to obtain a corresponding needle application parameter set, a corresponding coverage area and a corresponding residual area. And so on, thereby realizing the purpose of simulating the ablation to accurately predict the needle application parameter set. By means of 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, the needle application mode of the ablation needle can be obtained in advance, reasonable arrangement and time setting of the ablation needle are assisted, and the purpose of accurately ablating the target ablation area is achieved.

It should be understood that the teachings of the present invention need not achieve all of the benefits set forth above, but rather that certain technical solutions may achieve certain technical effects, and that other embodiments of the present invention may also achieve benefits not set forth above.

Drawings

The above, as well as additional purposes, features, and advantages of exemplary embodiments of the present invention will become readily apparent from the following detailed description when read in conjunction with the accompanying drawings. Several embodiments of the present 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 shows a schematic implementation flow diagram of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

fig. 2 shows a second implementation flow chart of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

fig. 3 shows a third implementation flow diagram of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

fig. 4 shows a schematic implementation flow diagram of a simulated ablation method based on residual fitting according to an embodiment of the present invention;

FIG. 5 shows a schematic diagram of an implementation module of a simulated ablation device based on residual fitting in accordance with an embodiment of the present invention;

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

Detailed Description

The principles and spirit of the present invention will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are presented merely to enable those skilled in the art to better understand and practice the invention and are not intended to limit the scope of the 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 scheme of the invention is further elaborated below with reference to the drawings and specific embodiments.

Fig. 1 shows a schematic implementation flow diagram of a simulated ablation method based on residual fitting according to an embodiment of the present invention.

Referring to fig. 1, according to a first aspect of an embodiment of the present invention, there is provided a simulated ablation method based on residual fitting, the method comprising: operation 101, simulating a corresponding target ablation zone according to target ablation tissue; operation 102, performing temperature field simulation according to a target ablation region to obtain a temperature distribution value; an 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 needle application parameter set, a first coverage region and a first residual region; and operation 104, performing ablation simulation on the second ablation needle according to the first residual region to obtain a second needle application parameter set, a second coverage area and a second residual region.

According to the simulated ablation method based on residual fitting, provided by the embodiment of the invention, the target ablation area corresponding to the target ablation tissue is simulated, and the temperature distribution value can be obtained by simulating the temperature field of the target ablation area; the ablation simulation is carried out on the ablation needle through the temperature distribution value, a corresponding needle application parameter set, a corresponding coverage area and a corresponding residual area can be obtained through the determination simulation, and the like, so that the purpose of accurately predicting the needle application parameter set is achieved by simulating the ablation. 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 method operation 101, a target ablated tissue is tissue of a living being that needs to be ablated, such as tumor tissue or the like. 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 a target ablation region. Furthermore, the method can simulate the shape and the size of the target ablation tissue in a Cartesian coordinate system to obtain the target ablation region. The method can acquire the three-dimensional image of the target ablation tissue through the image acquisition device so as to obtain the three-dimensional image of the target ablation tissue, and the three-dimensional image is simulated so as to obtain the target ablation region. In one embodiment, a three-dimensional image of the target ablated tissue may be obtained by CT slicing and three-dimensional reconstruction of the target ablated tissue.

In operation 102 of the method, by performing temperature field simulation on the target ablation region, a temperature change rule of the target ablation region can be known, and a temperature distribution value corresponding to the target ablation region can be obtained. Specifically, the method carries out temperature field simulation based on the thermal transfer process of the ablation needle in the target ablation tissue along with time change so as to know the corresponding temperature distribution value under the condition of target ablation area ablation. Further, according to practical situations, when the ablation process of the present application is transferred in a heat conduction manner, the temperature distribution value can be simulated by a heat conduction equation, and a numerical solution method can be generally adopted. When the ablation process of the present application is delivered in a convective heat transfer regime, the temperature profile can be modeled by an energy profile if the velocity profile is known. Namely, according to the specific implementation scene of the target ablation region, the temperature distribution value can be simulated by selecting a corresponding temperature field simulation equation.

In the present method operation 103, with the known temperature distribution value and the target ablation region, the first ablation region corresponding to the first ablation needle first ablation may be determined by performing a simulation of a needle application parameter set for characterizing 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, etc. The first needle application parameter set is used to characterize parameter information of the first ablation needle application. The first ablation region is used to characterize a theoretical ablation region corresponding to the first needle application parameter set. The first coverage area is used to characterize the area covered by the first ablation area when the first ablation area is covered by the target simulation area. That is, when the first ablation zone is fully overlaid onto the target simulation zone, the first ablation zone is the same size as the first coverage zone, and when the first ablation zone is partially overlaid onto the target simulation zone, the first ablation zone is different size from the first coverage zone. Depending on the first needle application parameter set, the first coverage area may be the same size as the target ablation area, and the first coverage area may also be different size from the target ablation area. Specifically, the first coverage area may be sized beyond the target ablation area, and the first coverage area may also be sized smaller than the target ablation area. Typically, since the shape of the ablation zone with which the ablation is directed is fixed and the shape of the target ablation zone is not fixed, the method typically sets the size of the first coverage zone to not exceed the target ablation zone in order to avoid ablation to other remaining areas. The first residual region is used to characterize an area of the target ablation region that is not covered by the first ablation region.

In the method operation 104, in order to be able to sufficiently ablate the target ablation region, the method needs to perform a second ablation simulation on the first residual region to achieve ablation of the first residual region. In one embodiment, the method may redetermine the first residual region as a target ablation region, perform temperature field simulation on the target ablation region and perform ablation simulation on the second ablation needle, and implement ablation on the first residual region. In another embodiment, since the target ablation region corresponds to the same target ablation tissue, the temperature distribution value obtained by the temperature field simulation of the target ablation region is already known in operation 102, and the method can directly perform ablation simulation on the first residual region through the temperature distribution value obtained in operation 102, so as to realize ablation on the first residual region.

It should be added that, according to the actual situation of the target ablation area, when the method can complete ablation of the target ablation area once, that is, when the size of the first coverage area is equal to the target ablation area and the size of the first residual area 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 target ablation area and the second residual area is 0, the ablation is not required to be continued, otherwise, the ablation simulation … is performed on other ablation needles according to other residual areas, the ablation simulation of the fourth ablation needle and the ablation simulation of the fifth ablation needle can be performed, and details are omitted below until the 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 merely used as distinction between different ablation needles to facilitate understanding, and have no substantial distinction.

To facilitate a further understanding of the above embodiments, a specific implementation scenario is provided below for illustration. In this scenario, the method is applied to an ablation medical device that receives a CT slice from a CT machine, the CT slice containing an image corresponding to a target ablated tissue. Firstly, three-dimensional recombination is carried out on CT slices, and a three-dimensional image corresponding to target ablation tissue is obtained. Then, simulation is performed in a Cartesian coordinate system according to the three-dimensional image, and a target ablation region is obtained. And then, performing temperature field simulation according to 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 performing simulation adjustment on the needle application parameter set of the first ablation needle so that the size of the ablation region can be as close to or consistent with the size of the target ablation region as possible. Then, after the first ablation simulation is completed, an ablation region left after the first ablation needle is ablated is obtained. And simulating the next ablation needle according to the rest ablation area and the temperature distribution value according to the same method, and so on until the complete ablation of the target ablation area is completed.

Fig. 2 shows a second implementation flow chart of a simulated ablation method based on residual fitting according to an embodiment of the invention.

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

In the method operation 201, in one implementation scenario, the method may determine the simulated ablation region by direct accumulation 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, a simulated ablation area is determined according to an actual overlapping range of the first coverage area and the second coverage area in a coordinate system, and specifically, the simulated ablation area can be obtained by subtracting the actual overlapping range after accumulating the first coverage area and the second coverage area. The specific determination mode of the simulated ablation area can be adaptively adjusted according to the actual scene. After determining the simulated ablation area, the method determines subsequent operation by judging whether the simulated ablation area meets preset indexes.

In one implementation scenario, the preset index 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 yet another implementation scenario, the preset index 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-mentioned multiple preset indexes. The specific index content of the preset index can be adaptively adjusted according to the actual scene.

In operation 202 of the method, in the case that the simulated ablation region meets the preset index, it may be that the simulated ablation region completely covers the target ablation region, and the needle application simulation scheme is determined according to the needle application parameter set. Specifically, a required needle application parameter set is determined according to the simulation times, and a 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 needle application parameter set exists in the actual situation, determining a corresponding needle application simulation scheme according to the first needle application parameter set. If the third needle application parameter set exists, a corresponding needle application simulation scheme is required to be determined according to the first needle application parameter set, the second needle application parameter set and the third needle application parameter set. Hereinafter, details will not be described. In a subsequent operation, the medical device may output a corresponding protocol report according to the needle delivery simulation protocol to provide to the staff.

After operation 201, the method further comprises, according to an embodiment of the present invention: an operation 203, performing ablation simulation on other ablation needles according to the second residual error region under the condition that the simulated ablation region does not meet the preset index, so as to obtain other needle application parameter sets, other coverage areas and other residual error regions; operation 204 redetermines the simulated ablation region based on the first coverage area, the second coverage area, and the other coverage areas.

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

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

According to an embodiment of the present invention, operation 102, performing temperature field simulation according to a target ablation region to obtain a temperature distribution value, includes: firstly, performing heat conduction simulation on a target ablation region according to a biological heat conduction equation, and determining a region simulation heat capacity, a region simulation heat conductivity and an ablation heat quantity value; then, a temperature distribution value is determined from the area simulated heat capacity, the area simulated heat conductivity, and the ablation heat quantity 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. Thus, a numerical simulation of the temperature field of thermal transfer between biological tissues at different times for each ablation needle is obtained.

Specifically, the method adopts a biological heat conduction equation based on Pennes to simulate a temperature field, and the specific formula of the biological heat conduction equation is as follows:

wherein, C is used for characterizing 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 for representing the time corresponding to the ablation needle; k is used for representing the heat conductivity coefficient corresponding to the target ablation tissue; x is used to characterize each point on the target ablated tissue, which can be characterized in terms of 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 in time.

C _b is used to characterize the blood heat capacity corresponding to the target ablated tissue; omega _b is used to characterize the blood flow perfusion volume corresponding to the target ablated tissue; q _m is used to characterize the amount of heat generation of the effective metabolism corresponding to the target ablated tissue, and T _a is used to characterize the arterial temperature corresponding to the target ablated tissue.

In the simulation of the method, the biological heat conduction equation is optimized, specifically, the optimized mode may be omitted or replaced by a constant, and then the optimized heat conduction equation is:

Wherein, The heat required for ablation corresponding to each point of the target ablation zone can be characterized; The self heat corresponding to each point of the target ablation zone can be characterized, and Qr can be used for characterizing the heat required to be provided by the ablation needle corresponding to each point of the target ablation zone.

When the heat capacity C of the tissue, the heat conductivity k and the heat Q _r of the ablation needle are known, the temperature distribution value of the target ablation region, that is, the distribution of the temperature value corresponding to each coordinate point of the target ablation region can be obtained through analysis by the heat conduction equation.

Fig. 3 shows a schematic implementation flow diagram 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, operation 103 performs 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 region, including: 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 simulated data set comprises a simulated needle application parameter set and a corresponding simulated ablation region; operation 1032, screening the simulated ablation region according to the target ablation region, and determining a region coverage value; an operation 1033 of determining the simulated ablation region with the largest corresponding region coverage value as the first coverage region; an operation 1034 of determining a simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set; operation 1035, integrating the target ablation region with the first coverage region, determines a first residual region.

After the temperature distribution value is determined, the method can simulate the simulated ablation area corresponding to the ablation needle according to the temperature distribution value.

In particular, the set of needle delivery parameters of the ablation needle can be correlated to the ablation needle heat Q _r. Specifically, the heat Q _r of the ablation needle is related to the power W of the ablation needle, the size d of the ablation needle and the ablation time t, and 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 ellipsoid, the standard volume formula of the ellipsoid in the cartesian coordinate system can be obtained:

Wherein, (x, y, z) is used to represent the boundary coordinate point of the ellipsoid with the geometrical center of the ellipsoid at the origin, and abc is used to represent the radius of the ellipsoid.

Correspondingly, the simulated ablation area corresponding to the ablation needle is a set of all points in the ellipsoid, and the ablation center of the ablation needle, namely the vertex of the ablation needle, is the geometric center of the ellipsoid, so that the set of all points in the simulated ablation area can be characterized as:

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

Mapping by needle application parameter setBy the expansion of the ablation needle power W, the ablation needle size d and the ablation time tRemapping is performed.

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

The corresponding fifth-order polynomial f ₁,f₂,f₃ can be constructed by combining the constraint of the heat Q _r of the ablation needle on the power W of the ablation needle, the size d of the ablation needle and the ablation time t, thereby performing the following functionsRemapping 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

Fitting the fifth-order polynomial f ₁,f₂,f₃ according to the actual situation, for example, fitting the fifth-order polynomial f1, f2, f3 according to the coordinate point set of the actual ablation area of the ablation needle under different parameter settings in the actual situation, so as to obtain the corresponding parameters of the polynomial, so as to determine the target mapping function, specifically, the simulated ablation area of each ablation needle in the static state can be characterized as follows:

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

Wherein S is used to characterize the simulated ablation region corresponding to each ablation needle.

Based on the above formula, by using the ablation needle to reach the ablation temperature region under different ablation needle power W, ablation needle size d and ablation time t, specifically, the first ablation region corresponding to the first ablation needle may be characterized as S ₁, the second ablation region corresponding to the second ablation needle may be characterized as S ₂, and so on, which will not be described in detail below.

By means of preset stepping traversal of the ablation needle power W, the ablation needle size d, the ablation time t deflection angle theta 1 and the deflection angle theta 2, multiple groups of simulation needle application parameter sets and corresponding simulation ablation areas corresponding to the ablation needles can be obtained, and integration of the multiple groups of simulation needle application parameter sets and the corresponding simulation ablation areas is a simulation data set.

The simulated ablation regions in the simulated data set are screened through the target ablation region to determine the simulated ablation region, and in one implementation, the screening criteria may be to determine the simulated ablation region that is capable of maximizing coverage of the target ablation region. Wherein whether to maximize coverage of the target ablation zone may be manifested by a zone coverage value.

Can be specifically characterized as ：S₁＝f_a(W₁,d₁,t₁,θ₁₁,θ₂₁)＝argmax(S_A∩S₁)

Wherein S _A is used to characterize the target ablation region and S ₁ is used to characterize the simulated ablation region corresponding to the ablation needle. S _A∩S₁ is used to characterize the area coverage value, argmax (S _A∩S₁) is used to characterize the area coverage value with the largest value.

The method further 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; then, ablation center coordinates corresponding to the first coverage area are determined from the set of coordinate points.

Specifically, under the condition that the maximum area coverage value is determined by argmax (S _A∩S₁), the actual parameter corresponding to the ablation needle and the simulated ablation area S ₁ can be obtained, and further, the coordinate of the ablation center can be determined by calculating the coordinate point in the cartesian coordinate system corresponding to the geometric center of the simulated ablation area S ₁, wherein the coordinate point is the point of the temperature source generated by the ablation needle, and is usually the vertex of the ablation needle.

Specifically, the searching method of the geometric center (X, Y, Z) is that in the formed simulated ablation region S ₁ point set (X, Y, Z), the searching method is determined by the following formula:

x＝(Xmax-Xmin)

y＝(Ymax-Ymin)

z＝(Zmax-Zmin)

wherein Xmax is a maximum X-axis coordinate, and Xmin is a 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 the 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 _r1 is used to characterize the first residual region, it can be understood that if S _r1 =0, it is indicated that the first ablation needle can already achieve the purpose of ablating the target ablation tissue, and 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 region, S _r1 may be redetermined as the target ablation region, the second ablation parameters corresponding to the second ablation needle are obtained through the ablation simulation provided in operation 103, and so on, until S _r1 =0. It should be added that if S _r1 < 0, the solution may damage the target reserved area, and the solution is adopted or discarded according to the actual situation of the target reserved area.

Fig. 4 shows a schematic implementation flow diagram of a simulated ablation method based on residual fitting according to an embodiment of the invention.

Referring to fig. 4, after simulating a corresponding target ablation zone according to the target ablation tissue in operation 101, the method further includes: operation 401, simulating target reserved tissue according to the target ablation region to obtain a target reserved region; operation 402, screening the simulated ablation region according to the target reserved region, and determining a region punishment value; and operation 403, 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.

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 target ablation area. In the method, during simulation, the target ablation tissue is simulated firstly, and then the target reserved tissue is simulated, so that the relative position between the target ablation region and the target reserved region is determined, and a basis is provided for the purpose of reducing the target reserved region as much as possible while the target ablation region is ablated.

In one implementation scenario, assuming that the target ablated tissue is tumor tissue, the positive direction of the X-axis may be kept consistent with the front of the living body when a 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 that is to be ablated by the simulated ablation region, i.e., a region penalty value, and in particular, 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, the corresponding first coverage area is determined by taking into account the region penalty value and the region coverage value. Specifically, the region penalty value may be used to modify the region coverage value by means of difference calculation, so as to determine a simulated ablation region that maximally covers the target ablation region, 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 _H is used to characterize the target reserve area. (S _H∩S₁) for characterizing the region penalty value.

Furthermore, the tolerance simulation can be performed on the target reserved area to determine a tolerance value of the target reserved area, a punishment weight value of the target reserved area and the punishment value of the target reserved area can be determined according to the tolerance value, and the punishment value of the target reserved area is weighted through the punishment weight value, so that the punishment value of the target reserved area is more suitable for an actual scene. Wherein, the weight range can be (0, in particular, + -infinity ], if the tolerance value is higher, the penalty weight can be set to a larger value, if the tolerance value is low, the penalty weight may be set to a smaller value, such as any positive number less than 1.

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

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

Where λ 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 positions between the target ablation region and the target retention region are far enough that the target retention region is completely impossible to intersect with the simulated ablation region, the region penalty value is null, i.e.

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

When the method is applied to a medical device with an ablation function, after the first needle application parameter set is obtained according to the above embodiment, a corresponding control instruction may be generated according to the first needle application parameter set to instruct the first ablation needle to perform a specific operation, and 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 needle application parameter set.

According to an embodiment of the present invention, after obtaining the temperature distribution value in operation 102, the method further includes: firstly, performing cold ablation simulation on a third ablation needle according to a temperature distribution value and a 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 area and the fourth ablation area.

In one implementation scenario, since various shapes may exist in the target ablation area, by matching the cold ablation simulation and the hot ablation simulation through the third ablation needle and the fourth ablation needle, overlapping portions of the third ablation needle and the fourth ablation needle can be offset, so that a concave surface is formed on the third coverage area, the fourth coverage area is not limited to an ellipsoid shape, but can be formed into various shapes with concave surfaces, so that the target ablation area is more suitable for the purpose of ablation, and the ablation effect of the completely conformal ablation target ablation area can be achieved by adjusting and controlling the position relationship between the first ablation area and the second ablation area.

It should be understood that the third ablation needle and the fourth ablation needle are only used to facilitate the distinction in terms of expression, 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.

The high temperature generated by the fourth ablation needle can be neutralized with the low temperature generated by the third ablation needle through the intersection part of the third ablation area and the fourth ablation area to form a temperature field neutralization, so that an ice ball is not formed in the intersection ablation area corresponding to the third ablation needle, namely the ablation area with the temperature neutralization of the third ablation needle and the fourth ablation needle cannot ablate the target ablation area. The third coverage area is an area obtained by removing the temperature neutralization area of the third ablation area by the fourth ablation area, and the overlapping area of the area coverage and the target ablation area is the third coverage area. The third coverage area is used to characterize the actual ablation area of the third and fourth ablation needles in the target ablation area.

It should be appreciated that the third coverage area may be sized beyond the target ablation area and the third coverage area may be sized smaller than the target ablation area. The method generally sets the size of the third coverage area to not exceed the target ablation zone.

By simulating the temperature neutralization between the third ablation area and the fourth ablation area, as the ablation areas corresponding to the ablation needles are ellipsoidal, under the condition that the two ellipsoidal ablation areas are intersected, the overlapped part can be neutralized, so that the rest ablation areas can form a concave shape, namely a third coverage area with a concave surface, and the concave surface can better fit the boundary of the target ablation area, thereby achieving the purpose of conforming to the boundary of the ablation target ablation area.

According to an embodiment of the present invention, after obtaining the first residual 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 area and a fifth residual error region; the positions of a first ablation center corresponding to the first needle application parameter set and a fifth ablation center corresponding to the fifth needle application parameter set are different.

In an implementation scenario, the method can control the first ablation needle to perform movement ablation simulation, and can enable the first ablation needle to dynamically move in the target ablation area according to the first ablation center and the fifth ablation center so as to ablate the target ablation area at different positions, so that the area ablated by the first ablation needle is not limited to an ellipsoid, but can form various shapes, and the method is suitable for the purpose of ablation of the target ablation area. The method can realize multiple ablations through one ablation needle, and reduces the usage amount of the ablation needle.

Specifically, in order to sufficiently ablate the target ablation region, the method needs to perform ablation simulation again on the first residual region so as to ablate the first residual region. According to the method, the first residual error area can be redetermined as the target ablation area, and the temperature field simulation and the fifth ablation simulation aiming at the first ablation needle are carried out on the target ablation area, so that the ablation of the first residual error area is realized. In another embodiment, since the target ablation region corresponds to the same target ablation tissue, the temperature distribution value obtained by the temperature field simulation of the target ablation region is already known in operation 102, and the method can directly perform ablation simulation on the first residual region through the temperature distribution value obtained in operation 102, so as to realize ablation on the first residual region. Specifically, the method can realize the ablation of the first residual error area by controlling the first ablation needle to move from the first ablation center to the fifth ablation center and performing ablation at the fifth ablation center, namely the fifth ablation center is determined according to the position of the first residual error area. It is understood 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 tangential state, or a separated state. It should be appreciated that the present method may be applied to any ablation needle, not limited to the first ablation needle, depending on the circumstances.

Fig. 5 shows a schematic diagram of 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 also provided a simulated ablation device based on residual fitting, according to a second aspect of an embodiment of the invention, the device comprising: the region simulation module 501 is configured to simulate a corresponding target ablation region according to the target ablation tissue; the temperature field simulation module 502 is configured to perform temperature field simulation according to the target ablation area, and 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 area; the ablation simulation module 503 is further configured to perform ablation simulation on the second ablation needle according to the first residual region, so as to obtain a second needle application parameter set, a second coverage area, 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 area 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 region to obtain other needle application parameter sets, other coverage areas and other residual regions when the simulated ablation region does not meet a preset index; the determining module 504 is further configured to redefine a simulated ablation area 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 method comprises the steps of performing heat conduction simulation on the target ablation region according to a biological heat conduction equation, and determining a region simulation heat capacity, a region simulation heat conductivity and an ablation heat quantity value; and determining a temperature distribution value according to the area simulation heat capacity, the area simulation heat conductivity and the ablation heat quantity value.

According to an embodiment of the present invention, the ablation simulation module 503 includes: a simulation sub-module 5031, configured to perform parameter simulation on the first ablation needle according to the temperature distribution value, and 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 region; 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, as the first ablation region, a simulated ablation region with a corresponding region coverage value having a largest value; the determining submodule 5033 is further configured to determine a simulated needle application parameter set corresponding to the first coverage area as a first needle application parameter set; an integration sub-module 5034 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 area simulation module 501 is further configured to simulate the target reserved tissue according to the target ablation area, so as to obtain a target reserved area; the apparatus further comprises: the screening module 505 is configured to screen the simulated ablation area according to the target reserved area, and determine an area penalty value; and a correction 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 area simulation module 501 is further configured to determine a first ablation area corresponding to a first ablation parameter, simulate the first ablation area according to the target ablation area, and obtain a coordinate point set corresponding to the first coverage area; the determining module 504 is further configured to determine an ablation center coordinate corresponding to the first coverage area according to the coordinate point set.

According to an embodiment of the invention, the device further comprises: a generating module 507 is configured to generate a control instruction according to the first needle application parameter set, where the control instruction is configured 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; and determining a third coverage area according to 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 mobile 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 area, and a fifth residual region; the positions of a first ablation center corresponding to the first needle application parameter set and a fifth ablation center corresponding to the fifth needle application parameter set are different.

It should be noted here that: the above description of an embodiment of a simulated ablation device based on residual fitting, which is similar to the description of the method embodiment shown in fig. 1-5 and has similar advantageous effects as the method embodiment shown in fig. 1-5, is not repeated. For technical details not disclosed in the embodiment of the simulated ablation device based on residual fitting of the present invention, please refer to the description of the method embodiments of fig. 1 to 5 for the sake of economy, and therefore, details are not repeated.

Referring to fig. 6, there is also provided an apparatus according to a third aspect of the present invention, the apparatus comprising: one or more processors; and 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 residual fitting-based simulated ablation methods described above.

At the hardware level the device comprises a processor 601, optionally together with 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 (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 interconnected by an internal bus 603, which internal bus 603 may be an ISA (Industry Standard Architecture ) bus, a PCI (PERIPHERAL COMPONENT INTERCONNECT, peripheral component interconnect standard) bus, or an EISA (Extended Industry Standard Architecture ) bus, or the like. The buses may be divided into address buses, data buses, control buses, etc. For ease of illustration, only one double-headed arrow is shown in the figures, but not only one bus or one type of bus.

Memory 602 for storing execution instructions. In particular, a computer program that executes instructions may be executed. The memory 602 may include memory and non-volatile storage and provides the processor with instructions and data for execution.

In one possible implementation, the processor 601 reads the corresponding execution instruction from the nonvolatile memory into the memory and then executes the corresponding execution instruction, and may also acquire the corresponding execution instruction from other devices, so as to form a simulated ablation method based on residual fitting on a logic level. The processor executes the execution instructions stored in the memory to implement the simulated ablation method based on residual fitting provided in any embodiment of the invention by executing the execution instructions.

The method performed by the residual fitting-based simulated ablation method provided in the embodiment of the invention shown in fig. 6 can be applied to the processor 601 or implemented by the processor 601. The processor 601 may be an integrated circuit chip with signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in the processor 601 or instructions in the form of software. The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; but may also be a digital signal Processor (DIGITAL SIGNAL Processor, DSP), application SPECIFIC INTEGRATED Circuit (ASIC), field-Programmable gate array (Field-Programmable GATE ARRAY, FPGA) or other Programmable logic device, discrete gate or transistor logic device, 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 a memory, and the processor 601 reads the information in the memory 602, and in combination with its hardware, performs the steps of the method described above.

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 for performing any of the above-described simulated ablation methods based on residual fitting when executed.

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 one … …" 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 by 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 only illustrative, e.g. the division of the units is only one logical function division, and there may be other divisions in practice, such as: multiple units or components may be combined or may be integrated into another system, or some features may be omitted, or not performed. In addition, the various components shown or discussed may be coupled or directly coupled or communicatively coupled to each other via some interface, whether indirectly coupled or communicatively coupled to devices or units, whether electrically, mechanically, or otherwise.

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

In addition, each functional unit in each embodiment of the present invention may be integrated in one processing unit, or each unit may be separately used as one unit, or two or more units may be integrated in one unit; the integrated units may be implemented in hardware or in hardware plus software functional units.

Those of ordinary skill in the art will appreciate that: all or part of the steps for implementing the above method embodiments may be implemented by hardware related to program instructions, and the foregoing program may be stored in a computer readable storage medium, where the program, when executed, performs steps including the above method embodiments; and the aforementioned storage medium includes: a mobile storage device, a Read Only Memory (ROM), a magnetic disk or an optical disk, or the like, which can store program codes.

Or the above-described integrated units of the invention may be stored in a computer-readable storage medium if implemented in the form of software functional modules and sold or used as separate products. Based on such understanding, the technical solutions of the embodiments of the present invention may be embodied in essence or a part contributing to the prior art in the form of a software product stored in a storage medium, including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the methods described in the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a removable storage device, a ROM, a magnetic disk, or an optical disk.

The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

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

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

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

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 region;

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

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

when the simulated ablation area meets a preset index, determining a needle application simulation scheme according to the first needle application parameter set and the second needle application parameter set;

The first coverage area is used for representing an area covered by the first ablation area when the first ablation area is covered on the target ablation area; the first residual region is used for representing a region which is not covered by the first ablation region in the target ablation region; the first ablation region is used for representing a theoretical ablation region corresponding to the first needle application parameter set; the Shi Zhen parameter set is used to characterize parameter information associated with the ablation needle and capable of affecting the ablation region.

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

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

And redefining a simulated ablation area according to the first coverage area, the second coverage area and the other coverage areas.

3. The method of claim 1, wherein performing a temperature field simulation from the target ablation zone to obtain temperature distribution values comprises:

Performing heat conduction simulation on the target ablation region according to a biological heat conduction equation, and determining a region simulation heat capacity, a region simulation heat conductivity and an ablation heat quantity value;

And determining a temperature distribution value according to the area simulation heat capacity, the area simulation heat conductivity and the ablation heat quantity value.

4. The method of claim 1, wherein performing an ablation simulation on the first ablation needle based on the temperature distribution value and the target ablation region to obtain a first needle placement parameter set, 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 region;

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

determining a simulated ablation region with the largest 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;

and integrating the target ablation region and the first coverage region, and determining the first residual region.

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

Simulating target reserved tissues according to the target ablation region to obtain a target reserved region;

screening the simulated ablation region according to the target reserved region, and determining a region punishment 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.

6. The method of claim 1, wherein the Shi Zhen set of parameters includes ablation needle power, ablation needle size, ablation time, and ablation needle angle.

7. The method according to claim 1, wherein 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 coordinates corresponding to the first coverage area according to the coordinate point set.

8. The method according to claim 1, wherein the method further comprises:

And generating control instructions according to the first needle application parameter set, wherein the control instructions are used for instructing the first ablation needle to execute specific operations.

9. The method according to claim 1, wherein 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;

and determining a third coverage area according to the third ablation area and the fourth ablation area.

10. 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 region to obtain a fifth needle application parameter set, a fifth coverage region and a fifth residual region;

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

11. A simulated ablation device based on residual fitting, 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 performing temperature field simulation according to the target ablation area to obtain a temperature distribution value;

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

the ablation simulation module is further used for performing ablation simulation on the second ablation needle according to the temperature distribution value and the first residual region to obtain a second needle application parameter set, a second coverage region and a second residual region;

a determining module for determining a simulated ablation region according to 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 area meets a preset index;

12. An apparatus, the apparatus comprising:

One or more processors;

storage means for storing one or more programs,

When executed by the one or more processors, causes the one or more processors to implement the method of any of claims 1-10.

13. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the method according to any one of claims 1-10.