A method and apparatus for transient analysis of myocardial microspheres

Abstract

The application provides a myocardial microsphere transient analysis method and device. The method comprises the following steps: obtaining a myocardial microsphere video, using a segmentation method to obtain a myocardial microsphere region of each frame in the myocardial microsphere video, determining the position of all myocardial microspheres of each frame according to the myocardial microsphere region obtained by each frame, finding corresponding myocardial microspheres between frames according to the nearest distance of myocardial microspheres between frames, selecting a reference frame and calculating a jumping strength and a jumping speed by using a frame difference method and a maximum value projection method, and generating a jumping strength curve and a jumping speed curve respectively, determining each wave peak according to the jumping curve, calculating a frequency according to the wave peak, and performing transient analysis on each wave peak. The myocardial microsphere region in the video is extracted by using the segmentation method to be analyzed separately, and the effect of simultaneously analyzing multiple myocardial microsphere videos can be achieved. The maximum value projection of the frame difference method result is obtained first, then the reference frame is automatically selected, and then the frame difference method result is multiplied by the maximum value projection, so that the noise can be suppressed.

Description

Technical Field

[0001] This invention belongs to the field of image processing of in vitro three-dimensional cell models, specifically relating to a transient analysis method and device for myocardial microspheres. Background Technology

[0002] Cardiotoxicity is a major cause of new drug development failures and the withdrawal of marketed drugs. Traditional preclinical studies primarily rely on cell screening, animal models, and non-cardiac cells expressing heterogeneous myocardial ion channels to evaluate drug safety. However, because the research subjects are all non-human cardiomyocytes, the pharmacodynamics and toxicology of candidate drugs cannot be precisely elucidated. Furthermore, animal models have limitations such as high cost and low throughput, and due to species differences, many drugs that pass animal experiments fail to pass human trials. A survey by the U.S. Food and Drug Administration (FDA) shows that 92% of drugs that have proven safe and effective in animal experiments fail in clinical human trials and cannot be truly put into production and market. This not only increases the financial consumption and experimental cycle of drug development but also inevitably causes huge waste, creating an awkward situation of high investment, high risk, and low output in the field of new drug development.

[0003] There is widespread interest in the development of in vitro models that can mimic tissues, organs, and their functional behavior, for example, to allow for non-invasive, rapid, economical, and regenerative testing and / or screening of new drug, chemical, and food candidates.

[0004] Directed differentiation of human pluripotent stem cells (hiPSCs) is crucial for producing cardiomyocytes expressing conventional cardiomyocyte-specific genes, as well as various contractile proteins and ion channels. Cardiomyocyte microspheres obtained by 3D culture of 2D cardiomyocytes differentiated from hiPSCs possess more mature structure, electrophysiological activity, and the ability to periodically contract and relax compared to ordinary 2D cardiomyocytes. Measuring the contraction and electrophysiological parameters of these cardiomyocyte microspheres is therefore essential.

[0005] Therefore, there is an urgent need to quantify cardiac contractile capacity in in vitro models at multiple levels, but existing research methods lack a quantitative evaluation method for myocardial contractile parameters in in vitro models. Summary of the Invention

[0006] The present invention aims to solve at least one of the technical problems existing in the prior art, and to provide a method and apparatus for transient analysis of myocardial microspheres.

[0007] One aspect of the present invention provides a transient analysis method for myocardial microspheres, the method comprising the following steps:

[0008] Acquire a video of myocardial microspheres and use a segmentation method to obtain the myocardial microsphere region in each frame of the video;

[0009] Based on the myocardial microsphere region obtained in each frame, the position of all myocardial microspheres in each frame is determined, and the corresponding myocardial microspheres between frames are found based on the nearest distance between myocardial microspheres between frames.

[0010] Using the frame difference method and the maximum value projection method, a reference frame is selected and the jitter intensity and jitter velocity are calculated to generate jitter intensity curves and jitter velocity curves respectively;

[0011] Each peak is determined based on the fluctuation curve, the frequency is calculated based on the peak, and transient analysis is performed on each peak.

