Simulation method, simulation device and storage medium for releasing internal flow disturbance device of aneurysm

Abstract

The present application discloses a simulation method, simulation device and storage medium for releasing an internal spoiler device for an aneurysm. The method comprises: constructing a spring mass system for the internal spoiler device according to the selection of the internal spoiler device for the aneurysm, wherein the spring mass system comprises a plurality of mass points connected by a plurality of spring edges; setting a target constraint for the spring mass system in simulating the release of the internal spoiler device into the aneurysm; and optimizing the position of the target mass point in the spring mass system based on the target constraint to simulate the release of the internal spoiler device. Through the scheme of the present application, the process of releasing the internal spoiler device can be simulated efficiently and accurately, so as to improve the success rate of selecting and releasing the internal spoiler device during surgery.

Description

Technical Field

[0001] The present application generally relates to the field of simulation technology, and more specifically, to a simulation method, device and computer-readable storage medium for releasing an internal flow disturbance device of an aneurysm. Background Art

[0002] With the advancement of medical technology, intracranial lesions can be identified through quantitative analysis of medical images. Using the brain data of intracranial lesions for the preparation of preoperative instruments and equipment and intraoperative process planning can effectively improve the quality of surgery, shorten the duration of surgery, and reduce potential risks during surgery.

[0003] Vascular intervention surgery has become the mainstream surgical method for the clinical treatment of vascular diseases such as aneurysms and vascular occlusions. Compared with traditional craniotomy, patients have better prognosis and fewer complications. In recent years, intratumoral flow disturbance devices ("WEB") have achieved remarkable results in the treatment of aneurysms at the apex and bifurcation. In current vascular intervention surgery, clinical doctors generally select and release interventional consumables by simply manually measuring blood vessels or aneurysms before surgery, and rely on experience. However, as a new type of consumable, WEB has different evaluations of consumable selection and post-operative release effects by doctors with different experiences. Failure to select and inappropriately release WEB will not only prolong the operation time, but also increase the risk of poor prognosis, causing patients to bear unnecessary additional economic burdens, and even worsen the quality of surgery, affecting the patient's prognosis and causing recurrence. The existing WEB release simulation method is to reconstruct blood vessels and aneurysms from images, and use Finite Element Analysis (Finite Element Analysis, "FEA") to simulate the WEB release process. Although this method is slightly more accurate, it takes a long time to model and build, requiring at least a few hours, and at most dozens or even hundreds of hours, making it difficult to apply in actual clinical practice.

[0004] In view of this, there is an urgent need to provide a simulation method for releasing an internal spoiler device for an aneurysm, so as to efficiently and accurately simulate the process of releasing the internal spoiler device and improve the success rate of selecting and releasing the internal spoiler device during surgery. Summary of the invention

[0005] In order to at least solve one or more of the above-mentioned technical problems, the present application proposes a simulation solution for releasing an internal flow disrupting device for an aneurysm in multiple aspects.

[0006] In a first aspect, the present application provides a simulation method for the release of an internal spoiler device for an aneurysm, comprising: constructing a spring mass system for the internal spoiler device according to the selection of the internal spoiler device of the aneurysm, wherein the spring mass system comprises a plurality of mass points connected by a plurality of spring edges; in simulating the release of the internal spoiler device into the aneurysm, setting a target constraint for the spring mass system; and optimizing the position of the target mass point in the spring mass system based on the target constraint to simulate the release of the internal spoiler device.

[0007] In some embodiments, the multiple spring edges and the multiple mass points are connected to form a triangular network including multiple triangles, the target constraints include stretching constraints, bending constraints and internal pressure constraints, and setting the target constraints for the spring mass system includes: setting the stretching constraints on the edges of the triangles in the spring mass system; setting the bending constraints on adjacent triangles in the spring mass system; and setting the internal pressure constraints on the triangular network in the spring mass system.

[0008] In some other embodiments, the stretch constraint is set by the following formula: stretch (p1, p2) = |p1-p2|-l0, where C stretch (P1, P2) represents the stretching constraint, P1 and P2 represent the two vertices of the side of the triangle connecting the spring-mass system, and l0 represents the length of the initial side.

[0009] In some further embodiments, the bending constraint is set by the following formula:

[0010] Among them, C bend (p1, p2, p3, p4) represents the bending constraint, p1, p2, p3, p4 represent the vertices of two adjacent triangles, Indicates the initial angle.

[0011] In some further embodiments, the internal pressure constraint is set by the following formula:

[0012] Among them, C(P1,…,P N ) represents the internal pressure constraint, n triangles represents the number of triangles in the triangulated network, and represents the vertices in each triangle, k pressure represents the pressure factor and V0 represents the initial volume.

[0013] In some further embodiments, the method further includes setting a unidirectional constraint on a target vertex of a triangular network in the spring-mass system.

[0014] In some further embodiments, the unidirectional constraint includes a first unidirectional constraint and a second unidirectional constraint, and setting the unidirectional constraint on the target vertex of the triangular network in the spring mass system includes: in response to the target vertex of the triangular network in the spring mass system passing through the corresponding triangle, setting the first unidirectional constraint on the target vertex of the triangular network in the spring mass system; or in response to the target vertex of the triangular network in the spring mass system merging from the bottom of the corresponding triangle, setting the second unidirectional constraint on the target vertex of the triangular network in the spring mass system.

[0015] In some further embodiments, the first unidirectional constraint is set by the following formula:

[0016] Among them, C(q, P1, P2, P3) represents the first unidirectional constraint, q represents the target vertex, P1, P2, P3 represent the vertices of the corresponding triangles, and h represents the thickness of the triangular network in the spring-mass system.

[0017] In some further embodiments, the second unidirectional constraint is set by the following formula:

[0018] Among them, C(q, P1, P2, P3)' represents the second unidirectional constraint, q represents the target vertex, P1, P2, P3 represent the vertices of the corresponding triangle, and h represents the thickness of the triangular network in the spring-mass system.

[0019] In some other embodiments, the selection of the internal spoiler device of the aneurysm is obtained by the following operations: reconstructing the aneurysm surface and extracting the aneurysm neck; obtaining geometric parameters of a target geometric shape close to the aneurysm morphology based on fitting the aneurysm surface and the aneurysm neck; and calculating the size of the aneurysm based on the geometric parameters to obtain the selection of the internal spoiler device of the aneurysm.

[0020] In some other embodiments, the target geometric shape includes an ellipsoid, the geometric parameters include major axis parameters of the ellipsoid, and the size of the aneurysm includes the height and width of the aneurysm.

[0021] In some other embodiments, the method further includes: setting a fitting constraint; and obtaining the geometric parameters of a target geometric shape close to the aneurysm morphology based on fitting the aneurysm surface and the aneurysm neck under the fitting constraint.

[0022] In a second aspect, the present application provides a simulation device for releasing an internal spoiler device for an aneurysm, comprising: a processor; and a memory on which computer instructions for releasing an internal spoiler device for an aneurysm are stored, and when the computer instructions are executed by the processor, one or more embodiments of the aforementioned first aspect are implemented.

[0023] In a third aspect, the present application provides a computer-readable storage medium having stored thereon computer program instructions for releasing an internal spoiler device for an aneurysm, wherein when the computer program instructions are executed by one or more processors, one or more embodiments of the aforementioned first aspect are implemented.

[0024] Through a simulation method, simulation device and storage medium for releasing an internal spoiler device for an aneurysm as provided above, the embodiment of the present application constructs a spring mass system for the internal spoiler device according to the selection of the internal spoiler device of the aneurysm, and sets a target constraint on the spring mass system to optimize the position of the target mass in the spring mass system to achieve the release simulation of the internal spoiler device. Based on this, the embodiment of the present application constructs a spring mass system for the internal spoiler device to approximate the internal spoiler device to an inflatable cloth with limited expansion and contraction, and controls the size and deformation range of the internal spoiler device by setting the target constraint of the cloth. Since the internal spoiler device will contact the surface of the aneurysm when released, the embodiment of the present application corrects the position of the target mass in the spring mass system based on the target constraint to achieve accurate simulation of the release of the internal spoiler device, thereby improving the success rate of the selection and release of the internal spoiler device during the operation. Further, the embodiment of the present application does not require complex modeling, thereby greatly improving the simulation efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] By reading the detailed description below with reference to the accompanying drawings, the above and other purposes, features and advantages of the exemplary embodiments of the present application will become easy to understand. In the accompanying drawings, several embodiments of the present application are shown in an exemplary and non-limiting manner, and the same or corresponding reference numerals represent the same or corresponding parts, wherein:

[0026] Figure 1 An exemplary flow chart showing a simulation method for releasing an internal flow disrupting device for an aneurysm according to some embodiments of the present application;

[0027] Figure 2 An exemplary flow chart showing the operation of obtaining the selection of the internal spoiler device in some embodiments of the present application;

[0028] Figure 3 An exemplary schematic diagram showing reconstruction of an aneurysm surface and extraction of an aneurysm neck according to some embodiments of the present application;

[0029] Figure 4An exemplary schematic diagram of a fitted ellipsoid according to some embodiments of the present application is shown;

[0030] Figure 5 Shown are exemplary schematic diagrams of two different types of spring-mass systems for some embodiments of the present application;

[0031] Figure 6 An exemplary schematic diagram showing the simulated effect of releasing the inner spoiler device of some embodiments of the present application;

[0032] Figure 7 An overall exemplary flow chart of simulating the release of an internal spoiler device according to some embodiments of the present application is shown;

[0033] Figure 8 A schematic block diagram of a simulation device for releasing an internal flow disrupting device for an aneurysm according to some embodiments of the present application is shown. DETAILED DESCRIPTION

[0034] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present application.

[0035] It should be understood that the terms "include" and "comprising" used in the specification and claims of the present application indicate the presence of described features, wholes, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or collections thereof.

[0036] It should also be understood that the terms used in this application specification are only for the purpose of describing specific embodiments and are not intended to limit the present application. As used in this application specification and claims, unless the context clearly indicates otherwise, the singular forms of "a", "an" and "the" are intended to include plural forms. It should also be further understood that the term "and / or" used in this application specification and claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.

[0037] As used in this specification and claims, the term "if" may be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" may be interpreted as meaning "upon determination" or "in response to determining" or "upon detection of [described condition or event]" or "in response to detecting [described condition or event]," depending on the context.

[0038] Before further describing the content and scope of the internal spoiler release simulation of the present application, it is first necessary to clarify the relevant basic concepts and their uses.

[0039] Aneurysms are abnormal expansions of the blood vessel wall that may lead to rupture and serious consequences. Internal spoiler devices are an interventional treatment that slows down blood flow in the aneurysm by placing a device inside the blood vessel, thereby reducing the risk of rupture. The simulation method for the release of internal spoiler devices aims to evaluate the blood flow after the release of the internal spoiler devices through numerical simulation and other means, so as to guide clinical practice and improve treatment effects.

[0040] As can be seen from the above background technology description, in current vascular interventional surgeries, doctors manually measure blood vessels or aneurysms before surgery, and then rely on experience to select internal spoilers and release the internal spoilers. Doctors with different experiences have different evaluations of consumables selection and post-operative release effects. Therefore, this method has many problems, such as prolonging the operation time, deteriorating the quality of the operation, affecting the patient's prognosis, causing recurrence, etc. The existing WEB release simulation method is to simulate the release through the FEA method. This method is slightly more accurate, but time-consuming and difficult to apply in actual clinical practice.

[0041] Based on this, the present application proposes a simulation method for the release of an internal spoiler device for an aneurysm. By constructing a spring-particle system for the internal spoiler device, the internal spoiler device is approximated as an inflatable fabric with limited expansion and contraction. At the same time, the accuracy and efficiency of the simulation of the release of the internal spoiler device are taken into account, making it more suitable for actual clinical use.

[0042] The specific implementation of the present application is described in detail below with reference to the accompanying drawings.

[0043] Figure 1 FIG. 1 is an exemplary flow chart of a simulation method 100 for releasing an inner spoiler device for an aneurysm according to some embodiments of the present application. Figure 1 As shown in FIG. 1 , at step S101 , a spring mass point system for the internal flow disturbance device is constructed according to the selection of the internal flow disturbance device of the aneurysm, wherein the spring mass point system includes a plurality of mass points connected by a plurality of spring edges.

[0044] It can be understood that the spring-mass system is a classic continuum mechanics model used to describe the interaction between springs (elastic elements) and mass points (mass elements). Springs are usually modeled as linear springs, and their force and displacement satisfy Hooke's law. The mass point represents the mass point in the system, which is subject to the combined force of the forces from the adjacent springs and the action force.

[0045] In the spring-mass system, multiple particles connected by multiple spring edges form a complex mechanical structure. The connection relationship between these spring edges and particles can be represented and modeled by a triangular network on the surface of the inner spoiler. The triangular network is composed of multiple triangles, each of which is connected by three nodes (or vertices) and edges. Thus, multiple spring edges and multiple particles of the embodiment of the present application are connected to form a triangular network comprising multiple triangles, and the spring edges correspond to the edges in the triangular network, and the particles correspond to the vertices in the triangular network.

[0046] Before constructing a spring-mass system for an internal flow disturbance device, it is first necessary to obtain the selection of the internal flow disturbance device of the aneurysm. In one embodiment, the surface of the aneurysm can be reconstructed and the neck of the aneurysm can be extracted, and then the geometric parameters of the target geometric shape close to the aneurysm morphology can be obtained based on the aneurysm surface and neck fitting, and then the size of the aneurysm can be calculated according to the geometric parameters to obtain the selection of the internal flow disturbance device of the aneurysm. In some embodiments, the aforementioned target geometric shape can include an ellipsoid, the aforementioned geometric parameters include the main axis parameters of the ellipsoid, and the size of the aneurysm includes the height and width of the aneurysm, wherein the aforementioned main axis of the ellipsoid includes three, corresponding to the three directions of x, y, and z respectively.

[0047] That is, the embodiment of the present application first approximates the aneurysm morphology by fitting an ellipsoid whose main axis is perpendicular to the aneurysm neck based on the segmented aneurysm, and then calculates the height of the internal spoiler device by the length of the main axis (for example, the axis in the z direction) and the compression ratio of the internal spoiler device, and calculates the diameter of the internal spoiler device by the length of the other two axes (for example, the axes in the x and y directions). In some embodiments, in order to ensure that one of the axes is perpendicular to the plane of the aneurysm neck and to avoid the problem that the internal spoiler device cannot completely fit the aneurysm wall after release due to the large difference between the long and short main axes, the embodiment of the present application also proposes setting fitting constraints, and under the fitting constraints, obtains the geometric parameters of the target geometric shape close to the aneurysm morphology based on the fitting of the aneurysm surface and the aneurysm neck. This will be discussed later in conjunction with Figure 2 The aforementioned ellipsoid fitting process and calculation of the height and diameter of the inner spoiler are described in detail.

[0048] Next, at step S102, in the simulation of releasing the internal spoiler into the aneurysm, a target constraint is set for the spring mass system. In some embodiments, the target constraint may include at least a stretch constraint, a bending constraint, and an internal pressure constraint. Among them, the stretch constraint can be used to limit the elongation of the spring, the bending constraint can keep the shape and angle of the spring stable, and the internal pressure constraint can adjust the simulated pressure distribution and balance. It should be understood that the surface of the internal spoiler is modeled as a set of particles connected by springs by using a spring mass system, wherein the spring represents structure, shear and bending constraints. In the simulation of releasing the internal spoiler into the aneurysm, the triangular mesh of the internal spoiler (corresponding to the spring mass system) is accepted as input, and the triangular mesh is restricted to a manifold structure, that is, each edge is shared by at most two triangles. Thus, the embodiment of the present application sets the target constraint so that each node becomes a vertex in the simulated release, and the position of the target particle in the spring mass system is corrected in the simulation to achieve an accurate simulation of the release of the internal spoiler.

[0049] Specifically, in one implementation scenario, a stretch constraint is set on the edge of a triangle in a spring mass system, a bending constraint is set on adjacent triangles in the spring mass system, and an internal pressure constraint is set on the triangular network in the spring mass system. In one implementation scenario, a stretch constraint can be set on the edge of a triangle in a spring mass system by the following formula:

[0050] C stretch (P1, P2) = |P1-P2|-l0 (1)

[0051] Among them, C stretch (P1, P2) represents the stretch constraint, P1 and P2 represent the two vertices of the triangle edge connecting the spring-mass system, and l0 represents the length of the initial edge between the two vertices.

[0052] In one implementation scenario, the bending constraints can be set for adjacent triangles in the spring-mass system by the following formula:

[0053]

[0054] Among them, C bend (p1, p2, p3, p4) represents the bending constraint, p1, p2, p3, p4 represent the vertices of two adjacent triangles, for example, the adjacent triangle pairs (p1, p3, p2) and (p1, p2, p4), Represents the initial angle between adjacent triangle pairs.

[0055] By setting the above-mentioned stretching constraint and bending constraint, the internal spoiler can be regarded as a cloth with low stretching stiffness but high resistance to bending. In addition, since the internal spoiler is modeled as a closed triangular cloth mesh, the pressure of the mesh can be easily modeled for the closed triangular mesh, which is similar to the inflation process of a hot air balloon. Therefore, an internal pressure constraint can be set for the triangular network in the spring-mass system by adding an equality constraint on all N vertices of the mesh. In an implementation scenario, the aforementioned internal pressure constraint can be set by the following formula:

[0056]

[0057] Among them, C(P1,…,P N ) represents the internal pressure constraint, n triangles represents the number of triangles in the triangulated network, and represents the three vertices in each triangle, k pressure represents the pressure factor, and V0 represents the initial volume of the triangulated network.

[0058] Furthermore, after the internal spoiler is subjected to stretching constraint, bending constraint and internal pressure constraint, a unidirectional constraint may be set for the target vertex of the triangular network in the spring mass system. The aforementioned unidirectional constraint includes a first unidirectional constraint and a second unidirectional constraint. In one embodiment, in response to the target vertex of the triangular network in the spring mass system passing through the corresponding triangle, a first unidirectional constraint is set for the target vertex of the triangular network in the spring mass system. Or in response to the target vertex of the triangular network in the spring mass system merging from the bottom of the corresponding triangle, a second unidirectional constraint is set for the target vertex of the triangular network in the spring mass system. That is, the self-collision inside the cloth can be handled by an additional unidirectional constraint.

[0059] In one implementation scenario, for a target vertex q that passes through a triangle formed by vertices p1, p2, and p3, a first unidirectional constraint may be set for the target vertex of the triangular network in the spring-mass system by the following formula:

[0060]

[0061] Among them, C(q, P1, P2, P3) represents the first unidirectional constraint, q represents the target vertex, P1, P2, P3 represent the vertices of the corresponding triangles, and h represents the thickness of the triangular network in the spring-mass system.

[0062] In another implementation scenario, for a target vertex that flows in from the bottom of a triangle, a second one-way constraint can be set on the target vertex of the triangular network in the spring-mass system relative to the triangle normal by the following formula:

[0063]

