A three-dimensional dynamic biological simulation method based on formalized modeling

Abstract

The application discloses a three-dimensional dynamic biological simulation method based on formula modeling, and the method comprises the following steps: initializing a three-dimensional simulation space and initial positions, velocities and morphological parameters of a biological body, and constructing a parameterized biological morphological model; dynamically updating a two-dimensional morphological reference of the biological body and positions and postures of the biological body in a three-dimensional space based on morphological modulation parameters and motion velocities; for different structural parts of the biological body, an equal-angle slicing or random-direction generation equal-difference strategy is adopted to convert the two-dimensional morphological reference into a three-dimensional entity structure with symmetry or randomness, and linear interpolation is performed to enhance continuity; finally, the three-dimensional entity structure is rendered into volume data, superimposed intensity projection is performed along a preset direction, continuous frame images are generated, and the continuous frame images are output as multi-page TIFF format simulation data after a preset number of frames is reached.

Description

Technical Field

[0001] This invention relates to the field of computer graphics simulation technology, and in particular to a three-dimensional dynamic biological simulation method based on formulaic modeling. Background Technology

[0002] 3D vision will be the mainstream trend in the future field of vision. Currently, 3D vision is widely used in biological and medical research related to perception and decision-making. Since mainstream sensors are two-dimensional, while most application scenarios with research value are three-dimensional and dynamic, encoding the scene through optical design (light field imaging, superlenses, etc.) and then reconstructing it through computation is the main approach to dynamic 3D perception. However, computational reconstruction based on optimization is often slow and has low resolution, while reconstruction based on artificial intelligence methods suffers from the lack of 3D dynamic ground truth and difficulty in training.

[0003] Therefore, using computer simulation technology to generate realistic three-dimensional dynamic data plays an important role in algorithm verification and theoretical model testing. Summary of the Invention

[0004] The present invention aims to at least partially solve one of the technical problems in the related art.

[0005] This invention proposes a three-dimensional dynamic biological simulation method based on formulaic modeling to solve the problem of generating three-dimensional dynamic biological simulation data with complex morphological changes and free movement capabilities.

[0006] Another objective of this invention is to propose a three-dimensional dynamic biological simulation system based on formulaic modeling.

[0007] To achieve the above objectives, this invention proposes a three-dimensional dynamic biological simulation method based on formulaic modeling, comprising:

[0008] A two-dimensional morphological reference for an organism is constructed based on a pre-built parametric biological morphological model. The two-dimensional morphological reference and the position and posture of the organism in three-dimensional space are updated according to the morphological modulation parameters and the initial motion velocity vector in the biological morphological model.

[0009] The updated two-dimensional morphological benchmark is converted into a three-dimensional solid structure, wherein different three-dimensional construction strategies are used for different structural parts of the organism.

[0010] The generated 3D solid structure is rendered as 3D volume data, and intensity projection is superimposed along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data.

[0011] The three-dimensional dynamic biological simulation method based on formulaic modeling according to an embodiment of the present invention may also have the following additional technical features:

[0012] In one embodiment of the present invention, before constructing a two-dimensional morphological baseline of an organism based on a pre-built parametric biological morphology model, the method further includes:

[0013] Initialize the three-dimensional simulation space parameters, construct a parameterized biological morphology model, and set the initial center position coordinates, initial motion velocity vector, and total number of simulation frames for the parameterized biological model.

[0014] In one embodiment of the present invention, the three-dimensional simulation space parameters are initialized, a parameterized biological morphology model is constructed, and the initial center position coordinates, initial motion velocity vector, and total number of simulation frames of the parameterized biological model are set, including:

[0015] Create a three-dimensional Cartesian coordinate system of a preset size as the activity space of the organism;

[0016] Construct a parameterized biological morphology model, and initialize the initial center position coordinates of the organism as the center point of the coordinate system;

[0017] Set the initial motion velocity vector, the total number of simulation frames, the density parameters of the head umbrella-shaped body, and the density of the tentacles in the body.

[0018] In one embodiment of the present invention, a two-dimensional morphological reference for an organism is constructed based on a pre-built parametric biological morphological model. The two-dimensional morphological reference and the position and orientation of the organism in three-dimensional space are updated according to the morphological modulation parameters and initial motion velocity vector in the biological morphological model, including:

[0019] Set a time variable, and update the two-dimensional morphological reference coordinates according to the time variable to obtain a new two-dimensional morphological reference.

[0020] The organism's center position is updated based on the initial motion velocity vector, and the distances from the organism's center position to the three-dimensional space boundary are determined in the x, y, and z directions respectively. If the distance is less than the preset distance, a boundary collision is detected, and the motion velocity vector is reversed in the corresponding direction.

[0021] Calculate the target direction vector of the head joint nodes, normalize the target direction vector, and update the matrix of the joint node direction vector and the matrix of the joint node position using a smooth propagation mechanism.

[0022] In one embodiment of the present invention, the updated two-dimensional morphological reference is converted into a three-dimensional solid structure, wherein different structural parts of the organism employ different three-dimensional construction strategies, including:

[0023] For structural parts requiring symmetry, a complete three-dimensional symmetrical structure is generated using an equal-angle slicing method;

[0024] For structural parts with random or disordered characteristics, three-dimensional modeling is performed by generating discrete filamentary structures in random directions;

[0025] Based on the three-dimensional solid structure of organisms, the three-dimensional coordinates of the three-dimensional solid structure are rotated and translated.

[0026] Linear interpolation is performed on the spatially discrete parts to generate intermediate connection points between adjacent structural points.

[0027] In one embodiment of the present invention, the generated three-dimensional solid structure is rendered as three-dimensional volume data, and superimposed intensity projection is performed along a preset direction to generate frame data. When a preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data, including:

[0028] Initialize a three-dimensional volume data matrix with dimensions consistent with the three-dimensional simulation space;

[0029] The components of an organism are intensified, with the intensity decreasing gradually from the center to the edge.

[0030] The three-dimensional volume data matrix is ​​subjected to intensity superposition projection along a preset direction to generate a projected image, and the projected images of consecutive frames are stored in a frame sequence array;

[0031] When the number of frames in the frame sequence array reaches the preset total number of frames, all frame data are merged and output as multi-page TIFF image simulation data.

[0032] To achieve the above objectives, another aspect of the present invention proposes a three-dimensional dynamic biological simulation system based on formulaic modeling, comprising:

[0033] An organism update module is used to construct a two-dimensional morphological reference of an organism based on the parameterized biological morphology model. The biological morphology model includes morphological modulation parameters. The module updates the two-dimensional morphological reference of the organism and the position and posture of the organism in three-dimensional space according to the morphological modulation parameters and the initial motion velocity vector.

[0034] The three-dimensional structure generation module is used to convert the updated two-dimensional morphological benchmark into a three-dimensional solid structure, and different three-dimensional construction strategies are used for different structural parts of the organism.

[0035] The simulation image output module is used to render the generated 3D solid structure into 3D volume data, and perform superimposed intensity projection along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data. Attached Figure Description

[0036] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein:

[0037] Figure 1 This is a flowchart of a three-dimensional dynamic biological simulation method based on formulaic modeling according to an embodiment of the present invention;

[0038] Figure 2 It is a two-dimensional dynamic simulation image of a jellyfish organism according to an embodiment of the present invention;

[0039] Figure 3 It is a rendered three-dimensional simulation image according to an embodiment of the present invention;

[0040] Figure 4 This is an architectural diagram of the technical solution according to an embodiment of the present invention;

[0041] Figure 5 This is a structural diagram of a three-dimensional dynamic biological simulation system based on formulaic modeling according to an embodiment of the present invention. Detailed Implementation

[0042] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other. The present invention will now be described in detail with reference to the accompanying drawings and embodiments.

[0043] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0044] The following description, with reference to the accompanying drawings, describes a three-dimensional dynamic biological simulation method and system based on formulaic modeling, according to an embodiment of the present invention.

[0045] Figure 1 This is a flowchart of a three-dimensional dynamic biological simulation method based on formulaic modeling according to an embodiment of the present invention, such as... Figure 1 As shown, the method includes:

[0046] S1. Construct a two-dimensional morphological reference for the organism based on a pre-built parametric biological morphological model. Update the two-dimensional morphological reference and the position and posture of the organism in three-dimensional space according to the morphological modulation parameters and initial motion velocity vector in the biological morphological model.

[0047] S2, convert the updated two-dimensional morphological benchmark into a three-dimensional solid structure, wherein different structural parts of the organism adopt different three-dimensional construction strategies;

[0048] S3 renders the generated 3D solid structure into 3D volume data and performs superimposed intensity projection along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data.

[0049] In one embodiment of the present invention, before constructing a two-dimensional morphological baseline of an organism based on a pre-built parametric biological morphology model, the method further includes:

[0050] Initialize the three-dimensional simulation space parameters, construct a parameterized biological morphology model, and set the initial center position coordinates, initial motion velocity vector, and total number of simulation frames for the parameterized biological model.

[0051] In one embodiment of the present invention, the three-dimensional simulation space parameters are initialized, a parameterized biological morphology model is constructed, and the initial center position coordinates, initial motion velocity vector, and total number of simulation frames of the parameterized biological model are set, including:

[0052] Create a three-dimensional Cartesian coordinate system of a preset size as the activity space of the organism;

[0053] Construct a parameterized biological morphology model, and initialize the initial center position coordinates of the organism as the center point of the coordinate system;

[0054] Set the initial motion velocity vector, the total number of simulation frames, the density parameters of the head umbrella-shaped body, and the density of the tentacles in the body.

[0055] In one embodiment of the present invention, a two-dimensional morphological reference for an organism is constructed based on a pre-built parametric biological morphological model. The two-dimensional morphological reference and the position and orientation of the organism in three-dimensional space are updated according to the morphological modulation parameters and initial motion velocity vector in the biological morphological model, including:

[0056] Set a time variable, and update the two-dimensional morphological reference coordinates according to the time variable to obtain a new two-dimensional morphological reference.

[0057] The organism's center position is updated based on the initial motion velocity vector, and the distances from the organism's center position to the three-dimensional space boundary are determined in the x, y, and z directions respectively. If the distance is less than the preset distance, a boundary collision is detected, and the motion velocity vector is reversed in the corresponding direction.

[0058] Calculate the target direction vector of the head joint nodes, normalize the target direction vector, and update the matrix of the joint node direction vector and the matrix of the joint node position using a smooth propagation mechanism.

[0059] In one embodiment of the present invention, the updated two-dimensional morphological reference is converted into a three-dimensional solid structure, wherein different structural parts of the organism employ different three-dimensional construction strategies, including:

[0060] For structural parts requiring symmetry, a complete three-dimensional symmetrical structure is generated using an equal-angle slicing method;

[0061] For structural parts with random or disordered characteristics, three-dimensional modeling is performed by generating discrete filamentary structures in random directions;

[0062] Based on the three-dimensional solid structure of organisms, the three-dimensional coordinates of the three-dimensional solid structure are rotated and translated.

[0063] Linear interpolation is performed on the spatially discrete parts to generate intermediate connection points between adjacent structural points.

[0064] In one embodiment of the present invention, the generated three-dimensional solid structure is rendered as three-dimensional volume data, and superimposed intensity projection is performed along a preset direction to generate frame data. When a preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data, including:

[0065] Initialize a three-dimensional volume data matrix with dimensions consistent with the three-dimensional simulation space;

[0066] The components of an organism are intensified, with the intensity decreasing gradually from the center to the edge.

[0067] The three-dimensional volume data matrix is ​​subjected to intensity superposition projection along a preset direction to generate a projected image, and the projected images of consecutive frames are stored in a frame sequence array;

[0068] When the number of frames in the frame sequence array reaches the preset total number of frames, all frame data are merged and output as multi-page TIFF image simulation data.

[0069] This invention provides a novel three-dimensional dynamic biological simulation method. This invention uses a parameterized mathematical model and a temporal dynamic update mechanism to generate biological simulation data with autonomous movement and dynamic changes in body shape (such as contraction and relaxation) within a preset three-dimensional space.

[0070] Furthermore, taking the simulation of a three-dimensional jellyfish as an example, the technical solution of this invention is as follows: Figure 2 , Figure 3 and Figure 4 As shown, the main steps include:

[0071] S10, Initialize simulation environment parameters;

[0072] S20, constructing a parameterized biological morphology model;

[0073] S30 enables dynamic movement updates of the organism's body structure;

[0074] S40 enables the transformation of biological structures from two-dimensional forms to three-dimensional entities;

[0075] S50 performs linear interpolation on jellyfish organisms to make their morphological structure more continuous;

[0076] S60 generates three-dimensional volume data based on simulation results, calculates the intensity distribution of different parts of the organism, and generates a superimposed intensity projection image along the z-axis;

[0077] S70, repeat steps S30 to S60, and output simulation data in the form of multi-page TIFF images.

[0078] Furthermore, step S1 specifically includes the following steps:

[0079] Create a 3D Cartesian coordinate system of a preset size (e.g., 1600×1600×800) as the organism's activity space, and initialize the organism's initial center position coordinates pos as the center point of the coordinate system. Set the initial motion velocity vector vel, the total number of frames num_frames, the density parameter N_slices_head of the head umbrella-shaped body, and the density N_slices_body of the body tentacles.

[0080] Furthermore, step 20 specifically includes the following steps:

[0081] S201. Generate the basic morphological parameter sequence of the organism, define the morphological index sequence vector i, and calculate the key morphological parameters k, e, and d, where k is the lateral fluctuation parameter, e is the longitudinal distribution parameter, and d is the morphological modulation parameter controlling the dynamic change amplitude of the organism. The specific calculation formula is as follows:

[0082]

[0083]

[0084]

[0085] S202. Calculate the two-dimensional morphological reference coordinates X and Y. These coordinates determine the outline of the jellyfish's bell-shaped body and tentacles. The specific calculation formula is as follows: (where scl is the scaling factor and t is the motion time variable).

[0086]

[0087]

[0088] S203. Define the joint node index sequence vector mi, initialize the joint node position matrix M and direction vector matrix T, and set the initial direction vectors to all [0, 1, 0]; calculate the reference position M0 of each joint node, where the x-coordinate of the reference position M0 is the mean of the morphological reference X-coordinate, and the y-coordinate is the maximum value of the morphological reference Y-coordinate in each joint segment. Calculate the distance matrix D between adjacent joint nodes, which is used to constrain the connection relationship between joints.

[0089] Furthermore, step S30 specifically includes the following steps:

[0090] S301. Update the time variable t, recalculate the two-dimensional morphological reference coordinates X and Y, obtain a new two-dimensional morphological reference, and realize the periodic dynamic change of the organism's morphology.

[0091] S302. Update the organism's center position pos = pos + vel based on the velocity vector. Determine the distance from the organism's center position to the 3D space boundary in the x, y, and z directions respectively. If the distance is too small, a boundary collision is detected, and the motion velocity vel is reversed in the corresponding direction.

[0092] S303. Calculate the target direction vector dir = pos - M(1,:) for the head joint nodes and normalize the direction vector. Update the matrix T representing the direction vectors of all joint nodes and the matrix M representing the positions of all joint nodes using a smooth propagation mechanism. The head node first performs smooth following, and each subsequent node performs smooth following relative to the previous node, thereby constraining the overall body shape of the organism.

[0093] Furthermore, step S40 specifically includes the following steps:

[0094] S401. Different 3D processing is performed according to the joint segment type: For the umbrella-shaped part of the head, a complete circumferential 3D model is created using an isometric slicing method to construct a complete and symmetrical 3D umbrella-shaped structure. The isometric slicing formula is as follows:

[0095]

[0096] For touch methods, a set of random angles is generated, and the specific generation formula is as follows:

[0097]

[0098] Discrete, filamentous tentacle structures are generated along these random directions, simulating the natural and disordered morphological characteristics of biological tentacles;

[0099] S402. For each generated pixel, calculate its 3D base coordinates and subtract the reference coordinates to obtain its relative coordinates. Then, based on the current motion posture of the organism, calculate the rotation matrix from the direction vector in the current relative coordinate system to the target vector in the world coordinate system. Apply the rotation matrix to each relative coordinate to perform coordinate transformation, obtain the rotated coordinates, and finally add them to the current joint node position coordinates to obtain the final 3D coordinates.

[0100] Furthermore, step S50 specifically includes the following steps:

[0101] For all relatively discrete points that require linear interpolation, the following linear interpolation method is applied: Taking p1 and p2 as an example, first calculate the distance between the two points, as shown in the following formula:

[0102]

[0103] After obtaining the distance, determine the number of interpolation steps, perform linear interpolation between the two points, and generate intermediate connection points. The specific formula is as follows:

[0104]

[0105] Where i_step represents the number of steps in the interpolation calculation between points p1 and p2. Finally, all interpolated points are added to the 3D point cloud data.

[0106] Furthermore, step S60 specifically includes the following steps:

[0107] S601. Initialize the three-dimensional volume data matrix V, with a size consistent with the simulation space dimensions;

[0108] S602. Intensity fill is applied to the 3D model of the head segment, setting the head intensity to a fixed value. Intensity fill is applied to the 3D model of the touch device, calculating the intensity change based on the radial position to achieve an intensity gradient distribution that decreases from the center to the edge.

[0109] S603. Perform intensity superposition projection on the three-dimensional volume data matrix V along the z-axis, store the projection result in the frame sequence array, and update the save counter;

[0110] Furthermore, step S70 specifically includes the following steps:

[0111] Repeat steps S30 to S60 to display the dynamic, continuous free movement of 3D jellyfish in space in the 3D display window. When the program reaches the preset total number of frames num_frames, call the writeMultipageTiff function to merge all frame data and output it as a multi-page TIFF file.

[0112] To achieve the above embodiments, such as Figure 5As shown, this embodiment also provides a three-dimensional dynamic biological simulation system 10 based on formulaic modeling, including:

[0113] The organism update module 100 is used to construct a two-dimensional morphological reference of an organism based on a parameterized biological morphology model. The biological morphology model includes morphological modulation parameters. Based on the morphological modulation parameters and the initial motion velocity vector, the two-dimensional morphological reference of the organism and the position and posture of the organism in three-dimensional space are updated.

[0114] The three-dimensional structure generation module 200 is used to convert the updated two-dimensional morphological benchmark into a three-dimensional solid structure, and different three-dimensional construction strategies are used for different structural parts of the organism.

[0115] The simulation image output module 300 is used to render the generated three-dimensional solid structure into three-dimensional volume data, and perform superimposed intensity projection along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data.

[0116] Furthermore, the organism renewal module 100 is also used for:

[0117] Set a time variable, and update the two-dimensional morphological reference coordinates according to the time variable to obtain a new two-dimensional morphological reference.

[0118] The organism's center position is updated based on the initial motion velocity vector, and the distances from the organism's center position to the three-dimensional space boundary are determined in the x, y, and z directions respectively. If the distance is less than the preset distance, a boundary collision is detected, and the motion velocity vector is reversed in the corresponding direction.

[0119] Calculate the target direction vector of the head joint nodes, normalize the target direction vector, and update the matrix of the joint node direction vector and the matrix of the joint node position using a smooth propagation mechanism.

[0120] Furthermore, the three-dimensional structure generation module 200 is also used for:

[0121] For structural parts requiring symmetry, a complete three-dimensional symmetrical structure is generated using an equal-angle slicing method;

[0122] For structural parts with random or disordered characteristics, three-dimensional modeling is performed by generating discrete filamentary structures in random directions;

[0123] Based on the three-dimensional solid structure of organisms, the three-dimensional coordinates of the three-dimensional solid structure are rotated and translated.

[0124] Linear interpolation is performed on the spatially discrete parts to generate intermediate connection points between adjacent structural points.

[0125] Furthermore, the simulation image output module 300 is also used for:

[0126] Initialize a three-dimensional volume data matrix with dimensions consistent with the three-dimensional simulation space;

[0127] The components of an organism are intensified, with the intensity decreasing gradually from the center to the edge.

[0128] The three-dimensional volume data matrix is ​​subjected to intensity superposition projection along a preset direction to generate a projected image, and the projected images of consecutive frames are stored in a frame sequence array;

[0129] When the number of frames in the frame sequence array reaches the preset total number of frames, all frame data are merged and output as multi-page TIFF image simulation data.

[0130] In summary, this invention achieves dynamic simulation of complex biological morphology: the simulation data generated by this invention is no longer limited to simple geometric shapes, but can simulate organisms with complex body structures (such as the umbrella-shaped body and tentacles of a jellyfish), and can dynamically display their morphological changes such as contraction and relaxation, thus enhancing the dynamic complexity of the simulation. It also achieves simulation of macroscopic biological behavior: this invention is no longer limited to static anatomical structures, but achieves simulation of macroscopic autonomous movement and non-rigid deformation of organisms, providing a new, high-quality simulation data source for research in fields such as biodynamics and behavioral science. Furthermore, it achieves three-dimensional dynamic biological simulation: the simulation data generated by this invention is no longer limited to a two-dimensional plane, but can reproduce the three-dimensional morphology and structural characteristics of organisms in three-dimensional space, providing higher-quality simulation data that can be directly applied to algorithm training and verification in fields such as three-dimensional perception.

[0131] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of the present invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0132] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise explicitly specified.

Claims

1. A three-dimensional dynamic biological simulation method based on formulaic modeling, characterized in that, include: A two-dimensional morphological reference for an organism is constructed based on a pre-built parametric biological morphological model. The two-dimensional morphological reference and the position and posture of the organism in three-dimensional space are updated according to the morphological modulation parameters and the initial motion velocity vector in the biological morphological model. The updated two-dimensional morphological benchmark is converted into a three-dimensional solid structure, wherein different three-dimensional construction strategies are used for different structural parts of the organism. The generated 3D solid structure is rendered as 3D volume data, and intensity projection is superimposed along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data. The updated two-dimensional morphological baseline is converted into a three-dimensional solid structure, wherein different structural parts of the organism employ different three-dimensional construction strategies, including: For structural parts requiring symmetry, a complete three-dimensional symmetrical structure is generated using an equal-angle slicing method; For structural parts with random or disordered characteristics, three-dimensional modeling is performed by generating discrete filamentary structures in random directions; Based on the three-dimensional solid structure of organisms, the three-dimensional coordinates of the three-dimensional solid structure are rotated and translated. Linear interpolation is performed on the spatially discrete parts to generate intermediate connection points between adjacent structural points; The generated 3D solid structure is rendered as 3D volume data, and intensity projection is superimposed along a preset direction to generate frame data. When the total number of simulation frames is reached, the frame data is merged and output as multi-page TIFF image simulation data, including: Initialize a three-dimensional volume data matrix with dimensions consistent with the three-dimensional simulation space; The components of an organism are intensified, with the intensity decreasing gradually from the center to the edge. The three-dimensional volume data matrix is ​​subjected to intensity superposition projection along a preset direction to generate a projected image, and the projected images of consecutive frames are stored in a frame sequence array; When the number of frames in the frame sequence array reaches the total number of simulation frames, all frame data are merged and output as multi-page TIFF image simulation data.

2. The method according to claim 1, characterized in that, Before constructing a two-dimensional morphological benchmark for an organism based on a pre-built parametric biological morphology model, the following steps are also included: Initialize the three-dimensional simulation space parameters, construct a parameterized biological morphology model, and set the initial center position coordinates, initial motion velocity vector, and total number of simulation frames for the parameterized biological model.

3. The method according to claim 2, characterized in that, Initialize the 3D simulation space parameters, construct a parameterized biological morphology model, and set the initial center position coordinates, initial motion velocity vector, and total number of simulation frames for the parameterized biological model, including: Create a three-dimensional Cartesian coordinate system of a preset size as the activity space of the organism; Construct a parameterized biological morphology model, and initialize the initial center position coordinates of the organism as the center point of the coordinate system; Set the initial motion velocity vector, the total number of simulation frames, the density parameters of the head umbrella-shaped body, and the density of the tentacles in the body.

4. The method according to claim 1, characterized in that, A two-dimensional morphological baseline of an organism is constructed based on a pre-built parametric biological morphological model. The two-dimensional morphological baseline and the organism's position and orientation in three-dimensional space are updated according to the initial center position coordinates and initial velocity vector in the biological morphological model, including: Set a time variable, and update the two-dimensional morphological reference coordinates according to the time variable to obtain a new two-dimensional morphological reference. The organism's center position is updated based on the initial center position coordinates and the initial motion velocity vector. The distances from the organism's center position to the three-dimensional space boundary in the x, y, and z directions are determined respectively. If the distance is less than a preset distance, a boundary collision is detected, and the motion velocity vector is reversed in the corresponding direction. Calculate the target direction vector of the head joint nodes, normalize the target direction vector, and update the matrix of the joint node direction vector and the matrix of the joint node position using a smooth propagation mechanism.

5. A three-dimensional dynamic biological simulation system based on formulaic modeling, characterized in that, include: The organism update module is used to construct a two-dimensional morphological reference of the organism based on a pre-built parametric biological morphological model, and update the two-dimensional morphological reference of the organism and the position and posture of the organism in three-dimensional space according to the morphological modulation parameters and the initial motion velocity vector in the biological morphological model. The three-dimensional structure generation module is used to convert the updated two-dimensional morphological benchmark into a three-dimensional solid structure, wherein different three-dimensional construction strategies are used for different structural parts of the organism. The simulation image output module is used to render the generated three-dimensional solid structure into three-dimensional volume data, and to perform superimposed intensity projection along a preset direction to generate frame data. When the preset total number of frames is reached, the frame data is merged and output as multi-page TIFF image simulation data. The 3D structure generation module is also used to generate complete 3D symmetrical structures for structural parts with symmetry requirements by using equal-angle slicing. For structural parts with random or disordered characteristics, three-dimensional modeling is performed by generating discrete filamentary structures in random directions; Based on the three-dimensional solid structure of organisms, the three-dimensional coordinates of the three-dimensional solid structure are rotated and translated. Linear interpolation is performed on the spatially discrete parts to generate intermediate connection points between adjacent structural points; The simulation image output module is also used to: initialize a three-dimensional volume data matrix with dimensions consistent with the three-dimensional simulation space; The components of an organism are intensified, with the intensity decreasing gradually from the center to the edge. The three-dimensional volume data matrix is ​​subjected to intensity superposition projection along a preset direction to generate a projected image, and the projected images of consecutive frames are stored in a frame sequence array; When the number of frames in the frame sequence array reaches the total number of simulation frames, all frame data are merged and output as multi-page TIFF image simulation data.

6. The system according to claim 5, characterized in that, The organism update module is also used for: Set a time variable, and update the two-dimensional morphological reference coordinates according to the time variable to obtain a new two-dimensional morphological reference. The organism's center position is updated based on the initial motion velocity vector, and the distances from the organism's center position to the three-dimensional space boundary are determined in the x, y, and z directions respectively. If the distance is less than the preset distance, a boundary collision is detected, and the motion velocity vector is reversed in the corresponding direction. Calculate the target direction vector of the head joint nodes, normalize the target direction vector, and update the matrix of the joint node direction vector and the matrix of the joint node position using a smooth propagation mechanism.

Images