[0012] In some preferred embodiments, the step of selecting a reference frame and calculating jitter intensity and jitter velocity using the frame difference method and the maximum value projection method, and generating jitter intensity curves and jitter velocity curves respectively, includes:

[0013] Using the frame difference method, each frame is subtracted from the frame two frames earlier, and the maximum value at each position is obtained by projecting the maximum value onto all the results.

[0014] The reference frame is automatically selected or manually selected according to the formula. Each frame is subtracted from the reference frame and multiplied pixel by pixel with the result of the maximum value projection method to obtain the average pixel value as the jitter intensity of the frame, and the jitter intensity curve is generated.

[0015] Each frame is subtracted from the frame two frames prior, and the average pixel value is calculated as the frame's jump speed, generating the jump speed curve.

[0016] In some alternative implementations, the automatic selection of the reference frame is formulated as follows:

[0017] m i =mean(f i -f i+2 )

[0018] f i For the i-th frame of the image, mean represents the average pixel value;

[0019]

[0020] Iterate through r from smallest to largest as r j Then the corresponding value is m j :

[0021]

[0022] Traverse u from smallest to largest as u i Then the corresponding value is m j :

[0023]

[0024] Then the minimum value of v is vk The corresponding j is the automatically selected reference frame coordinate.

[0025] In some optional implementations, the step of determining each peak based on the fluctuation curve, calculating the frequency based on the peaks, and performing transient analysis on each peak includes:

[0026] On the vibration intensity curve, if a value is the maximum of the five values ​​before and after it and is greater than the threshold, it is the peak, and the frequency is calculated based on the peak.

[0027] Find the time points 10%, 20%, ... 90% before and after each peak, and calculate the time from the peak. The time before the peak is the contraction time, and the time after the peak is the relaxation time.

[0028] In some preferred embodiments, determining the positions of all myocardial microspheres in each frame based on the myocardial microsphere region obtained in each frame, and finding corresponding myocardial microspheres between frames based on the nearest distance between myocardial microspheres between frames, includes:

[0029] The centroid coordinates are calculated based on the first and zero moments of the myocardial microsphere region contour to obtain the position of all myocardial microspheres in each frame.

[0030] Using the coordinates of the myocardial microspheres in the first frame as a reference, for each myocardial microsphere, the myocardial microsphere closest to it is found from the myocardial microspheres in other frames.

[0031] In some preferred embodiments, the moments and centroid coordinates of the myocardial microsphere region contour are calculated as follows:

[0032]

[0033]

[0034]

[0035]

[0036]

[0037] In the binary image, array(x,y) at the contour is always 1, therefore m 00 It is the number of contour points, m 10 It is the sum of the x-coordinates of all contour points, m 01 It is the sum of the y-coordinates of all contour points, and (X,Y) are the centroid coordinates.

[0038] In some preferred embodiments, when the bounce intensity curve shows an upward sloping trend, the method further includes, before determining each peak based on the bounce curve:

[0039] The troughs are determined based on the vibration intensity curve, and the troughs are used as reference frames to recalculate the vibration intensity and generate a new vibration intensity curve.

[0040] In another aspect, the present invention provides a transient analysis device for myocardial microspheres, the device comprising:

[0041] The segmentation module is used to acquire myocardial microsphere videos and use segmentation methods to obtain the myocardial microsphere region in each frame of the myocardial microsphere video;

[0042] The determination module is used to determine the position of all myocardial microspheres in each frame based on the myocardial microsphere region obtained in each frame, and to find the corresponding myocardial microspheres between frames based on the closest distance between myocardial microspheres between frames.

[0043] The calculation module is used to select a reference frame and calculate the jitter intensity and jitter velocity using the frame difference method and the maximum value projection method, and generate the jitter intensity curve and jitter velocity curve respectively.

[0044] The analysis module is used to determine each peak based on the fluctuation curve, calculate the frequency based on the peak, and perform transient analysis on each peak.

[0045] Another aspect of the present invention provides an electronic device comprising:

[0046] One or more processors;

[0047] A storage unit for storing one or more programs that, when executed by one or more processors, enable the one or more processors to implement the method described above.

[0048] In another aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, enables the implementation of the method described above.

[0049] In another aspect, the present invention provides an in vitro method for using the above-described transient analysis method of myocardial microspheres to detect the regulatory effect of drugs on the function of myocardial microspheres in vitro.

[0050] The transient analysis method and apparatus for myocardial microspheres in this invention utilizes a segmentation method to extract the myocardial microsphere region from a video for separate analysis, achieving the effect of simultaneously analyzing videos containing multiple myocardial microspheres. Using only the frame difference method will generate noise. This invention first obtains the maximum value projection of the frame difference method result, then automatically selects a reference frame, and finally multiplies the frame difference method result by the maximum value projection to suppress noise. Attached Figure Description

[0051] Figure 1 This is a schematic block diagram of the composition of an electronic device according to an embodiment of the present invention;

[0052] Figure 2 This is a flowchart of a transient analysis method for myocardial microspheres according to another embodiment of the present invention;

[0053] Figure 3 This is a flowchart of a transient analysis method for myocardial microspheres according to another embodiment of the present invention;

[0054] Figure 4 This is a schematic diagram of the first frame segmentation result according to another embodiment of the present invention;

[0055] Figure 5 This is a schematic diagram of the jitter curve after automatically selecting a reference frame according to another embodiment of the present invention;

[0056] Figure 6 This is a schematic diagram of transient analysis results according to another embodiment of the present invention;

[0057] Figure 7 This is a schematic diagram of the jitter curve after using the trough as a reference frame in another embodiment of the present invention;

[0058] Figure 8 This is a schematic diagram of the transient analysis device for myocardial microspheres according to another embodiment of the present invention. Detailed Implementation

[0059] To enable those skilled in the art to better understand the technical solution of the present invention, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0060] First, refer to Figure 1 This describes an example electronic device for implementing a transient analysis method and apparatus for myocardial microspheres according to embodiments of the present invention.

[0061] like Figure 1 As shown, the electronic device 200 includes one or more processors 210, one or more storage devices 220, one or more input devices 230, one or more output devices 240, etc., and these components are interconnected via a bus system 250 and / or other forms of connection mechanisms. It should be noted that... Figure 1 The components and structures of the electronic devices shown are merely exemplary and not limiting; other components and structures may be used as needed.

[0062] The processor 210 may be a central processing unit (CPU), or may be a processing unit consisting of multiple processing cores, or other forms of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device 200 to perform desired functions.

[0063] Storage device 220 may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, which a processor may execute to implement the client functions (implemented by the processor) in the embodiments of the present invention described below, and / or other desired functions. Various applications and various data may also be stored in the computer-readable storage medium, such as various data used and / or generated by the applications.

[0064] The input device 230 may be a device used by a user to input commands, and may include one or more of a keyboard, mouse, microphone, and touch screen.

[0065] The output device 240 can output various information (such as images or sounds) to the outside (e.g., a user) and may include one or more of a display, a speaker, etc.

[0066] Below, we will refer to Figure 2 and Figure 3 A transient analysis method for myocardial microspheres according to another embodiment of the present invention is described.

[0067] like Figure 2 and Figure 3 As shown, a transient analysis method S100 for myocardial microspheres includes the following steps:

[0068] S110. Acquire a video of myocardial microspheres and use a segmentation method to obtain the myocardial microsphere region in each frame of the video.

[0069] Specifically, in this step, the Otsu thresholding method is used to separate the myocardial microsphere region from the background in each frame of the myocardial microsphere video. The Otsu thresholding method is a maximum inter-class variance method. After image binarization segmentation using the threshold obtained by this method, the inter-class variance between the foreground and background images is maximized. The segmentation result is as follows: Figure 4 As shown, Figure 4 The left image shows a single myocardial microsphere, while the right image shows multiple myocardial microspheres.

[0070] The Otsu threshold segmentation process is as follows:

[0071] For an image I(x,y), the segmentation threshold between the foreground (i.e., the target) and the background is denoted as T. The proportion of pixels belonging to the foreground in the entire image is denoted as ω0, and its average gray level is μ0. The proportion of pixels belonging to the background in the entire image is denoted as ω1, and its average gray level is μ1. The overall average gray level of the image is denoted as μ, and the inter-class variance is denoted as g.

[0072] Assuming the background of the image is dark and the image size is M×N, let N0 be the number of pixels with a gray value less than a threshold T, and N1 be the number of pixels with a gray value greater than the threshold T. Then:

[0073] ω0=N0 / (M×N) (1)

[0074] ω1=N1 / (M×N) (2)

[0075] N0+N1=M×N (3)

[0076] ω0+ω1=1 (4)

[0077] μ=ω0*μ0+ω1*μ1 (5)

[0078] g = ω0 * (μ0 - μ) 2 +ω1*(μ1-μ) 2 (6)

[0079] Substituting equation (5) into equation (6), we obtain the equivalent formula:

[0080] g = ω0 * ω1 * (μ0 - μ1) 2 (7)

[0081] The threshold T that maximizes the inter-class variance g is obtained by traversing the data.

[0082] S120. Based on the myocardial microsphere region obtained in each frame, determine the position of all myocardial microspheres in each frame, and find the corresponding myocardial microspheres between frames based on the closest distance between myocardial microspheres between frames.

[0083] Specifically, in this step, the centroid coordinates are calculated based on the first and zero moments of the myocardial microsphere region contour to obtain the positions of all myocardial microspheres in each frame. Using the myocardial microsphere coordinates of the first frame as a reference, for each myocardial microsphere, the closest myocardial microsphere is found from the myocardial microspheres in other frames, thus laying the foundation for individual analysis of each myocardial microsphere. The calculation methods for the contour moments and centroid coordinates are as follows.

[0084] m 00 =∑ x,y array(x,y) (8)

[0085] m 10 =∑x,y array(x,y)*x (9)

[0086] m 01 =∑ x,y array(x,y)*y (10)

[0087]

[0088]

[0089] In the binary image, array(x,y) at the contour is always 1, therefore m 00 It is the number of contour points, m 10 It is the sum of the x-coordinates of all contour points, m 01 It is the sum of the y-coordinates of all contour points, and (X,Y) are the centroid coordinates.

[0090] S130. Using the frame difference method and the maximum value projection method, select a reference frame and calculate the jitter intensity and jitter velocity, and generate the jitter intensity curve and jitter velocity curve respectively.

[0091] Specifically, in this step, the frame difference method is used to subtract each frame from the frame two frames prior, and the maximum value at each position is obtained by projecting the maximum value onto all results. A reference frame is automatically or manually selected according to a formula. Each frame is subtracted from the reference frame and multiplied pixel-by-pixel by the result of the maximum value projection method to calculate the average pixel value as the jitter intensity of that frame, and the jitter intensity curve is generated. Similarly, each frame is subtracted from the frame two frames prior, and the average pixel value is calculated as the jitter velocity of that frame, and the jitter velocity curve is generated. Figure 5 As shown, in Figure 5 In the middle, the left figure shows the beating intensity curve of the myocardial microspheres, the right figure shows the beating velocity curve of the myocardial microspheres, the top figure shows the curve of a single myocardial microsphere, and the bottom figure shows the curves of multiple myocardial microspheres.

[0092] The formula for automatically selecting the reference frame is as follows:

[0093] m i =mean(f i -f i+2 (13)

[0094] f i Let be the i-th frame of the image, and mean represent the average pixel value.

[0095]

[0096] Iterate through r from smallest to largest as r j Then the corresponding value is m j :

[0097]

[0098] Traverse u from smallest to largest as u i Then the corresponding value is m j :

[0099]

[0100] Then the minimum value of v is v k The corresponding j is the automatically selected reference frame coordinate.

[0101] S140. Determine each peak based on the fluctuation curve, calculate the frequency based on the peak, and perform transient analysis on each peak.

[0102] Specifically, in this step, on the bounce intensity curve, if a value is the maximum of the five values ​​before and after it and is greater than a threshold, it is considered a peak. The frequency is calculated based on the peak, and transient analysis is performed on each peak. Transient analysis refers to finding nine time points before and after each peak, at 10%, 20%,...90% of the time interval, with the former representing the contraction time and the latter the diastolic time. The results are saved as text, such as... Figure 6 As shown.

[0103] In this embodiment, the method can generate a good beating curve for a video of one or more myocardial microspheres, and can effectively analyze the contraction and relaxation time of each beating, and obtain the peak value, which can be used to analyze the effectiveness of the drug on the myocardial microspheres.

[0104] The transient analysis method for myocardial microspheres in this invention utilizes a segmentation method to extract the myocardial microsphere region from the video for separate analysis, achieving the effect of simultaneously analyzing videos containing multiple myocardial microspheres. Using only the frame difference method will generate noise. This invention first obtains the maximum value projection of the frame difference method result, then automatically selects a reference frame, and finally multiplies the frame difference method result by the maximum value projection to suppress noise.

[0105] In some alternative embodiments, when the bounce intensity curve shows an upward sloping trend, the method further includes, before determining the individual peaks based on the bounce curve:

[0106] The troughs are determined based on the vibration intensity curve, and the troughs are used as reference frames to recalculate the vibration intensity and generate a new vibration intensity curve.

[0107] Specifically, in this step, if the curve shows an upward sloping trend, this step is performed; otherwise, it is skipped. On the jitter intensity curve, if a value is the minimum of the five values ​​before and after it and is less than a threshold, it is considered a trough. For each frame, the nearest trough before that frame is used as the reference frame, and the jitter intensity is recalculated to generate a jitter intensity curve. The jitter curve is as follows: Figure 7 As shown, Figure 7The left figure shows the myocardial microsphere beating intensity curve after adjusting for the trough as the reference frame, and the right figure shows the myocardial microsphere beating velocity curve. It can be seen that the intensity curve no longer tilts upward.

[0108] The transient analysis method for myocardial microspheres in this embodiment of the invention addresses the problem that the myocardial microspheres undergo slight displacement during beating, and using a single reference frame would cause the beating curve to tilt upwards. This embodiment of the invention solves this problem by reselecting a reference frame after each beating.

[0109] Another aspect of the present invention, such as Figure 8 As shown, a transient analysis device 100 for myocardial microspheres is provided. This device 100 is applicable to the analysis methods described above, and details can be found in the relevant previous descriptions, which will not be repeated here. The device 100 includes:

[0110] The segmentation module 110 is used to acquire myocardial microsphere video and use a segmentation method to obtain the myocardial microsphere region in each frame of the myocardial microsphere video;

[0111] The determining module 120 is used to determine the position of all myocardial microspheres in each frame based on the myocardial microsphere region obtained in each frame, and to find the corresponding myocardial microspheres between frames based on the closest distance between myocardial microspheres between frames.

[0112] The calculation module 130 is used to select a reference frame and calculate the jitter intensity and jitter velocity using the frame difference method and the maximum value projection method, and generate the jitter intensity curve and the jitter velocity curve respectively.

[0113] Analysis module 140 is used to determine each peak according to the jumping curve, calculate the frequency according to the peak, and perform transient analysis on each peak.

[0114] The transient analysis device for myocardial microspheres in this invention utilizes a segmentation method to extract the myocardial microsphere region from a video for separate analysis, achieving the effect of simultaneously analyzing videos containing multiple myocardial microspheres. Using only the frame difference method will generate noise. This invention first obtains the maximum value projection of the frame difference method result, then automatically selects a reference frame, and finally multiplies the frame difference method result by the maximum value projection to suppress noise.

[0115] In another aspect, the present invention provides an in vitro method for using the above-described transient analysis method of myocardial microspheres to detect the regulatory effect of drugs on the function of myocardial microspheres in vitro.

[0116] Specifically, in this embodiment, a drug is pre-injected into in vitro myocardial microspheres, and then the transient analysis method of the myocardial microspheres described above is used to obtain the beating curve of the in vitro myocardial microspheres with injected drug. Based on the beating curve, the contraction and relaxation time of each beating can be analyzed, and the peak value can be obtained, thereby analyzing the functional regulatory effect of the drug on the in vitro myocardial microspheres.

[0117] Another aspect of the present invention provides an electronic device comprising:

[0118] One or more processors;

[0119] A storage unit for storing one or more programs that, when executed by one or more processors, enable the one or more processors to implement the methods described above.

[0120] In another aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, enables the implementation of the method described above.

[0121] The computer-readable medium may be included in the apparatus, device, or system of the present invention, or it may exist independently.

[0122] The computer-readable storage medium may be any tangible medium that contains or stores a program, and may be an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device. More specific examples include, but are not limited to, electrical connections having one or more wires, portable computer disks, hard disks, optical fibers, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0123] The computer-readable storage medium may also include data signals propagated in baseband or as part of a carrier wave, carrying computer-readable program code, specific examples of which include, but are not limited to, electromagnetic signals, optical signals, or any suitable combination thereof.

[0124] It is understood that the above embodiments are merely exemplary implementations used to illustrate the principles of the present invention, and the present invention is not limited thereto. For those skilled in the art, various modifications and improvements can be made without departing from the spirit and essence of the present invention, and these modifications and improvements are also considered to be within the scope of protection of the present invention.

Claims

1. A transient analysis method for myocardial microspheres, characterized in that, The myocardial microspheres are obtained by 3D culture of 2D cardiomyocytes, and the method includes the following steps: Acquire video of myocardial microspheres. The myocardial microsphere region in each frame of the myocardial microsphere video is obtained using a segmentation method; the step of obtaining the myocardial microsphere region in each frame of the myocardial microsphere video using a segmentation method includes: separating the myocardial microsphere region from the background in each frame of the myocardial microsphere video using the Otsu threshold segmentation method; Based on the myocardial microsphere region obtained in each frame, the position of all myocardial microspheres in each frame is determined, and the corresponding myocardial microspheres between frames are found based on the nearest distance between myocardial microspheres between frames. Using the frame difference method and the maximum value projection method, a reference frame is selected and the jitter intensity and jitter velocity are calculated to generate jitter intensity curves and jitter velocity curves, respectively; the jitter curve includes the jitter intensity curve and the jitter velocity curve; Each peak is determined based on the jump curve, the frequency is calculated based on the peak, and transient analysis is performed on each peak. When the vibration intensity curve shows an upward sloping trend, before determining each peak based on the vibration curve, the method further includes: Based on the bounce intensity curve, the trough is determined, and the trough is used as a reference frame to recalculate the bounce intensity and generate a new bounce intensity curve, specifically as follows: On the jitter intensity curve, if a value is the minimum of the five values ​​before and after it and is less than the threshold, it is a trough. For each frame, the jitter intensity is recalculated and the jitter intensity curve is generated with the nearest trough before that frame as the reference frame.

2. The method according to claim 1, characterized in that, The process of selecting a reference frame and calculating jitter intensity and jitter velocity using the frame difference method and maximum projection method, and generating jitter intensity curves and jitter velocity curves respectively, includes: Using the frame difference method, each frame is subtracted from the frame two frames earlier, and the maximum value at each position is obtained by projecting the maximum value onto all the results. The reference frame is automatically selected or manually selected according to the formula. Each frame is subtracted from the reference frame and multiplied pixel by pixel with the result of the maximum value projection method to obtain the average pixel value as the jitter intensity of the frame, and the jitter intensity curve is generated. Each frame is subtracted from the frame two frames prior, and the average pixel value is calculated as the frame's jump speed, generating the jump speed curve.

3. The method according to claim 2, characterized in that, The formula for automatically selecting the reference frame is as follows: For the i-th frame of the image, mean represents the average pixel value; Then the minimum value of v The corresponding j is the automatically selected reference frame coordinate.

4. The method according to claim 1, characterized in that, The step of determining each peak based on the fluctuation curve, calculating the frequency based on the peaks, and performing transient analysis on each peak includes: On the vibration intensity curve, if a value is the maximum of the five values ​​before and after it and is greater than the threshold, it is the peak, and the frequency is calculated based on the peak. Find the time points 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90% before and after each peak value, and calculate the time from the peak value. The time before the peak value is the contraction time, and the time after the peak value is the relaxation time.

5. The method according to claim 1, characterized in that, The step of determining the position of all myocardial microspheres in each frame based on the myocardial microsphere region obtained in each frame, and finding corresponding myocardial microspheres between frames based on the nearest distance between myocardial microspheres between frames, includes: The centroid coordinates are calculated based on the first and zero moments of the myocardial microsphere region contour to obtain the position of all myocardial microspheres in each frame. Using the coordinates of the myocardial microspheres in the first frame as a reference, for each myocardial microsphere, the myocardial microsphere closest to it is found from the myocardial microspheres in other frames.

6. The method according to claim 5, characterized in that, The first and zero moments of the myocardial microsphere region contour and the centroid coordinates are calculated as follows: In the binary image, array(x,y) at the contour is always 1, therefore m 00 It is the number of contour points, m 10 It is the sum of the x-coordinates of all contour points, m 01 It is the sum of the y-coordinates of all contour points, and (X,Y) are the centroid coordinates.

7. A transient analysis device for myocardial microspheres, characterized in that, The myocardial microspheres are obtained by 3D culture of 2D cardiomyocytes; the device includes: A segmentation module is used to acquire a video of myocardial microspheres and to obtain the myocardial microsphere region in each frame of the video using a segmentation method; the step of obtaining the myocardial microsphere region in each frame of the video using a segmentation method includes: separating the myocardial microsphere region from the background in each frame of the video using the Otsu threshold segmentation method; The determination module is used to determine the position of all myocardial microspheres in each frame based on the myocardial microsphere region obtained in each frame, and to find the corresponding myocardial microspheres between frames based on the closest distance between myocardial microspheres between frames. The calculation module is used to select a reference frame and calculate the jitter intensity and jitter velocity using the frame difference method and the maximum value projection method, and generate jitter intensity curve and jitter velocity curve respectively; the jitter curve includes the jitter intensity curve and the jitter velocity curve; The analysis module is used to determine each peak based on the fluctuation curve, calculate the frequency based on the peak, and perform transient analysis on each peak. When the vibration intensity curve shows an upward sloping trend, before determining each peak based on the vibration curve, the method further includes: Based on the bounce intensity curve, the trough is determined, and the trough is used as a reference frame to recalculate the bounce intensity and generate a new bounce intensity curve, specifically as follows: On the jitter intensity curve, if a value is the minimum of the five values ​​before and after it and is less than the threshold, it is a trough. For each frame, the jitter intensity is recalculated and the jitter intensity curve is generated with the nearest trough before that frame as the reference frame.

8. An electronic device, characterized in that, include: One or more processors; A storage unit for storing one or more programs that, when executed by one or more processors, enable the one or more processors to implement the method according to any one of claims 1 to 6.

9. An in vitro method for detecting the regulatory effect of drugs on myocardial microsphere function using a transient analysis method as described in any one of claims 1-6, characterized in that, The in vitro method includes the following steps: In vitro myocardial microspheres were obtained by 3D culture of 2D cardiomyocytes. Injecting drugs into extracorporeal myocardial microspheres; The transient analysis method for myocardial microspheres as described in any one of claims 1-6 was used to obtain the in vitro beating curves of the myocardial microspheres. The functional regulatory effect of the drug on in vitro myocardial microspheres was obtained based on the described oscillation curve analysis.

Images