[0064] Among them, C(q,P1,P2,P3)' represents the second unidirectional constraint, q represents the target vertex, P1,P2,P3 represent the vertices of the corresponding triangle, and h represents the thickness of the triangular network in the spring-mass system.

[0065] After setting the above target constraint for the spring mass system, at step S103, the position of the target mass in the spring mass system is optimized based on the target constraint to simulate the release of the inner spoiler. It can be understood that the size and deformation range of the inner spoiler can be controlled by the aforementioned stretching constraint, bending constraint, internal pressure constraint and unidirectional constraint. When the inner spoiler is released, it will come into contact with the surface of the aneurysm. By detecting the collision of the triangular facets of the spring mass system, the mass position on the inner spoiler is corrected to simulate the release of the inner spoiler.

[0066] In combination with the above description, it can be seen that the embodiment of the present application constructs a spring mass system for the internal spoiler according to the selection of the internal spoiler of the aneurysm, so as to approximate the internal spoiler as a limited and retractable inflatable fabric. Then, a target constraint is set for the spring mass system to optimize the position of the target mass in the spring mass system to control the size and deformation range of the internal spoiler. When the internal spoiler is released and contacts the surface of the aneurysm, the position of the target mass in the spring mass system is corrected to achieve accurate simulation of the release of the internal spoiler, thereby improving the success rate of the selection and release of the internal spoiler during surgery. In addition, the embodiment of the present application does not require complex modeling, thereby greatly improving the simulation efficiency.

[0067] Figure 2 FIG. 2 is an exemplary flow chart of operation 200 for obtaining the selection of the internal spoiler device in some embodiments of the present application. Figure 2 As shown in FIG. 2 , at step S201, the aneurysm surface is reconstructed and the aneurysm neck is extracted. It should be understood that an aneurysm is an abnormal expansion of a blood vessel wall, forming a balloon-like structure, in which the most critical part is the aneurysm neck, i.e., the narrow area connected to the blood vessel. In medical image processing, reconstructing the aneurysm surface and accurately extracting the aneurysm neck are of great significance for diagnosis and treatment.

[0068] In one implementation scenario, a vascular image may be first acquired and then enhanced (eg Figure 3(a)), and then use methods such as threshold processing, region growing or deep learning to obtain masks of blood vessels and aneurysm areas. Among them, the aforementioned threshold processing refers to binarizing the blood vessels and aneurysms in the image by setting a suitable grayscale value range for subsequent segmentation and reconstruction. The aforementioned region growing is an algorithm based on seed point growth, which gradually expands the region according to the similarity between pixels, thereby achieving accurate positioning of blood vessels and aneurysm areas. The aforementioned deep learning can realize automatic recognition and segmentation of different structures in the image by training models such as convolutional neural networks, thereby improving accuracy and efficiency.

[0069] Based on the mask obtained above, the Marching Cube algorithm can be used to reconstruct the surface of the blood vessel and aneurysm according to the masks of the blood vessel and aneurysm region to obtain a reconstructed blood vessel model (eg Figure 3 (b)) and aneurysm models (e.g. Figure 3 (c)). Among them, the Marching Cube algorithm is a commonly used 3D surface reconstruction algorithm, which realizes surface reconstruction of complex structures by converting the scalar field in the volume data into a polygonal mesh representation. The core of the algorithm is to perform operations such as interpolation and depiction on the scalar field in the volume data to determine the boundary conditions of each voxel and finally generate a triangular mesh representation of the surface, which can effectively handle irregular shapes and curved surfaces and is suitable for the reconstruction of vascular structures such as aneurysms. This surface reconstruction method based on the Marching Cube algorithm can help medical radiologists and clinicians observe and analyze the morphological characteristics of aneurysms more clearly, and provide an important reference for surgical planning and treatment.

[0070] After obtaining the aneurysm surface and the aneurysm neck, in step S202, geometric parameters of a target geometric shape close to the aneurysm morphology are obtained based on the aneurysm surface and the aneurysm neck fitting. As can be seen from the foregoing, the target geometric shape can be, for example, an ellipsoid, and the geometric parameters include the three main axis parameters of the ellipsoid (e.g. Figure 4 The size of an aneurysm includes the height and width of the aneurysm.

[0071] Specifically, in an implementation scenario, the above ellipsoid model can be represented by the following method:

[0072]

[0073] In some embodiments, the ellipsoid can be fitted to approximate the aneurysm morphology by, for example, a gradient descent method, and the goal of the fitting is to minimize the ellipsoid function represented by the above formula (6) by finding the parameters a, b and c. The ellipsoid function is defined as the distance between a point on the aneurysm surface and the ellipsoid surface, that is, (x, y, z) in the above formula (6) represents a point on the aneurysm surface, and a, b and c represent the three main axis parameters of the ellipsoid.

[0074] In order to ensure that one of the axes is perpendicular to the plane of the tumor neck, an additional fitting constraint can be introduced to make the direction of the ellipsoid's main axis parallel to the normal vector of the plane. In an implementation scenario, a direction fitting constraint for the main axis direction can be expressed as follows:

[0075]

[0076] Wherein, N represents the normal vector of the aneurysm neck plane, and V = [1 / a, 1 / b, 1 / c]. In addition, during the actual release of the internal spoiler, if an aneurysm with a relatively uneven width is encountered, the internal spoiler will not be able to fully adapt to the aneurysm and will deform. Therefore, the embodiment of the present application also sets a fitting constraint for the ratio of the two main axes of the aneurysm cross section. In an implementation scenario, the proportional fitting constraint for the ratio of the two main axes of the aneurysm cross section can be expressed as:

[0077] C ratio =(b-kc) 2 (8)

[0078] Where k represents the preset ratio of the major diameter to the minor diameter. Furthermore, the embodiment of the present application introduces the Lagrange multiplier to add the direction constraint and the ratio constraint as a penalty term to the ellipsoid function of the above formula (6) to obtain a, b and c. Where a is the height perpendicular to the plane of the aneurysm neck, and b and c are the lengths of the main axes parallel to the cross section of the aneurysm neck, so that the fitted ellipsoid tends to be perpendicular to the plane of the aneurysm neck, and there will be no large size difference in the cross-sectional direction of the aneurysm.

[0079] According to the geometric parameters a, b and c obtained above, at step S203, the size of the aneurysm is calculated according to the geometric parameters to obtain the selection of the internal flow disturbance device of the aneurysm. Specifically, in one embodiment, the size of the aneurysm can be calculated based on the following formula:

[0080]

[0081] H=a (9)

[0082] Wherein, W represents the width of the aneurysm, and H represents the height of the aneurysm. In some embodiments, the width of the flow-disturbing device within the aneurysm is obtained by rounding the value of W to an integer and then adding 1; the height of the flow-disturbing device within the aneurysm is obtained by rounding the value of H to an integer and then subtracting 1, so as to obtain the selection of the internal flow-disturbing device.

[0083] Figure 3 An exemplary schematic diagram of reconstructing the aneurysm surface and extracting the aneurysm neck of some embodiments of the present application is shown. Figure 3 Figure (a) shows the vascular enhancement image. Figure 3 Figure (b) shows the reconstructed blood vessel model. Figure 3 Figure (c) shows the reconstructed aneurysm model. As mentioned above, firstly, a method such as threshold processing, region growing or deep learning can be used to obtain a mask of the blood vessel and aneurysm region, and then based on the mask, the blood vessel and aneurysm surface is reconstructed using, for example, the Marching Cube algorithm to obtain the reconstructed blood vessel model and aneurysm model.

[0084] Figure 4 FIG. 2 shows an exemplary schematic diagram of a fitting ellipsoid in some embodiments of the present application. Figure 4 The parameters a, b and c shown in the figure represent the three main axis parameters of the ellipsoid. In the implementation scenario, the three main axis parameters of the ellipsoid similar to the aneurysm morphology are obtained through the above-mentioned fitting based on the aneurysm surface and aneurysm neck, and then the selection of the internal flow disturbance device of the aneurysm is obtained based on the main axis parameters. For more details about the ellipsoid fitting, please refer to the above Figure 2 The description of the internal spoiler device will not be repeated in this application. According to the foregoing, according to the selection of the internal spoiler device, a spring mass system for the internal spoiler device can be constructed, and by setting target constraints on the spring mass system to optimize the position of the target mass in the spring mass system, the release of the internal spoiler device can be accurately simulated. In some embodiments, for constructing a spring mass system for the internal spoiler device, the original braided structure of the internal spoiler device can be first reconstructed to obtain a cloth model, and then the spring mass system can be constructed based on the reconstructed cloth model, for example Figure 5 shown.

[0085] Figure 5 Schematic diagrams of two different types of spring mass systems for some embodiments of the present application are shown. Figure 5The upper and lower rows in the middle are, from left to right, the original weaving structure of the internal spoiler device of different selections, the reconstructed cloth model, and the corresponding spring mass system. As can be seen from the figure, the spring mass system includes multiple mass points connected by multiple spring edges, and multiple spring edges and multiple mass points are connected to form a triangular network including multiple triangles. In an implementation scenario, by setting stretching constraints on the edges of the triangles in the spring mass system based on the above formulas (1) to (5), setting bending constraints on the adjacent triangles in the spring mass system, setting internal pressure constraints on the triangular network in the spring mass system, and setting unidirectional constraints on the target vertices of the triangular network in the spring mass system, the internal spoiler device can be approximated as a finitely retractable inflatable cloth. The size and deformation range of the internal spoiler device are controlled by the bending and stretching constraints, internal pressure constraints, and unidirectional constraints of the cloth, so that when the internal spoiler device is released and contacts the surface of the aneurysm, the position of the mass points on the internal spoiler device is corrected by detecting the collision of the triangular facets, thereby achieving accurate simulation of the internal spoiler device. For more details on the above constraint settings, please refer to the above Figure 1 The description of is not repeated here in this application.

[0086] Figure 6 An exemplary schematic diagram showing the simulated effect of releasing the inner spoiler device of some embodiments of the present application is shown. Figure 6 The effect diagram of the internal spoiler device after being released in the aneurysm is shown in FIG. Utilizing the solution of the present application, the process of releasing the internal spoiler device can be simulated efficiently and accurately, so as to improve the success rate of selecting and releasing the internal spoiler device during surgery.

[0087] Figure 7 FIG. 1 shows an exemplary flow chart of the overall process of simulating the release of the internal spoiler device in some embodiments of the present application. Figure 7 As shown in , at step S701, blood vessel segmentation and reconstruction are performed. Specifically, the blood vessel enhanced image described above, the mask of the blood vessel and aneurysm area is obtained, and the blood vessel and aneurysm surface are reconstructed to obtain the reconstructed blood vessel model and aneurysm model. Then, at step S702, the aneurysm is segmented and measured. As can be seen from the foregoing, the geometric parameters of the ellipsoid similar to the aneurysm morphology can be obtained by fitting the ellipsoid, so as to obtain the height and width of the aneurysm according to the geometric parameters (such as a, b, c mentioned above).

[0088] Further, at step S703, WEB model selection is performed. For example, by rounding the W value to an integer and then adding 1, the width value of the intra-tumor spoiler device is used; the H value is rounded to an integer and then subtracted by 1, as the height value of the intra-tumor spoiler device, to obtain the selection of the internal spoiler device. Based on this selection, at step S704, WEB modeling is performed. That is, a spring mass system is constructed. Finally, at step S705, a WEB release simulation is performed. Specifically, the position of the target mass can be optimized by setting a target constraint on the spring mass system to achieve an accurate simulation of the release of the internal spoiler device.

[0089] Figure 8 FIG. 8 is a schematic block diagram of a simulation device 800 for releasing an inner spoiler device for an aneurysm according to some embodiments of the present application. Figure 8 As shown, the simulation device 800 may include a processor 801 and a memory 802, wherein the processor 801 and the memory 802 communicate with each other via a bus. The memory 802 stores computer instructions for simulating the release of the internal spoiler device of the aneurysm. When the computer instructions are executed by the processor 801, the device 800 implements the method steps described in the above text in combination with the accompanying drawings: constructing a spring mass point system for the internal spoiler device according to the selection of the internal spoiler device of the aneurysm, wherein the spring mass point system includes a plurality of mass points connected by a plurality of spring edges; setting a target constraint for the spring mass point system in simulating the release of the internal spoiler device into the aneurysm; and optimizing the position of the target mass point in the spring mass point system based on the target constraint to simulate the release of the internal spoiler device.

[0090] According to the above description in combination with the accompanying drawings, those skilled in the art can also understand that the embodiments of the present application can also be implemented by a software program. Therefore, the present application also provides a computer-readable storage medium. The computer-readable storage medium stores thereon computer-readable instructions for simulating the release of the inner spoiler device of the aneurysm. When the computer-readable instructions are executed by one or more processors, the present application in combination with the accompanying drawings can be implemented. Figure 1 A simulated method for the deployment of an endovascular flow disruptor device for an aneurysm is described.

[0091] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0092] It should be noted that although the operations of the method of the present application are described in a specific order in the accompanying drawings, this does not require or imply that the operations must be performed in this specific order, or that all the operations shown must be performed to achieve the desired results. On the contrary, the steps depicted in the flow chart can be performed in a different order. Additionally or alternatively, some steps can be omitted, multiple steps can be combined into one step, and / or one step can be decomposed into multiple steps.

