CN113178001A - Silicone oil filling simulation method for pore-source retinal detachment and electronic equipment - Google Patents
Silicone oil filling simulation method for pore-source retinal detachment and electronic equipment Download PDFInfo
- Publication number
- CN113178001A CN113178001A CN202110353522.XA CN202110353522A CN113178001A CN 113178001 A CN113178001 A CN 113178001A CN 202110353522 A CN202110353522 A CN 202110353522A CN 113178001 A CN113178001 A CN 113178001A
- Authority
- CN
- China
- Prior art keywords
- silicone oil
- oil filling
- eyeball
- interaction
- liquid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 229920002545 silicone oil Polymers 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000004088 simulation Methods 0.000 title claims abstract description 60
- 206010038848 Retinal detachment Diseases 0.000 title claims abstract description 26
- 230000004264 retinal detachment Effects 0.000 title claims description 15
- 210000005252 bulbus oculi Anatomy 0.000 claims abstract description 62
- 239000007788 liquid Substances 0.000 claims abstract description 57
- 239000012530 fluid Substances 0.000 claims abstract description 52
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 48
- 230000003993 interaction Effects 0.000 claims abstract description 42
- 210000001508 eye Anatomy 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 230000000694 effects Effects 0.000 claims abstract description 18
- 208000019793 rhegmatogenous retinal detachment Diseases 0.000 claims abstract description 11
- 230000016507 interphase Effects 0.000 claims abstract description 10
- 230000000007 visual effect Effects 0.000 claims abstract description 9
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims abstract description 6
- 238000005859 coupling reaction Methods 0.000 claims abstract description 6
- 239000002245 particle Substances 0.000 claims description 53
- 238000004458 analytical method Methods 0.000 claims description 13
- 230000033001 locomotion Effects 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 230000015654 memory Effects 0.000 claims description 10
- 230000005514 two-phase flow Effects 0.000 claims description 10
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 229910052710 silicon Inorganic materials 0.000 claims description 9
- 239000010703 silicon Substances 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 239000003921 oil Substances 0.000 claims description 7
- 229920001296 polysiloxane Polymers 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 230000008846 dynamic interplay Effects 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000003709 image segmentation Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000003909 pattern recognition Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 16
- 230000002980 postoperative effect Effects 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000012800 visualization Methods 0.000 abstract description 2
- 238000010276 construction Methods 0.000 abstract 1
- 239000012071 phase Substances 0.000 description 23
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 238000005520 cutting process Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 239000011521 glass Substances 0.000 description 5
- 210000001525 retina Anatomy 0.000 description 5
- 230000002207 retinal effect Effects 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 238000011282 treatment Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 4
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 210000004127 vitreous body Anatomy 0.000 description 3
- 208000002367 Retinal Perforations Diseases 0.000 description 2
- 210000000695 crystalline len Anatomy 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000004410 intraocular pressure Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 208000010412 Glaucoma Diseases 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 208000002158 Proliferative Vitreoretinopathy Diseases 0.000 description 1
- 201000007737 Retinal degeneration Diseases 0.000 description 1
- 206010038897 Retinal tear Diseases 0.000 description 1
- 206010038934 Retinopathy proliferative Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 210000002159 anterior chamber Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000001742 aqueous humor Anatomy 0.000 description 1
- 239000008346 aqueous phase Substances 0.000 description 1
- 239000003855 balanced salt solution Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 230000003412 degenerative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000004945 emulsification Methods 0.000 description 1
- 210000005081 epithelial layer Anatomy 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 235000003642 hunger Nutrition 0.000 description 1
- 238000012482 interaction analysis Methods 0.000 description 1
- 206010023365 keratopathy Diseases 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 230000008338 local blood flow Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 208000021971 neovascular inflammatory vitreoretinopathy Diseases 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000000049 pigment Substances 0.000 description 1
- 230000006785 proliferative vitreoretinopathy Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000037351 starvation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Computer Graphics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Prostheses (AREA)
Abstract
The invention discloses a silicone oil filling simulation method and electronic equipment for rhegmatogenous retinal detachment, wherein the method comprises the following steps: acquiring continuous medical scanning images of eyeballs, and constructing an eyeball three-dimensional model based on the medical scanning images; based on the eyeball three-dimensional model, carrying out stress simulation on the interaction boundary of the liquid and the intraocular cavity; and analyzing the interphase coupling state of the silicone oil and the water in the eye cavity, expressing the interaction between fluids by using a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation by combining stress simulation results of the silicone oil and the interaction boundary of the water and the eye cavity to obtain the theoretically optimal silicone oil filling amount. According to the invention, the theoretical optimal silicone oil filling amount is obtained by customized eyeball model construction and non-compressible non-divergence hydrodynamic calculation combined with surface tension, so that the method can play an auxiliary role in the surgical process, reduce the occurrence probability of postoperative silicone oil complications, and has the advantages of high calculation precision, vivid visualization effect and the like.
Description
Technical Field
The invention relates to the technical field of computer graphics fluid simulation, in particular to a silicone oil filling simulation method and electronic equipment for pore-source retinal detachment.
Background
Retinal Detachment (RRD) is a common vitreoretinal disease, usually due to Retinal atrophy degeneration or vitreous traction to form a Retinal neuroepithelial full-thickness hole, causing degenerative liquefaction of vitreous body to enter the subretinal space through the hole, resulting in separation between the Retinal neuroepithelial layer and the pigment epithelial layer. One of the main treatments for rhegmatogenous retinal detachment is through vitreous cutting in combination with silicone intraocular tamponade to achieve the purpose of retinal reattachment. The vitreous is a colorless transparent gel, approximately 4.5ml in adult vitreous volume, located in the vitreous cavity between the lens and retina, and occupies approximately 80% of the intraocular volume. The vitreous body is a component of an eye dioptric medium, has three characteristics of viscoelasticity, permeability and transparency, and has supporting and nutritional effects on peripheral tissues such as crystalline lens, retina and the like.
Among various vitreous substitutes, silicone oil has been widely used as a safe and effective intraocular filling agent due to its advantages of stable physicochemical properties, good biological tolerance, etc. However, in the combined operation of vitreous cutting and silicone oil filling, silicone oil needs to fill the whole vitreous cavity, which is easy to cause a plurality of complications, such as proliferative vitreoretinopathy silicone oil entering into the anterior chamber to cause keratopathy due to silicone oil emulsification, pupil-blocking glaucoma due to excessive filling amount, and the like. Therefore, how to achieve the best therapeutic effect with the minimum amount of silicone oil in vitreoretinal surgery is a problem which needs to be solved urgently.
After a patient with porogenic retinal detachment receives the treatment of vitreous cutting and silicone oil filling, because the density of silicone oil is slightly lower than that of water, the patient usually needs to keep a face-down position in order to press retinal holes and prevent water from entering under the retina and effectively repair the retina. At the moment, water and oil in the vitreous cavity are layered, the water is positioned below, and the silicone oil is positioned above and completely covers the retinal hole.
At present, the prior art has no effective solution for how to quantitatively simulate and analyze and visualize the state of an intraocular cavity and determine the ideal silicone oil filling amount in the operation in the vitreous cutting and silicone oil filling combined operation for treating the porogenic retinal detachment.
Disclosure of Invention
The invention provides a silicone oil filling simulation method and electronic equipment for pore-derived retinal detachment, and solves the problems that in the prior art, the vitreous body cutting and silicone oil filling operation for treating the pore-derived retinal detachment cannot be used for quantitatively simulating, analyzing and visualizing the state of an intraocular cavity and determining the ideal silicone oil filling amount in the operation.
In order to solve the technical problems, the invention provides the following technical scheme:
on one hand, the invention provides a silicone oil filling simulation method facing porogenic retinal detachment, which comprises the following steps:
acquiring continuous medical scanning images of eyeballs, and constructing an eyeball three-dimensional model based on the medical scanning images;
based on the eyeball three-dimensional model and the liquid physical characteristics, carrying out stress simulation on the interaction boundary of the liquid and the eye cavity; wherein the liquid comprises water and silicone oil in the ocular lumen;
and analyzing the interphase coupling state of the silicone oil and the water in the eye cavity, expressing the interaction between fluids by using a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation by combining stress simulation results of the silicone oil and the interaction boundary of the water and the eye cavity to obtain the theoretically optimal silicone oil filling amount.
Further, constructing an eyeball three-dimensional model based on the medical scanning image, comprising the following steps:
preprocessing the medical scanning image, wherein the preprocessing comprises smoothing denoising and gamma correction;
performing image segmentation on the preprocessed medical scanning image based on a preset region segmentation algorithm, performing pattern recognition on an eye structure, and completely describing a segmentation interval and eyeball tissues to obtain an eye structure image;
completing three-dimensional reconstruction of an eyeball tissue structure based on the eye structure image;
and calculating the volume of the communicated part of the reconstructed three-dimensional model, and removing the part of the communicated part with the volume smaller than a preset threshold value to finish the correction of the reconstructed three-dimensional model to obtain the eyeball three-dimensional model.
Further, based on the eyeball three-dimensional model and the liquid physical characteristics, the method for performing force-bearing simulation on the interaction boundary of the liquid and the eye cavity comprises the following steps:
based on the physical characteristics of the liquid, the three-dimensional eyeball model is discretized into boundary particles, the boundary particles are added into the approximate calculation of the liquid density, and the stress simulation is carried out on the interaction boundary of the liquid and the intraocular cavity.
Further, based on the physical characteristics of the liquid, the three-dimensional eyeball model is discretized into boundary particles, the boundary particles are added into the approximate calculation of the density of the liquid, and the interaction boundary of the liquid and the intraocular cavity is subjected to stress simulation, wherein the stress simulation comprises the following steps:
uniformly sampling the eyeball three-dimensional model by using a three-dimensional point cloud sampling algorithm;
giving the particle volume and mass to the sampling point;
the method comprises the steps of calculating the weight value of eyeball boundary particles borne by fluid particles in an eye cavity in physical calculation, considering the adhesion effect generated by interaction among different materials, calculating the pressure and viscosity of the interaction boundary of liquid and the eye cavity based on Newton's second law, and carrying out stress simulation on the interaction boundary of the liquid and the eye cavity.
Further, performing phase-to-phase coupling state analysis on the silicone oil and the water in the inner cavity of the eye, expressing fluid interaction by using a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation by combining stress simulation results of the silicone oil and the interaction boundary of the water and the inner cavity of the eye to obtain the theoretically optimal silicone oil filling amount, wherein the visual simulation results comprise:
performing single-phase fluid motion simulation on water and silicon oil in the eye cavity by using a smooth particle fluid dynamics method based on a Navier-Stokes equation;
the viscous force, the dragging force and the intersolubility of different flow phases are approximately calculated by using a kernel function, the interaction state of water and a silicon oil mixed liquid in a silicon oil filling state is simulated, and the silicon oil-water dynamic interaction calculation under the property of approximately incompressible non-dispersive fluid is carried out;
and (3) combining with an Euler dynamics method, continuously calculating the size of the surface tension direction of the space in a grid space based on a volume fraction method, and performing surface curvature analysis combining with the interphase surface tension of the two-phase flow to realize the simulation of the operation state under different silicone oil filling amounts, so as to generate the visual simulation result of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation, thereby obtaining the theoretically optimal silicone oil filling amount.
Further, the approximate silicone oil-hydrodynamic interaction calculation under incompressible non-dispersive fluid properties comprises:
performing density field calculation on the fluid distribution state in the discrete time interval to obtain the current flow field compression state;
calculating the physical field quantity of each macroscopic fluid particle position by using a smooth particle fluid dynamics method according to the fluid description in a compression state;
and calculating the pressure field quantity at each moment by using a pressure force suppression fluid compression characteristic principle and by explicitly solving a pressure state equation or an implicit iteration method.
Further, the surface curvature analysis in combination with the surface tension between the two-phase flow phases comprises:
extracting a surface boundary by setting a threshold value by adopting a color field method, calculating the surface curvature of the particle, and acquiring a unit normal vector;
the shrinkage degree of the interface is controlled by the tension coefficients of different substances, the normal direction of the interface is the stress direction, and the tension is in direct proportion to the surface curvature;
simulating the intermolecular attraction and repulsion characteristics of the fluid, setting an action threshold value, enabling the particles to attract each other when the distance between the macro discrete particles is lower than the threshold value, and enabling the particles to repel each other when the distance between the particles exceeds the threshold value.
In yet another aspect, the present invention also provides an electronic device comprising a processor and a memory; wherein the memory has stored therein at least one instruction that is loaded and executed by the processor to implement the above-described method.
In yet another aspect, the present invention also provides a computer-readable storage medium having at least one instruction stored therein, the instruction being loaded and executed by a processor to implement the above method.
The technical scheme provided by the invention has the beneficial effects that at least:
1) according to the invention, after related parameters of the eyeball (such as intraocular pressure, diameter of the eyeball, upper and lower diameters of a vitreous cavity, a horizontal diameter, a front and back diameter and the like) are obtained through a series of special examinations, a set of standard three-dimensional eyeball model is established, and the high-dimensional medical image information is realized.
2) The invention applies the method for calculating the hydrodynamics numerical value to the actual medical auxiliary process, and realizes interdisciplinary development. The relation between the filling amount of the silicone oil and the contact area of the retina is analyzed by combining the factors of intraocular pressure, the surface tension of the silicone oil, the specific gravity of a vitreous cavity balanced salt solution, the specific gravity of the silicone oil and the like, so that the minimum amount of the silicone oil required to be filled for repairing the porogenic retinal detachment is obtained, the occurrence of postoperative silicone oil complications is reduced to the greatest extent, and a safer and more effective novel operation treatment method is provided for treating the porogenic retinal detachment and other vitreoretinal diseases.
3) The technology of the invention has strong portability. The method for analyzing the multiphase fluid dynamics of the intraocular cavity can be transplanted to the related calculation and analysis of fluids in other biological tissues, such as local blood circulation, contrast agent injection, cerebrospinal fluid buffer effect simulation and the like.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a silicone oil filling simulation method for rhegmatogenous retinal detachment provided in an embodiment of the present invention;
fig. 2 is a schematic view of an intra-eyeball solid-liquid interaction analysis provided in the embodiment of the present invention;
FIG. 3 is a schematic diagram of interaction between two-phase liquid silicone oil and water provided by an embodiment of the invention;
FIG. 4 is a graph illustrating the effect of surface tension on various parameters provided by an embodiment of the present invention;
FIG. 5 is a graph illustrating the surface area variation for different frames according to an embodiment of the present invention;
FIG. 6 is a diagram illustrating an interaction simulation effect of two-phase liquid inside an eyeball according to an embodiment of the present invention;
fig. 7 is a schematic overall implementation flow diagram of a silicone oil filling simulation method for rhegmatogenous retinal detachment provided in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
First embodiment
In order to observe the structure of an eyeball and verify the angle and the detachment state of a retinal fissure at multiple angles, the simulation of the interaction environment of silicone oil-water two-phase liquid in an eye cavity after vitrectomy is needed, therefore, the embodiment provides a silicone oil filling simulation method for hole-source retinal detachment. Specifically, the execution flow of the silicone oil filling simulation method for rhegmatogenous retinal detachment of the embodiment is shown in fig. 1 and 7, and includes the following steps:
s1, acquiring continuous medical scanning images of the eyeball, and constructing a three-dimensional eyeball model based on the medical scanning images;
it should be noted that in the present embodiment, in the whole scheme, the image is preprocessed in a smooth denoising and gamma correction manner, then the region of interest is segmented by using an image segmentation algorithm based on watershed, three-dimensional reconstruction is performed on the segmented region, then the reconstructed three-dimensional model is trimmed, the volume of the connected portion is calculated, the smaller volume portion and the irrelevant content are removed, and the region of interest is retained. The continuous medical scanning image of the eyeball comprises an X-ray beam, gamma rays, ultrasonic waves and nuclear magnetic resonance used by the electronic computer tomography or all images which can be used for carrying out medical imaging after the scanning of human organs.
Specifically, the implementation process of S1 includes the following steps:
s11, obtaining eyeball related images through eye Computed Tomography (CT)/Magnetic Resonance Imaging (MRI) and the like, processing eye scanning images by using a noise reduction algorithm and gamma correction, and performing smooth correction processing;
s12, constructing a region segmentation algorithm, carrying out pattern recognition on the eye structure, and completely describing a segmentation interval and eyeball tissues;
and S13, analyzing the connectivity and distribution state of each segmented region of the image, reconstructing an eyeball model and acquiring the three-dimensional structure of the eyeball.
S2, performing stress simulation on the interaction boundary of the liquid and the eye cavity based on the eyeball three-dimensional model and the liquid physical characteristics; wherein the liquid comprises water and silicone oil in the ocular lumen;
it should be noted that, in this embodiment, based on the physical characteristics of the liquid, the solid model is discretized into boundary particles, and the boundary particles are added to the approximate calculation of the density of the liquid to perform the force simulation on the interaction boundary between the liquid and the intraocular cavity; as shown in fig. 2. Sampling the surface of a rigid object by using boundary particles, so that the rigid object can cope with eyeball characteristics of different shapes, wherein the rigid object comprises a low-dimensional rigid body formed by one layer or one row of boundary particles; while mitigating the problem of particle starvation near the boundary caused by the SPH approach, thereby preventing density (and hence pressure) discontinuities and particle adhesion artifacts at the boundary. Nearby boundary particles are considered in calculating the sum force of the fluid particles to avoid unrealistic adhesion phenomena.
Specifically, the implementation process of S2 includes the following steps:
s21, uniformly sampling the eyeball model by using a three-dimensional point cloud sampling algorithm, and endowing the sampling points with particle volume and mass;
s22, calculating the weight value of eyeball boundary particles borne by the intraocular fluid particles in physical calculation, considering the adhesion effect generated by the interaction between different materials, and calculating the boundary pressure and the viscosity based on the Newton' S second law.
And S3, performing interphase coupling state analysis on the silicone oil and the water in the eye cavity, expressing fluid interaction in a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation by combining stress simulation results of the silicone oil and the interaction boundary of the water and the eye cavity to obtain the theoretically optimal silicone oil filling amount.
It should be noted that, for the above S3, as shown in fig. 3, the implementation process is as follows:
a numerical value discretization model is constructed, dynamic analysis is carried out on the substances in a particle form through an SPH method, a linear iteration solving algorithm is designed for multiphase fluid motion, passive expression of a velocity field is achieved for incompressible viscous Newton fluid motion, and high-order motion precision of Lagrange particle expression is guaranteed. And analyzing the interactive movement between the silicone oil and the water two-phase liquid, marking the distribution condition and the movement state of different liquid phases, and constructing a miscible dragging model based on convection movement. And modeling by using volume fraction ratio to interfacial surface tension of the two-phase liquid to realize continuous and stable tangential tension expression. A constraint simulation model is designed aiming at the liquid surface tension effect under the eyeball small Weber number environment, and the comprehensive modeling of the Rayleigh instability phenomenon is realized.
Specifically, the implementation process of S3 includes the following steps:
s31, performing single-phase fluid motion simulation on water and silicone oil in the inner cavity of the eye based on a Navier-Stokes Equations (by using a Smooth Particle Hydrodynamics (SPH) method;
s32, using a kernel function to approximately calculate viscosity, drag force and intersolubility of different flow phases, simulating the interaction state of the mixed liquid in the filling state, and performing silicone oil-water dynamic interaction calculation under the property of approximate incompressible non-divergence fluid;
and S33, combining with an Euler dynamics method, continuously calculating the size of the space surface tension direction in a grid space based on a volume fraction method, and performing surface curvature analysis combining the two-phase flow interphase surface tension to realize the operation state simulation under different injection quantities.
It should be noted that, in the coupling process between silicone oil and water, two key technologies are adopted in this embodiment: firstly, approximate silicone oil-water dynamic interactive calculation under the property of incompressible non-dispersive fluid; and secondly, combining surface curvature analysis of the surface tension between two-phase flow phases. Wherein the content of the first and second substances,
the silicon oil-water dynamic interaction calculation under the property of approximate incompressible non-divergence fluid comprises the following steps:
s321, compressibility calculation: performing density field calculation on the fluid distribution state in the discrete time interval to obtain the current flow field compression state;
s322, calculating a physical field: calculating the physical field quantity of each macroscopic fluid particle position by using a smooth particle fluid dynamics method according to the fluid description in a compression state;
s323, pressure calculation: and calculating the pressure field quantity at each moment by using a pressure force suppression fluid compression characteristic principle and by explicitly solving a pressure state equation or an implicit iteration method.
Hydrodynamic methods using smooth particles based on fluid description under compression
Surface curvature analysis in combination with two-phase flow interphase surface tension comprising the steps of:
s331, interface extraction: extracting a surface boundary by setting a threshold value by adopting a color field method, calculating the surface curvature of the particle, and acquiring a unit normal vector;
s332, surface area minimization tension calculation: the shrinkage degree of the interface is controlled by the tension coefficients of different substances, the normal direction of the interface is the stress direction, and the tension is in direct proportion to the surface curvature;
s333, calculating cohesive tension: simulating fluid intermolecular attractive and repulsive force characteristics, and setting an action threshold value, wherein when the distance between macro discrete particles is lower than the threshold value, the particles are attracted to each other, and when the distance between the macro discrete particles exceeds the threshold value, the particles are repelled to each other.
In this embodiment, the effectiveness of the two-phase flow interaction model is verified through experiments. In the two-phase dam break experiment shown in fig. 3, the silicone oil and the water are two-phase fluids respectively, and at the beginning of the scene, the two fluids fall due to the action of gravity at the beginning of the simulation and start to contact and gradually mix in the process. In the process, the silicone oil phase continuously floats upwards due to the action of pressure force and interfacial force because of relatively small static density; while the aqueous phase with the greater inertia sinks continuously. The process is completed after the kinetic energy is completely dissipated, the two-phase fluid is completely layered, and has an obvious interface.
The experiment can express the eyeball silicone oil filling process in real time, when filling is carried out, silicone oil interacts with water in the eyeball to generate floating and sinking of fluid, and the fluid is divided into two layers in the eyeball to repair the retinal fracture.
In this example, the effectiveness of the surface tension method was verified through experiments. Figure 4 shows a schematic of the splash formed after a water column has fallen on a horizontal plate. The water column gradually horizontally spreads after falling on the plane, is scattered around due to the influence of residual kinetic energy, and finally stands still.
As can be seen in FIG. 4, the aggregation and non-uniform arrangement of the water droplets. When the surface tension coefficient is set to 0 (no surface tension effect), it can be observed that the liquid is uniformly distributed, covers most of the glass plate, and shows that the water drops are not aggregated and have no height protrusion. In the process of gradually increasing the surface tension coefficient, taking frame 22 as an example, the liquid splashing range gradually becomes smaller, because the effect of the surface tension is larger, the liquid cohesion is also increased, so that the water drops tend to gather. From frame 503 (at rest), it can be seen that the glass sheet is covered with less and less area by the same volume of liquid due to the increased coefficient, which laterally illustrates that surface tension is the primary cause of increased height of liquid accumulation.
The surface tension test (liquid falling on glass plate) shows the broken line of the sum of the liquid surface area with time for different parameters as shown in fig. 5. The vertical axis represents surface area and the horizontal axis represents time frame. From 0-100 frames, the surface area value is increased sharply because the liquid falls on the glass plate and impacts, and then the surface area is reduced for a short time because the liquid rolls back to touch the four walls of the glass plate. In the absence of surface tension (γ ═ 0), no significant oscillation of the surface area values occurred during this time period, while in the case of surface tension (γ ═ 0.1, 0.5, 0.8), the surface area values appeared to noticeably oscillate like a "w", indicating a cancellation of the surface tension effects from other external forces, pressures, during this time.
The present example is illustrated by fig. 5, where after 100 frames the liquid motion starts to gradually level off, which is where surface tension effects start to act as one of the main forces. At γ ═ 0, the liquid was free to spread, resulting in a gradual increase in the sum of the surface areas; as the gamma parameter becomes larger, the surface tension force becomes more and more significant, and thus the sum of the surface areas of the liquids also tends to be smaller and also stabilizes more quickly.
In the embodiment of the invention, an eyeball water and silicone oil two-phase interaction simulation experiment is carried out, as shown in fig. 6, in three groups of experiments, white represents silicone oil, and transparent represents water. Wherein (a) is that both phases have no surface tension effect, (b) is that only silicone oil has surface tension effect, gamma silicone oil is 0.85, and (c) is that both phases have surface tension effect, gamma silicone oil is 0.85, gamma silicone oil is 0.5. It was found that the silicone oil floats up in (a) without the action of surface tension, but a phenomenon of low aggregation occurs, and the white silicone oil is in a dispersed state as compared with (b) (c). From (b), it can be found that under the condition that water has no surface tension, the boundary between water and silicone oil is not obvious. Under the condition that both phases are provided with proper surface tension, (c) the boundary between the two phases is obvious, and the silicone oil is in a stable aggregation state and has a protrusion phenomenon; meanwhile, the gathering condition of the water below is more reasonable.
In summary, the present embodiment proposes a silicone oil-water dynamic interaction calculation under the property of approximately incompressible non-dispersive fluid combined with surface curvature analysis of surface tension between two-phase flow phases, so as to simulate the postoperative intraocular dynamics process and visualization process in the operation, and assist the doctor to achieve the best treatment effect with the least amount of silicone oil and reduce the occurrence of postoperative silicone oil complications. The method of the embodiment can reduce the filling usage amount of the silicone oil in eyes as much as possible and reduce complications on the premise of effectively repairing the retinal detachment. Aiming at the problems that the silicone oil filling is difficult to control and the filling state in the eyeball is difficult to master, the embodiment adopts a fluid simulation method based on physics to simulate the silicone oil filling in the eyeball. The method comprises the steps of establishing an eyeball three-dimensional model by using medical detection data, sampling particles, modeling water-silicone oil two-phase flow and interaction according to different physical properties of water and silicone oil, and finally establishing a solid-liquid interaction model to realize interactive simulation of two-phase liquid (water and silicone oil) and a solid boundary (eyeball). The method of the embodiment can better present the effects of multiphase motion, surface tension, boundary treatment and the like in the eyeball.
Second embodiment
The present embodiment provides an electronic device, which includes a processor and a memory; wherein the memory has stored therein at least one instruction that is loaded and executed by the processor to implement the method of the first embodiment.
The electronic device may have a relatively large difference due to different configurations or performances, and may include one or more processors (CPUs) and one or more memories, where at least one instruction is stored in the memory, and the instruction is loaded by the processor and executes the method.
Third embodiment
The present embodiment provides a computer-readable storage medium, in which at least one instruction is stored, and the instruction is loaded and executed by a processor to implement the method of the first embodiment. The computer readable storage medium may be, among others, ROM, random access memory, CD-ROM, magnetic tape, floppy disk, optical data storage device, and the like. The instructions stored therein may be loaded by a processor in the terminal and perform the above-described method.
Furthermore, it should be noted that the present invention may be provided as a method, apparatus or computer program product. Accordingly, embodiments of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the present invention may take the form of a computer program product embodied on one or more computer-usable storage media having computer-usable program code embodied in the medium.
Embodiments of the present invention are described with reference to flowchart illustrations and/or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, embedded processor, or other programmable data processing terminal to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing terminal to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing terminal to cause a series of operational steps to be performed on the computer or other programmable terminal to produce a computer implemented process such that the instructions which execute on the computer or other programmable terminal provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It should also 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 terminal 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 terminal. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or terminal that comprises the element.
Finally, it should be noted that while the above describes a preferred embodiment of the invention, it will be appreciated by those skilled in the art that, once the basic inventive concepts have been learned, numerous changes and modifications may be made without departing from the principles of the invention, which shall be deemed to be within the scope of the invention. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all such alterations and modifications as fall within the scope of the embodiments of the invention.
Claims (8)
1. A silicone oil filling simulation method for rhegmatogenous retinal detachment is characterized by comprising the following steps:
acquiring continuous medical scanning images of eyeballs, and constructing an eyeball three-dimensional model based on the medical scanning images;
based on the eyeball three-dimensional model and the liquid physical characteristics, carrying out stress simulation on the interaction boundary of the liquid and the eye cavity; wherein the liquid comprises water and silicone oil in the ocular lumen;
and analyzing the interphase coupling state of the silicone oil and the water in the eye cavity, expressing the interaction between fluids by using a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation by combining stress simulation results of the silicone oil and the interaction boundary of the water and the eye cavity to obtain the theoretically optimal silicone oil filling amount.
2. The silicone oil filling simulation method for pore-derived retinal detachment according to claim 1, wherein constructing a three-dimensional eyeball model based on the medical scan image comprises:
preprocessing the medical scanning image, wherein the preprocessing comprises smoothing denoising and gamma correction;
performing image segmentation on the preprocessed medical scanning image based on a preset region segmentation algorithm, performing pattern recognition on an eye structure, and completely describing a segmentation interval and eyeball tissues to obtain an eye structure image;
completing three-dimensional reconstruction of an eyeball tissue structure based on the eye structure image;
and calculating the volume of the communicated part of the reconstructed three-dimensional model, and removing the part of the communicated part with the volume smaller than a preset threshold value to finish the correction of the reconstructed three-dimensional model to obtain the eyeball three-dimensional model.
3. The silicone oil filling simulation method for rhegmatogenous retinal detachment of claim 1, wherein the force simulation of the interaction boundary of the liquid and the eye cavity is performed based on the eyeball three-dimensional model and the liquid physical characteristics, and comprises the following steps:
based on the physical characteristics of the liquid, the three-dimensional eyeball model is discretized into boundary particles, the boundary particles are added into the approximate calculation of the liquid density, and the stress simulation is carried out on the interaction boundary of the liquid and the intraocular cavity.
4. The silicone oil filling simulation method facing rhegmatogenous retinal detachment of claim 3, wherein the eyeball three-dimensional model is discretized into boundary particles based on the physical characteristics of the liquid, the boundary particles are added into the approximate calculation of the liquid density, and the interaction boundary of the liquid and the intraocular cavity is subjected to stress simulation, which comprises the following steps:
uniformly sampling the eyeball three-dimensional model by using a three-dimensional point cloud sampling algorithm;
giving the particle volume and mass to the sampling point;
the method comprises the steps of calculating the weight value of eyeball boundary particles borne by fluid particles in an eye cavity in physical calculation, considering the adhesion effect generated by interaction among different materials, calculating the pressure and viscosity of the interaction boundary of liquid and the eye cavity based on Newton's second law, and carrying out stress simulation on the interaction boundary of the liquid and the eye cavity.
5. The silicone oil filling simulation method for porogenic retinal detachment according to claim 1, wherein the method comprises the steps of performing interphase coupling state analysis on silicone oil and water in an eye cavity, expressing interaction between fluids by using a Lagrange mode, and generating visual simulation results of eyeball parameters under different silicone oil filling amounts in silicone oil filling by combining stress simulation results of interaction boundaries of the silicone oil and the water and the eye cavity to obtain the theoretically optimal silicone oil filling amount, wherein the method comprises the following steps:
performing single-phase fluid motion simulation on water and silicon oil in the eye cavity by using a smooth particle fluid dynamics method based on a Navier-Stokes equation;
the viscous force, the dragging force and the intersolubility of different flow phases are approximately calculated by using a kernel function, the interaction state of water and a silicon oil mixed liquid in a silicon oil filling state is simulated, and the silicon oil-water dynamic interaction calculation under the property of approximately incompressible non-dispersive fluid is carried out;
and (3) combining with an Euler dynamics method, continuously calculating the size of the surface tension direction of the space in a grid space based on a volume fraction method, and performing surface curvature analysis combining with the interphase surface tension of the two-phase flow to realize the simulation of the operation state under different silicone oil filling amounts, so as to generate the visual simulation result of eyeball parameters under different silicone oil filling amounts in the silicone oil filling operation, thereby obtaining the theoretically optimal silicone oil filling amount.
6. The method for simulating silicone oil filling for rhegmatogenous retinal detachment according to claim 5, wherein the silicone oil-hydrodynamics interaction calculation under the approximate incompressible non-dispersive fluid property comprises:
performing density field calculation on the fluid distribution state in the discrete time interval to obtain the current flow field compression state;
calculating the physical field quantity of each macroscopic fluid particle position by using a smooth particle fluid dynamics method according to the fluid description in a compression state;
and calculating the pressure field quantity at each moment by using a pressure force suppression fluid compression characteristic principle and by explicitly solving a pressure state equation or an implicit iteration method.
7. The method for simulating silicone oil filling for rhegmatogenous retinal detachment according to claim 5, wherein the surface curvature analysis combined with two-phase flow interphase surface tension comprises:
extracting a surface boundary by setting a threshold value by adopting a color field method, calculating the surface curvature of the particle, and acquiring a unit normal vector;
the shrinkage degree of the interface is controlled by the tension coefficients of different substances, the normal direction of the interface is the stress direction, and the tension is in direct proportion to the surface curvature;
simulating the intermolecular attraction and repulsion characteristics of the fluid, setting an action threshold value, enabling the particles to attract each other when the distance between the macro discrete particles is lower than the threshold value, and enabling the particles to repel each other when the distance between the particles exceeds the threshold value.
8. An electronic device comprising a processor and a memory; wherein the memory stores at least one instruction, and the instruction is loaded and executed by the processor to realize the silicone oil filling simulation method for rhegmatogenous retinal detachment according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110353522.XA CN113178001B (en) | 2021-04-01 | 2021-04-01 | Silicon oil filling simulation method for hole-origin retinal detachment and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110353522.XA CN113178001B (en) | 2021-04-01 | 2021-04-01 | Silicon oil filling simulation method for hole-origin retinal detachment and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113178001A true CN113178001A (en) | 2021-07-27 |
| CN113178001B CN113178001B (en) | 2023-07-04 |
Family
ID=76922923
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110353522.XA Active CN113178001B (en) | 2021-04-01 | 2021-04-01 | Silicon oil filling simulation method for hole-origin retinal detachment and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113178001B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110178784A1 (en) * | 2009-11-18 | 2011-07-21 | Daisuke Sato | Electro-oculography estimating device, electro-oculography calculating method, eye-gaze tracking device, wearable camera, head-mounted display, and electronic eyeglasses |
| CN102762220A (en) * | 2007-10-15 | 2012-10-31 | 森托科尔奥索生物科技公司 | Human anti-amyloid antibodies, compositions,methods and uses |
| CN107798198A (en) * | 2017-11-06 | 2018-03-13 | 北方工业大学 | Physical-based melting phenomenon realistic simulation method |
| CN109152657A (en) * | 2016-05-17 | 2019-01-04 | 诺华股份有限公司 | Automation viscous fluid control in vitrectomy |
-
2021
- 2021-04-01 CN CN202110353522.XA patent/CN113178001B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102762220A (en) * | 2007-10-15 | 2012-10-31 | 森托科尔奥索生物科技公司 | Human anti-amyloid antibodies, compositions,methods and uses |
| US20110178784A1 (en) * | 2009-11-18 | 2011-07-21 | Daisuke Sato | Electro-oculography estimating device, electro-oculography calculating method, eye-gaze tracking device, wearable camera, head-mounted display, and electronic eyeglasses |
| CN109152657A (en) * | 2016-05-17 | 2019-01-04 | 诺华股份有限公司 | Automation viscous fluid control in vitrectomy |
| CN107798198A (en) * | 2017-11-06 | 2018-03-13 | 北方工业大学 | Physical-based melting phenomenon realistic simulation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113178001B (en) | 2023-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Auer et al. | Reconstruction and finite element mesh generation of abdominal aortic aneurysms from computerized tomography angiography data with minimal user interactions | |
| CN104574503B (en) | contrast agent diffusion process simulation device | |
| Zhu et al. | Choroid neovascularization growth prediction with treatment based on reaction-diffusion model in 3-D OCT images | |
| CN106997594A (en) | The localization method and device of a kind of part tissue of eye | |
| Karimi et al. | Finite element modeling of the complex anisotropic mechanical behavior of the human sclera and pia mater | |
| CN114010157A (en) | Blood vessel plaque risk prediction method, device, equipment and medium | |
| Wang et al. | Finite element modeling of optic chiasmal compression | |
| Lee et al. | Evaluation of lacrimal tear drainage mechanism using dynamic fluoroscopic dacryocystography | |
| CN113178001B (en) | Silicon oil filling simulation method for hole-origin retinal detachment and electronic equipment | |
| Valencia et al. | Modeling of brain shift phenomenon for different craniotomies and solid models | |
| Yoshida et al. | A finite element analysis of the retinal hemorrhages accompanied by shaken baby syndrome/abusive head trauma | |
| Dai et al. | Simulating the effects of elevated intraocular pressure on ocular structures using a global finite element model of the human eye | |
| CN111047707B (en) | Bleeding simulation method of mixed particle blood model based on SPH | |
| Wei et al. | Fast ray-tracing of human eye optics on graphics processing units | |
| Szczerba et al. | Unstructured mesh generation from the Virtual Family models for whole body biomedical simulations | |
| Hötte et al. | Visualization of sliding and deformation of orbital fat during eye rotation | |
| Altintas | Effect of slice thickness variation on free vibration properties of micro-ct based trabecular bone models | |
| Chavoshnejad et al. | Mechanical hierarchy in the formation and modulation of cortical folding patterns | |
| Ledda et al. | Can wall shear-stress topology predict proliferative vitreoretinopathy localization following pars plana vitrectomy? | |
| Castro et al. | Hemodynamic differences in intracranial anerusysm blebs due to blood rheology | |
| Reilly et al. | Orbital fat swelling: A biomechanical theory and supporting model for spaceflight-associated neuro-ocular syndrome (SANS) | |
| Cebral et al. | Cerebral aneurysm hemodynamics modeling from 3d rotational angiography | |
| CN115019877B (en) | Method and device for modeling and updating biological tissue model and storage medium | |
| Kling et al. | Corneal Viscoelastic Properties from Finite-Element Analysis of | |
| JP2019518266A (en) | Biomechanical model generation for human or animal torso |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |