Cell motion trajectory curve generation method, device, equipment and storage medium

Abstract

This invention provides a method, apparatus, device, and storage medium for generating cell motion trajectory curves. The method includes: acquiring multiple diffraction ring images of microspheres engulfed by cells; determining multiple three-dimensional coordinates of the microspheres based on the multiple diffraction ring images; determining the multiple three-dimensional coordinates of the microspheres as multiple three-dimensional coordinates of the cell; and generating a three-dimensional motion trajectory curve of the cell based on the multiple three-dimensional coordinates of the cell. This allows for the acquisition of multiple diffraction ring images of microspheres engulfed by cells and the automatic determination of multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images. Because the cell engulfs the microspheres, the microspheres are subject to dual constraints from the cell's cytoskeleton and cell membrane, allowing the microspheres to effectively track and reproduce changes in cell motion. The three-dimensional motion trajectory curve of the cell more realistically reflects changes in cell motion in three-dimensional space, increasing the realism of the cell motion changes reflected by the three-dimensional motion trajectory curve.

Description

Technical Field

[0001] This invention relates to the field of biomedical materials research technology, and in particular to a method, apparatus, device and storage medium for generating cell movement trajectory curves. Background Technology

[0002] Currently, the study of cell motility changes in vitro mainly relies on optical detection methods. In optical detection, a traction microscope is used to calibrate the displacement information of fluorescent beads in contact with the cell membrane, and then the cell motility changes in response to different drugs are determined based on this displacement information. This is equivalent to indirectly determining cell motility changes based on the displacement information of the fluorescent beads.

[0003] However, due to the limitations of traction microscopes, displacement information can only reflect changes in cell movement in two-dimensional space, making it impossible for optical detection, an indirect measurement method, to accurately reflect changes in cell movement in vitro. Therefore, how to accurately reflect changes in cell movement is a problem that urgently needs to be solved. Summary of the Invention

[0004] This invention provides a method, apparatus, device, and storage medium for generating cell motion trajectory curves, which can solve the problem that optical detection, an indirect measurement method, cannot truly reflect the changes in cell motion in vitro.

[0005] According to a first aspect of the present invention, a method for generating a cell motion trajectory curve is provided, the method comprising:

[0006] Acquire multiple diffraction ring images of microspheres engulfed by cells;

[0007] Based on the multiple diffraction ring images, the multiple three-dimensional coordinates of the microsphere are determined;

[0008] The multiple three-dimensional coordinates of the microspheres are determined as the multiple three-dimensional coordinates of the cell;

[0009] The three-dimensional motion trajectory curve of the cell is generated based on multiple three-dimensional coordinates of the cell.

[0010] According to a second aspect of the present invention, a cell motion trajectory curve generation device is provided, the device comprising:

[0011] The image acquisition module is used to acquire multiple diffraction ring images of microspheres engulfed by cells;

[0012] The first determining module is used to determine multiple three-dimensional coordinates of the microsphere based on the multiple diffraction ring images;

[0013] The second determining module is used to determine the multiple three-dimensional coordinates of the microsphere as multiple three-dimensional coordinates of the cell;

[0014] The curve generation module is used to generate a three-dimensional motion trajectory curve of the cell based on multiple three-dimensional coordinates of the cell.

[0015] According to a third aspect of the present invention, an electronic device is provided, comprising a processor and a memory.

[0016] The memory is used to store code and related data;

[0017] The processor is used to execute code in the memory to implement the cell motion trajectory curve generation method as described in any of the embodiments of the present invention.

[0018] According to a fourth aspect of the present invention, a storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the method for generating cell motion trajectory curves as described in any of the embodiments of the present invention.

[0019] The present invention has the following beneficial effects and advantages:

[0020] In this embodiment of the invention, multiple diffraction ring images of microspheres engulfed by cells are acquired; multiple three-dimensional coordinates of the microspheres are determined based on the multiple diffraction ring images; the multiple three-dimensional coordinates of the microspheres are then used as multiple three-dimensional coordinates of the cells; and a three-dimensional motion trajectory curve of the cells is generated based on the multiple three-dimensional coordinates of the cells. That is, the technical solution of this invention can acquire multiple diffraction ring images of microspheres engulfed by cells and automatically determine multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images. Because the cells engulf the microspheres, the microspheres are subject to the dual constraints of the cell's cytoskeleton and cell membrane, allowing the microspheres to effectively track and reproduce the cell's motion changes. Therefore, the multiple three-dimensional coordinates of the microspheres can be used as multiple three-dimensional coordinates of the cells. Finally, a three-dimensional motion trajectory curve of the cells is generated based on the multiple three-dimensional coordinates of the cells. This three-dimensional motion trajectory curve more realistically reflects the cell's motion changes in three-dimensional space, increasing the realism of the cell's motion changes reflected by the three-dimensional motion trajectory curve, and achieving the goal of realistically reflecting the cell's motion changes in three-dimensional space. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a schematic flowchart of a method for generating cell motion trajectory curves provided in an embodiment of the present invention;

[0023] Figure 2This is a schematic diagram of the microsphere measuring device provided in an embodiment of the present invention;

[0024] Figure 3 This is another flowchart illustrating the method for generating cell motion trajectory curves provided in this embodiment of the invention;

[0025] Figure 4 This is a schematic diagram of the pixel points obtained from the initial two-dimensional coordinates of cell expansion provided in an embodiment of the present invention;

[0026] Figure 5 This is a schematic diagram of a target light intensity curve provided in an embodiment of the present invention;

[0027] Figure 6 This is a schematic diagram of a cell motion trajectory curve generation device provided in an embodiment of the present invention;

[0028] Figure 7 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0029] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. 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 are within the scope of protection of the present invention.

[0030] The terms “comprising” and “having”, and any variations thereof, in the specification, claims, and accompanying drawings of this invention are intended to cover a non-exclusive inclusion, for example, a process, method, system, product, or apparatus that includes a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such process, method, product, or apparatus.

[0031] The technical solution of the present invention will be described in detail below with reference to specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments.

[0032] Figure 1 This is a schematic flowchart illustrating a method for generating cell motion trajectory curves according to an embodiment of the present invention. This method can be executed by a cell motion trajectory curve generation device, which can be implemented using software and / or hardware. In a specific embodiment, the device can be integrated into an electronic device, such as a computer or server. The following embodiments will illustrate this using the integration of the device into an electronic device as an example. Figure 1The method may specifically include the following steps:

[0033] Step 101: Obtain multiple diffraction ring images of the microspheres engulfed by cells.

[0034] In this context, microspheres can be understood as particulate dispersion systems formed by drug dispersion or adsorption within a polymer matrix. Diffraction rings can be understood as bright and dark ring-shaped fringes with specific spacing produced by microspheres under laser illumination. Diffraction ring images can be understood as images of bright and dark ring-shaped fringes with specific spacing produced by microspheres at different locations under laser illumination.

[0035] In one alternative implementation, the cells that have engulfed the microspheres can be placed in a culture dish, and then the culture dish can be placed in... Figure 2 Below the light source, through Figure 2 When a light source illuminates cells, the microspheres at different distances produce different diffraction rings due to interference effects, which can be utilized... Figure 2 The camera lens captures diffraction rings generated by microspheres at different positions within cells in a culture dish, obtaining optical videos of a preset duration. These videos include videos of cardiomyocytes and diffraction rings. The optical videos are acquired from the camera via data transmission. Then, the diffraction ring videos are segmented to obtain multiple diffraction ring images of the microspheres at multiple time points. The three-dimensional coordinates of the microspheres and cells determined from these multiple diffraction ring images are all at the nanometer level, enabling nanometer-scale microsphere measurement. Finally, a three-dimensional motion trajectory curve of the cell is generated based on the multiple three-dimensional coordinates of the cell. This three-dimensional motion trajectory curve more realistically reflects the cell's motion changes in three-dimensional space, increasing the realism and accuracy of the cell motion changes reflected by the trajectory curve, achieving a high-precision and realistic reflection of cell motion changes in three-dimensional space.

[0036] Step 102: Determine multiple three-dimensional coordinates of the microsphere based on multiple diffraction ring images.

[0037] The diffraction ring images are numbered, and the earlier the diffraction ring image was captured, the smaller the image number.

[0038] In one optional implementation, the diffraction ring image corresponding to the smallest image number among multiple diffraction ring images can be determined as the current diffraction ring image; the initial two-dimensional center coordinates of the current diffraction ring image are obtained; the pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates are expanded to obtain multiple pixel points; the current diffraction ring image is segmented with the multiple pixel points as the segmentation origin to obtain multiple sets of segmented images; the target two-dimensional center coordinates of the current diffraction ring image are determined based on the multiple sets of segmented images; the three-dimensional coordinates of the microsphere at the corresponding time point of the current diffraction ring image are determined based on the target two-dimensional center coordinates and the current diffraction ring image; and the three-dimensional coordinates of the microsphere at the corresponding time point of the current diffraction ring image are determined based on the target two-dimensional center coordinates and the current diffraction ring image. When the image number equals a preset threshold, the step of determining multiple three-dimensional coordinates of the microspheres as multiple three-dimensional coordinates of the cell is triggered. When the image number of the current diffraction ring image is less than the preset threshold, the target two-dimensional center coordinates of the current diffraction ring image are determined as the initial two-dimensional center coordinates of the diffraction ring image corresponding to the next image number of the current diffraction ring image. The diffraction ring image corresponding to the next image number of the current diffraction ring image is determined as the current diffraction ring image, triggering the step of expanding the diffraction ring corresponding to the initial two-dimensional center coordinates with the initial two-dimensional center coordinates as the expansion origin to obtain an expanded image. Here, the preset threshold can be understood as the preset number of multiple diffraction ring images.

[0039] Specifically, a coordinate acquisition interface can be displayed so that the user can input the initial two-dimensional center coordinates on the interface, and then the initial two-dimensional center coordinates of the first diffraction ring image input by the user can be obtained.

[0040] For example, multiple diffraction ring images include P1, P2, and P3, where P1 was captured at time T1 and its image number is 1, P2 was captured at time T2 and its image number is 2, P3 was captured at time T3 and its image number is 3, T1 < T2 < T3, and the preset threshold is 3. The diffraction ring image corresponding to the smallest image number (1) among the multiple diffraction ring images is determined as the current diffraction ring image (P1), and the initial two-dimensional center coordinates of the current diffraction ring image (P1) are obtained; the pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates are expanded to obtain multiple pixel points; the current diffraction ring image is segmented with the multiple pixel points as the segmentation origin to obtain multiple sets of segmented images; based on the multiple sets of segmented images, the target two-dimensional center coordinates (X1, Y1) of the current diffraction ring image (P1) are determined; based on the target two-dimensional center coordinates and the current diffraction ring image, the three-dimensional coordinates of the microsphere at the corresponding time point of the current diffraction ring image (P1) are determined.

[0041] If the image number (1) of the current diffraction ring image is less than the preset threshold (3), the target two-dimensional center coordinates of the current diffraction ring image (P1) are determined as the initial two-dimensional center coordinates of the diffraction ring image (P2) corresponding to the next image number of the current diffraction ring image. The diffraction ring image corresponding to the next image number of the current diffraction ring image is determined as the current diffraction ring image (P2). The pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates are expanded to obtain multiple pixels. The current diffraction ring image is divided with the multiple pixels as the division origin to obtain multiple sets of segmented images. Based on the multiple sets of segmented images, the target two-dimensional center coordinates (X2, Y2) of the current diffraction ring image (P2) are determined. Based on the target two-dimensional center coordinates and the current diffraction ring image, the three-dimensional coordinates of the microsphere at the corresponding time point of the current diffraction ring image (P2) are determined.

[0042] If the image number (2) of the current diffraction ring image is less than the preset threshold (3), the target two-dimensional center coordinates of the current diffraction ring image (P2) are determined as the initial two-dimensional center coordinates of the diffraction ring image corresponding to the next image number of the current diffraction ring image. The diffraction ring image corresponding to the next image number of the current diffraction ring image is determined as the current diffraction ring image (P3). The pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates are expanded to obtain multiple pixels. The current diffraction ring image is segmented with the multiple pixels as the segmentation origin to obtain multiple sets of segmented images. Based on the multiple sets of segmented images, the target two-dimensional center coordinates (X3, Y3) of the current diffraction ring image (P3) are determined. Based on the target two-dimensional center coordinates and the current diffraction ring image, the three-dimensional coordinates of the microsphere at the corresponding time point of the current diffraction ring image (P3) are determined. If the image number (3) of the current diffraction ring image is equal to the preset threshold (3), step 103 is executed.

[0043] Step 103: Determine the multiple three-dimensional coordinates of the microspheres as multiple three-dimensional coordinates of the cell.

[0044] In this embodiment of the invention, since the cell engulfs the microsphere, the microsphere is doubly constrained by the cell's cytoskeleton and cell membrane. The microsphere can effectively track and reproduce the cell's movement changes, and the changes in the three-dimensional coordinates of the microsphere reflect the changes in the three-dimensional coordinates of the cell. This can more realistically reflect the cell's movement changes in three-dimensional space. Therefore, multiple three-dimensional coordinates of the microsphere can be determined as multiple three-dimensional coordinates of the cell. Finally, a three-dimensional motion trajectory curve of the cell is generated based on the multiple three-dimensional coordinates of the cell. The three-dimensional motion trajectory curve of the cell can more realistically reflect the cell's movement changes in three-dimensional space, increasing the realism of the cell's movement changes reflected by the cell's motion trajectory curve, and realizing a realistic reflection of the cell's movement changes in three-dimensional space.

[0045] Step 104: Generate the three-dimensional motion trajectory curve of the cell based on multiple three-dimensional coordinates of the cell.

[0046] In this embodiment of the invention, multiple diffraction ring images of microspheres engulfed by cells can be acquired, and multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images can be automatically determined. Since the cells engulf the microspheres, the microspheres are subject to the dual constraints of the cell's cytoskeleton and cell membrane. The microspheres can effectively track and reproduce the movement changes of the cells. Therefore, the multiple three-dimensional coordinates of the microspheres can be determined as the multiple three-dimensional coordinates of the cells. Finally, a three-dimensional motion trajectory curve of the cells is generated based on the multiple three-dimensional coordinates of the cells. The three-dimensional motion trajectory curve of the cells can more realistically reflect the movement changes of the cells in three-dimensional space, increasing the realism of the cell movement changes reflected by the three-dimensional motion trajectory curve of the cells, and achieving the goal of realistically reflecting the movement changes of cells in three-dimensional space.

[0047] In some embodiments, the outer surface of the microspheres is coated with a layer of mucin. During the process of co-culturing with cells in a static environment, the microspheres are autonomously phagocytosed into the cells, and the cells can maintain a certain level of activity. In this way, the microspheres are subject to the dual constraints of the cytoskeleton and cell membrane of the active cells, allowing for good tracking and reproduction of cell movement changes. Therefore, multiple three-dimensional coordinates of the microspheres can be used as multiple three-dimensional coordinates of the cells. Based on the multiple three-dimensional coordinates of the cells, a three-dimensional motion trajectory curve of the cells can be generated. The three-dimensional motion trajectory curve of the cells can more realistically reflect the movement changes of the cells in three-dimensional space, increasing the realism of the cell movement changes reflected by the three-dimensional motion trajectory curve, and realizing the realistic reflection of cell movement changes in three-dimensional space.

[0048] In this embodiment of the invention, the microspheres can be silica spheres; the cells can be cardiomyocytes to which the drug is added. The drug can be various test drugs such as isoproterenol or metoprolol. After the cardiomyocytes engulf the microspheres, the drug can be added to the cardiomyocytes, and then multiple diffraction ring images of the engulfed microspheres can be acquired. Multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images are automatically determined. Since the cardiomyocytes engulf the microspheres, the microspheres are subject to the dual constraints of the cardiomyocyte cytoskeleton and the cardiomyocyte membrane. The microspheres can effectively track and reproduce the movement changes of the cardiomyocytes after the drug is added. The movement trajectory curve of the cardiomyocytes can more realistically reflect the movement changes of the cardiomyocytes in three-dimensional space after the drug is added, increasing the realism of the movement changes of the cardiomyocytes reflected by the three-dimensional movement trajectory curve after the drug is added. This achieves the goal of realistically reflecting the movement changes of the cardiomyocytes after the drug is added in three-dimensional space, thereby realizing the detection of the drug effect of the test drug and providing a more accurate and convenient method for drug effect detection. When either isoproterenol or metoprolol is added, the system can realistically reflect the changes in myocardial cell movement after the addition of isoproterenol or metoprolol in three-dimensional space, thereby enabling the detection of the drug effect of isoproterenol or metoprolol.

[0049] The method for generating cell motion trajectory curves provided in the embodiments of the present invention is further described below, such as... Figure 3 As shown, Figure 3 This is another flowchart illustrating the method for generating cell motion trajectory curves provided in this embodiment of the invention, which may specifically include the following steps:

[0050] Step 201: Obtain multiple diffraction ring images of the microspheres engulfed by the cells.

[0051] Step 202: Determine the diffraction ring image corresponding to the smallest image number among multiple diffraction ring images as the current diffraction ring image.

[0052] Step 203: Obtain the initial two-dimensional center coordinates of the current diffraction ring image.

[0053] Step 204: Dilate the pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates to obtain multiple pixel points.

[0054] Step 205: Segment the current diffraction ring image using multiple pixels as the segmentation origin to obtain multiple sets of segmented images.

[0055] In one optional implementation, the current diffraction ring image can be divided into a preset number of sub-images along 360°, using the pixels as the segmentation origin. The preset number is a positive integer, such as 8. A larger preset number results in more sub-images and a higher accuracy in determining the target two-dimensional center coordinates of the current diffraction ring image based on multiple segmented images. The number of segmented image groups is the same as the number of pixels, and each segmented image group includes the same number of sub-images as the preset number.

[0056] Step 206: Determine the target two-dimensional center coordinates of the current diffraction ring image based on multiple sets of segmented images.

[0057] Each group of segmented images includes multiple sub-images.

[0058] If the segmentation origin is the center point of the diffraction ring, then the intensity curves of multiple sub-images in each group of segmented images coincide. Ideally, the multiple ordinates corresponding to the same abscissa of the intensity curves of multiple sub-images are the same, and the variance of the ordinates corresponding to the same abscissa of the intensity curves of multiple sub-images is zero. Therefore, when the variance of the ordinates corresponding to the same abscissa of the intensity curves of multiple sub-images is zero, the segmentation origin of the segmented image group containing the multiple sub-images is the center point of the current diffraction ring image. However, since noise in the image is unavoidable, the segmentation origin of the segmented image group corresponding to the smallest ordinate variance can be determined as the center point of the current diffraction ring image; the two-dimensional coordinates of the center point of the current diffraction ring image can be determined as the target two-dimensional center coordinates of the current diffraction ring image. Therefore, in an optional implementation, the intensity curves of multiple sub-images in multiple groups are obtained, and the intensity curves include the ordinates of multiple points, with the abscissas of each ordinate being the same; the target two-dimensional center coordinates of the current diffraction ring image are determined based on the ordinates of the multiple points.

[0059] Specifically, based on the ordinates of multiple points, the variance of the ordinates of the light intensity curves of multiple sets of segmented images is determined; based on the variance of the ordinates of the light intensity curves of multiple sets of segmented images, the target two-dimensional center coordinates of the current diffraction ring image are determined.

[0060] Furthermore, the mean of the ordinates of multiple segmented images can be determined based on the ordinates of multiple points; and the variance of the ordinates of the light intensity curves of multiple segmented images can be determined based on the mean of the ordinates and the ordinates of multiple points.

[0061] In this embodiment of the invention, the mean of the vertical coordinates of each group of segmented images and the multiple vertical coordinates corresponding to the same horizontal coordinate in the light intensity curves of multiple sub-images of each group of segmented images can be substituted into the vertical coordinate variance calculation formula to calculate the vertical coordinate variance of the light intensity curves of each group of segmented images.

[0062] For example, the current diffraction ring image is image P1. The pixels of the current diffraction ring image corresponding to the initial two-dimensional center coordinates of the diffraction ring are obtained as follows: Figure 4 The image shows 25 pixels. The current diffraction ring image is segmented using these 25 pixels as the origin, resulting in 25 groups of segmented images. The first group of segmented images includes three sub-images: p1, p2, and p3. The intensity curve of p1 is S1, the intensity curve of p2 is S2, and the intensity curve of p3 is S3. Intensity curve S1 includes point A with coordinates (x1, ay1), intensity curve S2 includes point B with coordinates (x1, by1), and intensity curve S3 includes point C with coordinates (x1, cy1). Here, x1 is the x-coordinate of points A, B, and C, ay1 is the y-coordinate of point A, by1 is the y-coordinate of point B, and cy1 is the y-coordinate of point C. The mean ordinate of the first group of segmented images is determined as Y = (ay1 + by1 + cy1) / 3. Substituting the mean ordinate of the first group of segmented images and the multiple ordinates corresponding to the same x-axis in the intensity curves of multiple sub-images of the first group of segmented images into the formula for calculating the ordinate variance, the ordinate variance of the intensity curve of the first group of segmented images is calculated, resulting in D1. The process for determining the ordinate variance of the intensity curves of the remaining 24 groups of segmented images is the same and will not be repeated. Assuming that the smallest ordinate variance among the 25 groups of segmented images is D3, the origin of the segmentation of the group of segmented images corresponding to the smallest ordinate variance D3 can be determined as the center point of the current diffraction ring image. The two-dimensional coordinates (X1, Y1) of the center point of the current diffraction ring image are then determined as the target two-dimensional center coordinates (X1, Y1) of the current diffraction ring image (P1).

[0063] The formula for calculating the variance of the ordinate is as follows:

[0064]

[0065] D represents the variance of the ordinate, Y represents the mean of the ordinate, and y i This represents the ordinate of each point on the light intensity curve, where n represents the number of ordinates.

[0066] Step 207: Determine the three-dimensional coordinates of the microsphere based on the two-dimensional center coordinates of the target and the current diffraction ring image.

[0067] In this embodiment of the invention, since the microsphere, located at different positions, will produce different diffraction ring images due to interference effects, and the center of the diffraction ring in each diffraction ring image is the center of the microsphere, the target two-dimensional center coordinates of each diffraction ring image can be determined as the abscissa and ordinate of the microsphere's three-dimensional coordinates. Furthermore, since the diffraction ring image is segmented with the center of the diffraction ring as the origin, the intensity curve of the diffraction ring image can be obtained, and the lowest point of the intensity curve is the ordinate of the diffraction ring image. Therefore, in an optional implementation, the current diffraction ring image is segmented with the target two-dimensional center coordinates as the origin to obtain the target intensity curve of the current diffraction ring image; the ordinates of multiple points on the target intensity curve are obtained; the minimum ordinate is determined as the ordinate of the current diffraction ring image; and the three-dimensional coordinates composed of the target two-dimensional center coordinates and the ordinate are determined as the three-dimensional coordinates of the microsphere.

[0068] For example, multiple diffraction ring images include image P1, image P2, and image P3. The current diffraction ring image is image P1, and the target two-dimensional center coordinates of image P1 are (X1, Y1). The current diffraction ring image is segmented with the target two-dimensional center coordinates as the segmentation origin, resulting in the following: Figure 5 The target intensity curve L1 of the current diffraction ring image is shown. The ordinates of multiple points on the target intensity curve L1 are obtained. The ordinate Z1 of the lowest point in the target intensity curve L1 is determined as the vertical coordinate of the current diffraction ring image. The coordinates composed of the target's two-dimensional center coordinates and the vertical coordinates of the current diffraction ring image are determined as the three-dimensional coordinates of the microsphere, i.e., the target three-dimensional center coordinates of image P1 are (X1, Y1, Z1). The process of determining the target three-dimensional center coordinates of images P2 and P3 is the same as that of P1 and will not be repeated. Finally, the target three-dimensional center coordinates (X1, Y1, Z1) of image P1, (X2, Y2, Z2) of image P2, and (X3, Y3, Z3) of image P3 are obtained.

[0069] Step 208: Determine whether the image number of the current diffraction ring image is equal to the preset threshold. If yes, proceed to step 211; otherwise, proceed to step 209.

[0070] Step 209: Determine the target two-dimensional center coordinates of the current diffraction ring image as the initial two-dimensional center coordinates of the diffraction ring image corresponding to the next image number of the current diffraction ring image.

[0071] Step 210: Determine the diffraction ring image corresponding to the next image number after the current diffraction ring image as the current diffraction ring image.

[0072] After executing step 210, return to execute step 203.

[0073] Step 211: Determine the multiple three-dimensional coordinates of the microspheres as multiple three-dimensional coordinates of the cell.

[0074] Step 212: Generate the three-dimensional motion trajectory curve of the cell based on multiple three-dimensional coordinates of the cell.

[0075] In this embodiment of the invention, multiple diffraction ring images of microspheres engulfed by cells are acquired; multiple three-dimensional coordinates of the microspheres are determined based on the multiple diffraction ring images; the multiple three-dimensional coordinates of the microspheres are then used as multiple three-dimensional coordinates of the cells; and a three-dimensional motion trajectory curve of the cells is generated based on the multiple three-dimensional coordinates of the cells. That is, the technical solution of this invention can acquire multiple diffraction ring images of microspheres engulfed by cells and automatically determine multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images. Because the cells engulf the microspheres, the microspheres are subject to the dual constraints of the cell's cytoskeleton and cell membrane, allowing the microspheres to effectively track and reproduce the cell's motion changes. Therefore, the multiple three-dimensional coordinates of the microspheres can be used as multiple three-dimensional coordinates of the cells. Finally, a three-dimensional motion trajectory curve of the cells is generated based on the multiple three-dimensional coordinates of the cells. This three-dimensional motion trajectory curve more realistically reflects the cell's motion changes in three-dimensional space, increasing the realism of the cell's motion changes reflected by the three-dimensional motion trajectory curve, and achieving the goal of realistically reflecting the cell's motion changes in three-dimensional space.

[0076] Figure 6 This is a schematic diagram of a cell motion trajectory curve generation device provided in an embodiment of the present invention. This device is suitable for executing the cell motion trajectory curve generation method provided in an embodiment of the present invention. Figure 6 As shown, the device may specifically include:

[0077] Image acquisition module 301 is used to acquire multiple diffraction ring images of microspheres engulfed by cells;

[0078] The first determining module 302 is used to determine multiple three-dimensional coordinates of the microsphere based on the multiple diffraction ring images;

[0079] The second determining module 303 is used to determine the multiple three-dimensional coordinates of the microsphere as the multiple three-dimensional coordinates of the cell;

[0080] The curve generation module 304 is used to generate a three-dimensional motion trajectory curve of the cell based on multiple three-dimensional coordinates of the cell.

[0081] Optionally, the diffraction ring image includes an image number, and the first determining module 302 is specifically used for:

[0082] The diffraction ring image corresponding to the smallest image number among the multiple diffraction ring images is determined as the current diffraction ring image;

[0083] Obtain the initial two-dimensional center coordinates of the current diffraction ring image;

[0084] Dilatation of the pixels in the current diffraction ring image corresponding to the initial two-dimensional center coordinates yields multiple pixel points;

[0085] The current diffraction ring image is segmented using the plurality of pixels as the segmentation origin to obtain multiple sets of segmented images;

[0086] Based on the multiple sets of segmented images, determine the target two-dimensional center coordinates of the current diffraction ring image;

[0087] The three-dimensional coordinates of the microsphere are determined based on the two-dimensional center coordinates of the target and the current diffraction ring image.

[0088] When the image number of the current diffraction ring image is equal to a preset threshold, the step of determining the multiple three-dimensional coordinates of the microsphere as the multiple three-dimensional coordinates of the cell is triggered.

[0089] When the image number of the current diffraction ring image is less than the preset threshold, the target two-dimensional center coordinates of the current diffraction ring image are determined as the initial two-dimensional center coordinates of the diffraction ring image corresponding to the next image number of the current diffraction ring image.

[0090] The diffraction ring image corresponding to the next image number of the current diffraction ring image is determined as the current diffraction ring image, triggering the step of expanding the diffraction ring corresponding to the initial two-dimensional center coordinates with the initial two-dimensional center coordinates as the expansion origin to obtain an expanded image.

[0091] Optionally, each set of segmented images includes multiple sub-images. The first determining module 302 determines the target two-dimensional center coordinates of the current diffraction ring image based on the multiple sets of segmented images, including:

[0092] Obtain the light intensity curves of the multiple sub-images in the multiple groups, wherein the light intensity curves include the ordinates of multiple points, and the abscissas of each ordinate are the same;

[0093] Based on the ordinates of the multiple points, the target two-dimensional center coordinates of the current diffraction ring image are determined.

[0094] Optionally, the first determining module 302 determines the target two-dimensional center coordinates of the current diffraction ring image based on the ordinates of the plurality of points, including:

[0095] Based on the ordinates of the multiple points, determine the variance of the ordinates of the light intensity curves of the multiple segmented images;

[0096] The target two-dimensional center coordinates of the current diffraction ring image are determined based on the variance of the ordinate of the light intensity curves of the multiple sets of segmented images.

[0097] Optionally, the first determining module 302 determines the variance of the ordinate of the light intensity curves of multiple sets of segmented images based on the ordinates of the multiple points, including:

[0098] Based on the ordinates of the multiple points, determine the mean ordinate of the multiple sets of segmented images;

[0099] Based on the mean of the ordinate and the ordinates of the multiple points, the variance of the ordinate of the light intensity curves of the multiple sets of segmented images is determined.

[0100] Optionally, the first determining module 302 determines the target two-dimensional center coordinates of the current diffraction ring image based on the variance of the ordinate of the light intensity curves of the multiple sets of segmented images, including:

[0101] The ordinate corresponding to the minimum ordinate variance is determined as the two-dimensional center ordinate of the current diffraction ring image;

[0102] The x-coordinate corresponding to the y-coordinate is determined as the two-dimensional center x-coordinate of the current diffraction ring image;

[0103] The two-dimensional coordinates composed of the two-dimensional center ordinate and the two-dimensional center abscissa are determined as the target two-dimensional center coordinates of the current diffraction ring image.

[0104] Optionally, the first determining module 302 determines the three-dimensional coordinates of the microsphere based on the target's two-dimensional center coordinates and the current diffraction ring image, including:

[0105] The current diffraction ring image is segmented using the two-dimensional center coordinates of the target as the segmentation origin to obtain the target light intensity curve of the current diffraction ring image;

[0106] Obtain the ordinates of multiple points on the target light intensity curve;

[0107] The minimum ordinate is determined as the vertical coordinate of the current diffraction ring image;

[0108] The three-dimensional coordinates, composed of the two-dimensional center coordinates of the target and the vertical coordinates, are determined as the three-dimensional coordinates of the microsphere.

[0109] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional modules is merely an example. In practical applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. The specific working process of the functional modules described above can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0110] The cell motion trajectory curve generation device provided in this invention can acquire multiple diffraction ring images of microspheres engulfed by cells and automatically determine multiple three-dimensional coordinates of the microspheres in the multiple diffraction ring images. Since the cells engulf the microspheres, the microspheres are subject to the dual constraints of the cell's cytoskeleton and cell membrane. The microspheres can effectively track and reproduce the cell's motion changes. Therefore, the multiple three-dimensional coordinates of the microspheres can be determined as the multiple three-dimensional coordinates of the cell. Finally, a three-dimensional motion trajectory curve of the cell is generated based on the multiple three-dimensional coordinates of the cell. The three-dimensional motion trajectory curve of the cell can more realistically reflect the cell's motion changes in three-dimensional space, increasing the realism of the cell's motion changes reflected by the three-dimensional motion trajectory curve, and achieving the goal of realistically reflecting the cell's motion changes in three-dimensional space.

[0111] Figure 7 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention.

[0112] Please refer to Figure 7 An electronic device 50 is provided, comprising:

[0113] Processor 51; and,

[0114] Memory 52 is used to store the executable instructions of the processor;

[0115] The processor 51 is configured to execute the methods described above by executing the executable instructions.

[0116] The processor 51 can communicate with the memory 52 via the bus 53.

[0117] This invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the methods described above.

[0118] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0119] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for generating a cell motion trajectory curve, characterized in that, The method includes: Acquire multiple diffraction ring images of microspheres engulfed by cells; Based on the multiple diffraction ring images, the multiple three-dimensional coordinates of the microsphere are determined; The multiple three-dimensional coordinates of the microspheres are determined as the multiple three-dimensional coordinates of the cell; The three-dimensional motion trajectory curve of the cell is generated based on multiple three-dimensional coordinates of the cell.

2. The method according to claim 1, characterized in that, The diffraction ring images include image numbers, and determining multiple three-dimensional coordinates of the microsphere based on the multiple diffraction ring images includes: The diffraction ring image corresponding to the smallest image number among the multiple diffraction ring images is determined as the current diffraction ring image; Obtain the initial two-dimensional center coordinates of the current diffraction ring image; Dilatation of the pixels in the current diffraction ring image corresponding to the initial two-dimensional center coordinates yields multiple pixel points; The current diffraction ring image is segmented using the plurality of pixels as the segmentation origin to obtain multiple sets of segmented images; Based on the multiple sets of segmented images, determine the target two-dimensional center coordinates of the current diffraction ring image; The three-dimensional coordinates of the microsphere are determined based on the two-dimensional center coordinates of the target and the current diffraction ring image. When the image number of the current diffraction ring image is equal to a preset threshold, the step of determining the multiple three-dimensional coordinates of the microsphere as the multiple three-dimensional coordinates of the cell is triggered. When the image number of the current diffraction ring image is less than the preset threshold, the target two-dimensional center coordinates of the current diffraction ring image are determined as the initial two-dimensional center coordinates of the diffraction ring image corresponding to the next image number of the current diffraction ring image. The diffraction ring image corresponding to the next image number of the current diffraction ring image is determined as the current diffraction ring image, triggering the step of expanding the diffraction ring corresponding to the initial two-dimensional center coordinates with the initial two-dimensional center coordinates as the expansion origin to obtain an expanded image.

3. The method according to claim 2, characterized in that, Each group of segmented images includes multiple sub-images. Determining the target two-dimensional center coordinates of the current diffraction ring image based on the multiple groups of segmented images includes: Obtain the light intensity curves of the multiple sub-images in the multiple groups, wherein the light intensity curves include the ordinates of multiple points, and the abscissas of each ordinate are the same; Based on the ordinates of the multiple points, the target two-dimensional center coordinates of the current diffraction ring image are determined.

4. The method according to claim 3, characterized in that, Determining the target two-dimensional center coordinates of the current diffraction ring image based on the ordinates of the plurality of points includes: Based on the ordinates of the multiple points, determine the variance of the ordinates of the light intensity curves of the multiple segmented images; The target two-dimensional center coordinates of the current diffraction ring image are determined based on the variance of the ordinate of the light intensity curves of the multiple sets of segmented images.

5. The method according to claim 4, characterized in that, The step of determining the variance of the ordinate of the light intensity curves of multiple sets of segmented images based on the ordinates of the multiple points includes: Based on the ordinates of the multiple points, determine the mean ordinate of the multiple sets of segmented images; Based on the mean of the ordinate and the ordinates of the multiple points, the variance of the ordinate of the light intensity curves of the multiple sets of segmented images is determined.

6. The method according to claim 4, characterized in that, Determining the target two-dimensional center coordinates of the current diffraction ring image based on the variance of the ordinate of the light intensity curves of the multiple sets of segmented images includes: The ordinate corresponding to the minimum ordinate variance is determined as the two-dimensional center ordinate of the current diffraction ring image; The x-coordinate corresponding to the y-coordinate is determined as the two-dimensional center x-coordinate of the current diffraction ring image; The two-dimensional coordinates composed of the two-dimensional center ordinate and the two-dimensional center abscissa are determined as the target two-dimensional center coordinates of the current diffraction ring image.

7. The method according to claim 2, characterized in that, Determining the three-dimensional coordinates of the microsphere based on the target's two-dimensional center coordinates and the current diffraction ring image includes: The current diffraction ring image is segmented using the two-dimensional center coordinates of the target as the segmentation origin to obtain the target light intensity curve of the current diffraction ring image; Obtain the ordinates of multiple points on the target light intensity curve; The minimum ordinate is determined as the vertical coordinate of the current diffraction ring image; The three-dimensional coordinates, composed of the two-dimensional center coordinates of the target and the vertical coordinates, are determined as the three-dimensional coordinates of the microsphere.

8. A device for generating the jumping trajectory of a cell, characterized in that, The device includes: The image acquisition module is used to acquire multiple diffraction ring images of microspheres engulfed by cells; The first determining module is used to determine multiple three-dimensional coordinates of the microsphere based on the multiple diffraction ring images; The second determining module is used to determine the multiple three-dimensional coordinates of the microsphere as multiple three-dimensional coordinates of the cell; The curve generation module is used to generate a three-dimensional motion trajectory curve of the cell based on multiple three-dimensional coordinates of the cell.

9. An electronic device, characterized in that, Including processor and memory, The memory is used to store code and related data; The processor is configured to execute code in the memory to implement the cell motion trajectory curve generation method according to any one of claims 1 to 7.

10. A storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for generating the motion trajectory curve of a cell according to any one of claims 1 to 7.

Images