[0093] It should be understood that when the terms "first", "second", "third" and "fourth" are used in the claims, the specification and the drawings of the present application, they are only used to distinguish different objects, rather than to describe a specific order. The terms "include" and "comprise" used in the specification and claims of the present application indicate the presence of the described features, wholes, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or their collections.

[0094] It should also be understood that the terms used in this application specification are only for the purpose of describing specific embodiments and are not intended to limit the present application. As used in this application specification and claims, unless the context clearly indicates otherwise, the singular forms of "a", "an" and "the" are intended to include plural forms. It should also be further understood that the term "and / or" used in this application specification and claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.

[0095] Although the implementation methods of the present application are as above, the contents described are only examples adopted to facilitate the understanding of the present application, and are not intended to limit the scope and application scenarios of the present application. Any technician in the technical field described in the present application can make any modifications and changes in the form and details of the implementation without departing from the spirit and scope disclosed in the present application, but the scope of patent protection of the present application shall still be subject to the scope defined in the attached claims.

[0096] In addition, the collection and acquisition of various data in this application complies with relevant laws and regulations and is authorized by the data provider. Any organization or individual that needs to obtain external data must obtain authorization in accordance with the law and ensure data security. It is not allowed to illegally collect, use, process, or transmit unauthorized or unprotected data, or to illegally buy, sell, provide, or disclose unauthorized or unprotected data.

Claims

1. A method for simulating the release of an internal flow disturbance device for an aneurysm, comprising: According to the selection of the internal flow disturbance device of the aneurysm, a spring mass point system for the internal flow disturbance device is constructed, wherein the spring mass point system comprises a plurality of mass points connected by a plurality of spring edges; In simulating the release of the internal spoiler into an aneurysm, setting a target constraint on the spring-mass system; as well as optimizing the position of the target mass point in the spring mass point system based on the target constraint to simulate the release of the internal spoiler, Wherein, constructing a spring mass system for the internal flow spoiler according to the selection of the internal flow spoiler of the aneurysm includes: The cloth model is obtained by reconstructing the original weaving structure of the internal spoiler device according to the selection of the internal spoiler device of the aneurysm; The plurality of spring edges and the plurality of mass points are connected based on the reconstructed cloth model to form a triangular network including a plurality of triangles, so as to construct a spring mass point system.

2. The simulation method according to claim 1, wherein the target constraint comprises a stretch constraint, a bending constraint and an internal pressure constraint, and setting the target constraint for the spring-mass system comprises: Setting the stretch constraint on the side of the triangle in the spring-mass system; Setting the bending constraint on adjacent triangles in the spring-mass system; as well as The internal pressure constraint is set on the triangular network in the spring-mass system.

3. The simulation method according to claim 2, wherein the stretching constraint is set by the following formula: in, represents the stretch constraint, and denote the two vertices connecting the sides of the triangle in the spring-mass system, Indicates the length of the initial edge.

4. The simulation method according to claim 2, wherein the bending constraint is set by the following formula: in, represents the bending constraint, represents the vertices of two adjacent triangles, Indicates the initial angle.

5. The simulation method according to claim 2, wherein the internal pressure constraint is set by the following formula: in, represents the internal pressure constraint, represents the number of triangles in the triangulated network, and Represents the vertices in each triangle ， represents the pressure factor, represents the initial volume.

6. The simulation method according to claim 2, further comprising: A unidirectional constraint is set on a target vertex of a triangular network in the spring-mass system.

7. The simulation method according to claim 6, wherein the one-way constraint comprises a first one-way constraint and a second one-way constraint, and setting the one-way constraint on the target vertex of the triangular network in the spring-mass system comprises: In response to a target vertex of the triangular network in the spring-mass system passing through a corresponding triangle, setting the first unidirectional constraint on the target vertex of the triangular network in the spring-mass system; or In response to a target vertex of the triangular network in the spring-mass system being introduced from below a corresponding triangle, the second unidirectional constraint is set on the target vertex of the triangular network in the spring-mass system.

8. The simulation method according to claim 7, wherein the first unidirectional constraint is set by the following formula: in, represents the first unidirectional constraint, represents the target vertex, denote the vertices of the corresponding triangle, Represents the thickness of the triangular network in the spring-mass system.

9. The simulation method according to claim 7, wherein the second unidirectional constraint is set by the following formula: in, represents the second unidirectional constraint, represents the target vertex, denote the vertices of the corresponding triangle, Represents the thickness of the triangular network in the spring-mass system.

10. The simulation method according to claim 1, wherein the selection of the internal flow disturbance device of the aneurysm is obtained by the following operations: Reconstruct the aneurysm surface and extract the aneurysm neck; Obtaining geometric parameters of a target geometric shape close to the aneurysm morphology based on the aneurysm surface and the aneurysm neck fitting; and The size of the aneurysm is calculated according to the geometric parameters to obtain the selection of the internal flow disturbance device of the aneurysm. 11 . The simulation method according to claim 10 , wherein the target geometric shape comprises an ellipsoid, the geometric parameters comprise the major axis parameters of the ellipsoid, and the size of the aneurysm comprises the height and width of the aneurysm.

12. The simulation method according to claim 10, further comprising: Set fitting constraints; as well as Under the fitting constraints, the geometric parameters of the target geometric shape close to the aneurysm morphology are obtained based on the aneurysm surface and the aneurysm neck fitting.

13. A simulation device for releasing an internal flow disturbance device for an aneurysm, comprising: processor; as well as A memory having computer instructions for releasing an internal spoiler device for an aneurysm stored thereon, wherein when the computer instructions are executed by a processor, the simulation method according to any one of claims 1-12 is implemented.

14. A computer-readable storage medium having stored thereon computer program instructions for releasing an internal flow disturbance device for an aneurysm, wherein when the computer program instructions are executed by one or more processors, the simulation method according to any one of claims 1 to 12 is implemented.

Images