Cell nuclear membrane regularity calculation method and device, equipment and storage medium

Abstract

The invention discloses a cell nuclear membrane regularity calculation method and device, equipment and a storage medium, and relates to the field of image processing. Sampling and preprocessing the original cell nucleus image data to obtain a complex contour sequence; performing fast Fourier transform and spectral analysis on the contour data of the complex contour sequence, and calculating the frequency component of the complex contour; calculating nuclear membrane regularity index data according to the set Fourier transform frequency threshold; the nuclear membrane regularity index data comprises positive frequency nuclear membrane regularity, negative frequency nuclear membrane regularity and comprehensive nuclear membrane regularity; and performing quantitative evaluation on the regularity of the nuclear membrane based on the nuclear membrane regularity index data, and outputting an evaluation result. According to the scheme, accurate quantitative evaluation on the regularity of the nuclear membrane is realized through a Fourier analysis technology, the limitation of a traditional method is overcome, and powerful technical support is provided for cell pathology diagnosis.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of image processing, in particular to a method and device for calculating the regularity of a cell nucleus membrane, an equipment and a storage medium. BACKGROUND

[0002] In cytopathology diagnosis, the morphological characteristics of the cell nucleus membrane are one of the important indicators for judging the degree of cell malignancy. The nucleus membrane of normal cells is usually regular, while the nucleus membrane of abnormal cells (such as tumor cells) often shows irregular morphological changes.

[0003] Traditional methods for evaluating the regularity of the nucleus membrane mainly include geometric morphological indicators such as circularity, ellipticity, and perimeter-to-area ratio. Although these methods are simple and easy to use, they have obvious limitations: 1. The traditional method for calculating circularity cannot effectively distinguish between regular nucleus membranes of elliptical shape and irregular nucleus membranes of circular shape. For example, a smooth ellipse and a nearly circular contour with local concave-convex may have similar circularity values, but their pathological significance is completely different. 2. The traditional method cannot quantitatively analyze the frequency domain characteristics of the nucleus membrane contour, and cannot separate the different effects of overall shape changes and local curvature changes on regularity. 3. The traditional method lacks intuitive visualization means, and doctors or researchers have difficulty understanding the specific influence of various morphological parameters on the evaluation results of the regularity of the nucleus membrane.

[0004] Therefore, there is an urgent need to develop a new quantitative evaluation method for the regularity of the nucleus membrane that can overcome the above-mentioned defects and provide a more accurate and intuitive quantitative tool for cytopathology diagnosis. SUMMARY

[0005] The present application provides a method and device for calculating the regularity of a cell nucleus membrane, an equipment and a storage medium to solve the technical problems of being unable to distinguish the regularity of the nucleus membrane contour and lacking frequency domain analysis capability in the prior art.

[0006] In one aspect, the present application provides a method for calculating the regularity of a cell nucleus membrane, which comprises: sampling and preprocessing the original cell nucleus image data to obtain a complex contour sequence; performing fast Fourier transform and spectral analysis on the contour data of the complex contour sequence to calculate the frequency components of the complex contour; calculating nucleus membrane regularity index data according to the set Fourier transform frequency threshold; the nucleus membrane regularity index data includes positive frequency nucleus membrane regularity, negative frequency nucleus membrane regularity, and comprehensive nucleus membrane regularity; quantitatively evaluating the regularity of the cell nucleus membrane based on the nucleus membrane regularity index data and outputting the evaluation results.

[0007] Specifically, the original cell nucleus image data is sampled and preprocessed, including: The original cell nucleus image data is sampled according to a sampling frequency, and the sampled contour points are two-dimensional coordinate points ; is a sampling index, is a maximum number of sampling points; The contour centroid coordinates of the nuclear membrane are calculated , and the calculation formula is as follows:

[0008] All contour points are converted into polar coordinate form, and the polar coordinate conversion is as follows: ,

[0009] wherein and respectively represent the angle and length of the contour point from the contour centroid point; All contour points are sorted according to the angle , and equal-interval resampling is performed to generate N standardized contour points; The standardized contour points are normalized to the nuclear membrane contour centroid and converted into the complex contour sequence, and the complex contour point is represented as follows:

[0010] wherein represents an imaginary unit.

[0011] Specifically, the sampling interval of the N standardized contour points is , and the angle of the th standardized contour point after resampling is .

[0012] Specifically, the fast Fourier transform calculates the frequency component , and is represented as follows:

[0013] wherein represents the th frequency component, is a frequency index, and the value range is , is a maximum number of frequency items, is a maximum number of sampling points, represents the th complex contour point, represents an imaginary unit.

[0014] Specifically, the nuclear membrane regularity index data includes: Setting a frequency threshold , calculating the nuclear membrane regularity index according to the frequency threshold; Positive frequency nuclear membrane regularity: ; Negative frequency nuclear membrane regularity: ; Comprehensive nuclear membrane regularity: ; Wherein And is a weight coefficient.

[0015] Specifically, the larger the regularity calculation result is, the more irregular the nuclear membrane is, and the smaller the calculation result is, the more regular the nuclear membrane is; The nuclear membrane regularity index data is used to quantitatively evaluate the regularity of the nuclear membrane, including: Setting a sensitivity threshold; When the calculated nuclear membrane regularity reaches the set sensitivity threshold, the quantitative evaluation result is directly outputted; When the calculated nuclear membrane regularity does not reach the set sensitivity threshold, the profile of the cell nuclear membrane is reconstructed based on the Fourier series.

[0016] Specifically, the profile of the cell nuclear membrane is reconstructed as follows: The frequency component of the profile of the cell nuclear membrane is multiplied by And , to form an Euler circle rotating with time ; It is expressed as follows:

[0017] Rewrite the frequency component as a polar coordinate vector form under the complex plane ; Wherein And are the radius and rotation starting phase of the Euler circle ; By changing the time parameter , the Euler circles are rotated and synthesized to obtain the reconstructed profile of the cell nuclear membrane.

[0018] On the other hand, the present application provides a cell nuclear membrane regularity calculation device, which comprises: A sampling processing module for sampling and preprocessing the original cell nuclear image data to obtain a complex profile sequence; A spectrum analysis module for performing fast Fourier transform and spectrum analysis on the profile data of the complex profile sequence to calculate the frequency component of the complex profile; ​The regularity calculation module is configured to calculate a nuclear membrane regularity index data according to a set Fourier transform frequency threshold; the nuclear membrane regularity index data includes a positive frequency nuclear membrane regularity, a negative frequency nuclear membrane regularity and a comprehensive nuclear membrane regularity. The quantitative evaluation module is configured to quantitatively evaluate the regularity of the nuclear membrane based on the nuclear membrane regularity index data and output an evaluation result.

[0019] In another aspect, the present application provides a computer device, which comprises a processor and a memory, and the memory stores at least one instruction, at least one program, a code set or an instruction set, which is loaded and executed by the processor to implement the nuclear membrane regularity calculation method in the above aspect.

[0020] In another aspect, the present application provides a computer readable storage medium, which stores at least one instruction, at least one program, a code set or an instruction set, which is loaded and executed by a processor to implement the nuclear membrane regularity calculation method in the above aspect.

[0021] The technical scheme provided by the embodiments of the present application has at least the following beneficial effects: by converting the nuclear membrane contour into a complex sequence for frequency domain analysis, combining Fourier transform and spectrum feature calculation to calculate the regularity index, the quantitative differentiation of the overall shape and local curvature change characteristics of the nuclear membrane is realized, which has the advantages of realizing quantitative evaluation of nuclear membrane morphology through frequency domain analysis, effectively distinguishing overall shape change and local irregular characteristics, and providing a visual reconstruction means to improve the diagnostic interpretability. BRIEF DESCRIPTION OF DRAWINGS

[0022] Figure 1 is a flowchart of the nuclear membrane regularity calculation method provided by the embodiments of the present application; Figure 2 is an algorithm flowchart of the nuclear membrane regularity calculation method provided by the embodiments of the present application; Figure 3 shows the reconstruction of the nuclear membrane contour with the maximum number of frequency terms shows the Fourier analysis spectrum distribution graph in the scene; Figure 4 shows the nuclear membrane contour graph directly generated by Fourier transform; Figure 5 lists the schematic diagram of constructing the Euler circle rotation trajectory to reconstruct the nuclear membrane contour; Figure 6 shows the schematic diagram of reconstructing the nuclear membrane contour with 5 significant frequency components; Figure 7 shows the schematic diagram of reconstructing the nuclear membrane contour with 15 significant frequency components; Figure 8 A schematic diagram showing dynamic demonstration of cell nucleus membrane profile with time and Euler circle number is shown. Figure 9 A schematic diagram showing dynamic demonstration of cell nucleus membrane profile with frequency threshold and significant frequency component is shown. Figure 10 A schematic diagram showing typical normal cell nucleus membrane profile analysis is shown. Figure 11 A schematic diagram showing typical abnormal cell nucleus membrane profile analysis is shown. Figure 12 A structural block diagram of a cell nucleus membrane regularity calculation device provided by an embodiment of the present application is shown. DETAILED DESCRIPTION

[0023] In order to make the purpose, technical scheme and advantages of the present application clearer, the embodiments of the present application will be further described in detail below with reference to the drawings.

[0024] “Multiple” mentioned in the present application refers to two or more. “And / or” describes the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B can represent the three cases of A alone, A and B together, and B alone. The character “ / ” generally represents that the associated objects before and after it are in an “or” relationship.

[0025] Figure 1 A flowchart of a cell nucleus membrane regularity calculation method provided by an embodiment of the present application is shown, which includes the following steps: S1, original cell nucleus image data is sampled and preprocessed to obtain a plurality of profile sequences; The sampling and preprocessing are to simplify the original cell nucleus image data for fast calculation and visual display. The preprocessing process involves dimensioning the sampling point data according to the coordinates, and forming a plurality of profile sequences from all the sampling points. The plurality of profile sequences refer to a set of profile points converted from two-dimensional space coordinates to complex numbers, which can specifically use the center point or the centroid point of each profile point to establish a coordinate system. This conversion converts the geometric shape into a complex signal sequence that can be analyzed in the frequency domain, establishing a mathematical basis for subsequent Fourier transform.

[0026] S2, the profile data of the plurality of profile sequences is subjected to fast Fourier transform and frequency spectrum analysis to calculate the frequency components of the plurality of profile sequences; The fast Fourier transform (FFT) refers to a process of performing a discrete Fourier transform on a complex sequence and calculating frequency components. Through the FFT algorithm, the contour data of the cell can be converted from the time domain to the frequency domain to obtain a series of frequency components. After performing amplitude spectrum analysis on these frequency components, the energy contribution of each frequency component can be evaluated, and the nuclear membrane regularity can be analyzed in a targeted and directional manner. Specifically, the fast Fourier transform algorithm can be used to achieve this, and the low-frequency components corresponding to the overall shape and the high-frequency components corresponding to the local details can be separated through spectral analysis.

[0027] S3, calculating the nuclear membrane regularity index data according to the set Fourier transform frequency threshold; the nuclear membrane regularity index data includes positive frequency nuclear membrane regularity, negative frequency nuclear membrane regularity, and comprehensive nuclear membrane regularity. The Fourier transform frequency threshold refers to the critical value of the frequency component used in the regularity calculation. Since low-frequency components mainly reflect the circularity of the nuclear membrane and do not reflect the regularity of the nuclear membrane, the threshold needs to be set to reduce the influence of circularity in the regularity calculation. Specifically, the preset frequency index value or the dynamic calculation of significant frequency components can be used to determine the threshold, and the frequency range that plays a dominant role in the overall shape and local details can be distinguished by setting the threshold.

[0028] The design of the nuclear membrane regularity index is used to judge the irregular state of the nuclear membrane contour. In the present application, the nuclear membrane regularity index data includes one or more combinations of positive frequency, negative frequency, and comprehensive regularity. When calculating the positive or negative frequency parameters, the high-frequency components exceeding the threshold can be accumulated and weighted summed to achieve this. The positive and negative frequencies correspond to the clockwise and counterclockwise rotating parts of the nuclear contour, respectively. The low-frequency nuclear membrane regularity reflects the overall shape deviation of the nuclear membrane contour, and the high-frequency nuclear membrane regularity corresponds to the local curvature change. The comprehensive nuclear membrane regularity integrates the high and low frequency information in the positive and negative frequencies through a weight coefficient, which reflects the irregularity of the cell nuclear membrane and provides a comprehensive evaluation that takes into account the overall shape and local details.

[0029] In some schemes, using a single index data may cause accuracy bias, so in the ideal state, the scheme of using multiple combined evaluations is optimal.

[0030] S4, quantitatively evaluating the regularity of the cell nuclear membrane based on the nuclear membrane regularity index data and outputting the evaluation result.

[0031] Quantitative evaluation refers to outputting a numerical evaluation result based on the regularity index. Specifically, the threshold setting comparison or the generation of a visualized reconstructed contour can be used to achieve this. The quantitative result can be used to objectively judge the nuclear membrane regularity and support pathological diagnosis decisions. In the present application, the focus is on calculating or evaluating whether the regularity of the nuclear membrane contour meets the evaluation standard based on the nuclear membrane regularity index data parameters. The specific schemes for evaluating cell health or tumor cells belong to the biomedical field, and the embodiments of the present application do not limit the description of this.

[0032] Because it is the way to extract high frequency and low frequency components using fast Fourier transform to evaluate the regularity of the nuclear membrane. And the greater the regularity calculation result indicates the more irregular the nuclear membrane, the smaller the calculation result the more regular the nuclear membrane. In some special scenarios, when the original cell nucleus image data sampling and preprocessing is not enough, or the preliminary Fourier transform calculation data cannot effectively separate the overall and local detail feature components, the evaluation result may affect the medical diagnosis. To avoid this phenomenon, the present application can also introduce the nuclear standard, that is, whether the sensitivity of the nuclear membrane irregularity is satisfied, the evaluation standard or sensitivity threshold of the nuclear membrane regularity index data can be set, or the regularity index is compared with the pre-established database to determine whether the contour reconstruction is needed.

[0033] Figure 2 The algorithm flow chart of the nuclear membrane regularity calculation method provided by the present application is to set the threshold value. When the calculated nuclear membrane regularity reaches the set sensitivity threshold, the quantitative evaluation result is directly output according to the related parameters and rules; when the calculated nuclear membrane regularity does not reach the set sensitivity threshold, the frequency components extracted in the foregoing are combined to reconstruct the contour of the nuclear membrane, and then the frequency components participating in the reconstruction are adjusted according to the requirements to recalculate the nuclear membrane regularity index data until the evaluation result is output after reaching a satisfactory degree.

[0034] In summary, the present application converts the two-dimensional contour problem into a one-dimensional complex sequence analysis problem through the original cell nucleus image data sampling and preprocessing, lays a foundation for subsequent Fourier analysis; through the frequency domain decomposition, the nuclear membrane contour is converted into a complex signal for Fourier analysis, a multi-dimensional evaluation system of positive and negative frequency component energy accumulation is constructed, and the whole shape and local details are dynamically separated combined with the frequency threshold, which is helpful for researchers to directional analysis. Through this method of frequency domain decomposition and multi-dimensional index construction, the technical solution effectively solves the problem that the traditional method cannot distinguish the irregularity degree of different types of nuclear membrane contour, and realizes the accurate quantification of the smoothness and concave-convex situation of the nuclear membrane.

[0035] In some embodiments, there is a problem of uneven distribution of contour point angles in the original contour data sampling process, and direct Fourier transform will cause distortion of the frequency spectrum analysis result. Non-equidistant angle distribution will introduce high-frequency noise interference, affecting the accuracy of the nuclear membrane regularity index. For this, the present application can sample and polar coordinate convert the contour according to the sampling frequency, sort according to the angle size and equidistantly resample to generate standardized contour points. The specific process can be realized through the following steps: A, sample the contour according to the sampling frequency of the original cell nucleus image data, and the sampled contour points are two-dimensional coordinate points ; is the sampling index, For the maximum number of sampling points; B, calculate the centroid coordinates of the nuclear membrane contour The calculation formula is as follows:

[0036] C, convert all contour points to polar coordinates, and the polar coordinate conversion is as follows: ,

[0037] Among them and respectively represent the angle and length of the contour point from the contour centroid point; D, sort all contour points according to the angle Size and carry out equidistant resampling to generate N standard contour points; This step is mainly to ensure the continuity of the contour, and the sorted expression is . Assuming that it is N standard contour points, the sampling interval is , and the angle of the th standard contour point after resampling is .

[0038] Since the resampled contour points are evenly distributed in angle, the contour distortion caused by inconsistent angle intervals is effectively avoided, thereby improving the accuracy and reliability of the calculation of the regularity of the nuclear membrane.

[0039] E, normalize the standard contour points to the nuclear membrane contour centroid and convert them to complex contour sequences, and the complex contour point is expressed as follows:

[0040] Among them is the imaginary unit.

[0041] For example, when the number of sampling points is 128, the angular interval of each standardized contour point is set to 2π / 128 ≈ 0.049 radians. During resampling, the radial distance values ​​of the original contour points are mapped to new equiangular positions using an interpolation algorithm. The normalized coordinate data is converted into complex form, with the real part corresponding to the x-axis offset and the imaginary part corresponding to the y-axis offset. This data processing method ensures that the subsequent Fourier transform accurately reflects the periodic characteristics of the contour shape, avoiding spectral leakage caused by uneven angular distribution. By establishing a standardized coordinate system based on the centroid, the influence of the positional differences of different cell nuclei on the analysis results is effectively eliminated, ensuring that the regularity index only reflects shape characteristics. At the same time, the normalization to the centroid eliminates the influence of spatial translation components on the complex sequence, enabling the Fourier transform results to more accurately reflect the contour shape characteristics. This standardization process provides input data that meets the requirements of the Fourier transform for subsequent frequency domain analysis, thereby improving the accuracy and reliability of the nuclear membrane regularity index calculation.

[0042] In some embodiments, the steps of performing fast Fourier transform and spectral analysis on the contour data of the complex contour sequence can be achieved through the following steps:

[0043]

[0044] Among them Indicates the first One frequency component, Indicates the amplitude component. This is a frequency index, with a value range of [value range missing]. , The number of terms with the highest frequency. The maximum number of sampling points, Indicates the first A complex number of contour points, It represents the imaginary unit.

[0045] In some embodiments, a maximum number of frequency terms is set to ensure computational stability. , among them As a significance frequency component threshold, in the application scenario of cell nuclear membrane contour, the preferred value is... By The value of is restricted to Within the range, the effective frequency range corresponding to the main morphological features of the nuclear membrane contour can be covered, while excluding high-frequency noise components that exceed this range.

[0046] In general, the maximum number of sampling points is preferred. Taking 60 as an example, this is the number of terms with the highest frequency. For example, when Time, At this time, the amplitude spectrum values of the first 5 frequency components are significantly higher than the subsequent components, indicating that the main energy is concentrated in the low frequency band, corresponding to the smooth profile characteristics of the regular nuclear membrane.

[0047] The amplitude component converts the complex frequency component into a real energy value, facilitating subsequent statistics and comparison of regularity indicators. Figure 3 The maximum number of frequency terms is shown The Fourier analysis spectrum distribution diagram under the scene, the first item amplitude is the highest 42.0159.

[0048] In the embodiments of the present application, for the calculation of the nuclear membrane regularity index data, the following method can be used for calculation: Set the frequency threshold , calculate the nuclear membrane regularity index according to the frequency threshold; Positive frequency nuclear membrane regularity: ; Negative frequency nuclear membrane regularity: ; Comprehensive nuclear membrane regularity: ; Wherein and are weight coefficients. The positive / negative frequency nuclear membrane regularity weight coefficients , are initialized to 1, and then set according to experience, which is not limited in the present application.

[0049] Specifically, by setting the frequency threshold , the low-frequency components representing the overall shape can be effectively filtered out, and the high-frequency components reflecting the local details can be retained. In the calculation process of the positive frequency nuclear membrane regularity , the amplitude spectrum of the positive frequency component is accumulated to quantify the cumulative intensity of the local curvature change in the counterclockwise direction of the nuclear membrane; the calculation of the negative frequency nuclear membrane regularity accumulates the amplitude spectrum of the negative frequency component to represent the cumulative intensity of the local curvature change in the clockwise direction. The comprehensive nuclear membrane regularity weights and fuses the local curvature changes in the positive and negative directions through the weight coefficients. This technical solution, combined with the pre-generated complex profile sequence method, uses the frequency components obtained by fast Fourier transform as the calculation basis, and realizes the accurate quantification of the local regularity of the nuclear membrane through threshold screening and direction separation.

[0050] In some embodiments, when the profile regularity constructed by preliminary Fourier transform does not meet the requirements, the reconstruction program needs to be started. Figure 4 The nuclear membrane profile generated directly by Fourier transform is shown, and its reconstruction process can be summarized as follows: A, the frequency component of the nuclear membrane profile With multiplication, constitute the time rotating Euler circle ; as follows:

[0051] B, rewrite the frequency component as the polar coordinate vector form under the complex plane ; Where and is the radius and rotation start phase of the Euler circle ; C, by the change of time parameter , the Euler circle rotation synthesis, get the reconstructed nuclear membrane contour trajectory.

[0052] The steps of reconstructing the nuclear membrane contour based on Fourier series can include various implementations. For example, the frequency component can be converted to polar coordinate vector form, the radius and rotation start phase of which are determined by the amplitude and phase angle. During the reconstruction synthesis process, the Euler circle corresponding to the positive frequency component (i.e. the Euler circle of ) can be set to rotate counterclockwise, and the Euler circle corresponding to the negative frequency component (i.e. the Euler circle of ) can be set to rotate clockwise, and the difference in the rotation direction of the two can reflect the physical meaning of high-frequency fluctuations and low-frequency shape changes. The dynamic change of the center position can be realized by vector superposition of the previous frequency component, for example, when processing the nth frequency component, its center position can be located at the end of the vector corresponding to the (n-1)th frequency component (i.e. the center of the circle of is located at the end of the polar coordinate vector corresponding to ). The value range of the time parameter t can be set to 0 to 2π, and the continuous rotation of the Euler circle is driven by the parameter change to synthesize the complete contour trajectory. The reconstructed contour trajectory can be displayed as a dynamic image, which is convenient for observing the superimposed effect of different frequency components on the overall shape.

[0053] Figure 5 The nuclear membrane contour reconstructed based on the Euler circle rotation trajectory constructed by the frequency component is listed, and by comparison with Figure 4 , it can be seen that the reconstruction retains the characteristics of the nuclear membrane contour, and each frequency component corresponds to a rotating Euler circle , and the superimposed trajectory of all Euler circles is the reconstructed nuclear membrane contour. The reconstructed contour pays more attention to local contour details and is more accurate in quantification.

[0054] In some embodiments, there is a problem in the reconstruction process that cannot intuitively show the dynamic influence of different frequency components on the contour shape, which may make it difficult for the operator to understand the difference in the contribution of high-frequency components and low-frequency components to the local regularity and overall shape change of the nuclear membrane, thereby affecting the interpretability of the regularity evaluation results.

[0055] To this end, the present application can target the significant frequency components by setting the maximum number of frequency terms wherein is the significant frequency component threshold. This process is mainly screened by amplitude size, for example, the present application sets Based on the determination of the amplitude component , the normalization processing is continued, which is expressed as follows:

[0056] Further, based on the significant frequency component threshold and the standard amplitude spectrum value, the significant frequency component is determined. The significant frequency component is an adjustable parameter item for the visualization interface display, and specifically an interactive control interface can be provided, and the operator can manually adjust the frequency component participating in the reconstruction of the nuclear membrane contour, and then observe the influence of each frequency component on the overall roundness and local change of the visual generated nuclear membrane contour, thereby determining the frequency threshold range.

[0057] Figure 6 and Figure 7 respectively show the schematic diagram of reconstructing the nuclear membrane contour with 5 and 15 significant frequency components, wherein the number of contour points , the frequency threshold is taken as an example for illustration. Figure 6 The black solid line in the middle is the original nuclear membrane contour, and the gray dashed line is the reconstructed nuclear membrane contour. As can be seen from the comparison, Figure 6 the low-frequency component mainly affects the overall shape, Figure 7 the high-frequency component mainly affects the local details. As the number of high-frequency components increases, the reconstructed nuclear membrane contour is closer to the input nuclear membrane contour.

[0058] As shown in Figure 8 , in some embodiments, the visualization interface can also set the time and the number of Euler circles to realize dynamic demonstration. Figure 9 Also shown is the setting of the frequency threshold and the number of significant frequency components (i.e. the number of Fourier terms), which facilitates the dynamic adjustment of the number of significant frequency components and the frequency threshold by the operator. Through the intuitive user interface, including spectrum display, Euler circle animation, parameter control and result comparison, etc. Functional modules. The operator can adjust the frequency threshold in real time through the slider to observe the change of the contour reconstruction effect, thereby determining the optimal analysis parameters.​

[0059] The following is illustrated with specific embodiments: Example 1: Analysis of normal nuclear membrane The profile of a normal cell's nuclear membrane was analyzed. 105 profile points were sampled, and the maximum frequency term number M was set to 20. The calculation result showed that the main energy of the profile was concentrated in the first 5 frequency components, and the amplitude of the high-frequency components was very small. The frequency threshold was set to , and the comprehensive regularity degree = 0.513579 was obtained, indicating that the nuclear membrane was relatively regular.

[0060] Example 2: Analysis of abnormal nuclear membrane The profile of an abnormal cell's nuclear membrane was analyzed. 120 profile points were sampled, and M was set to 20. The calculation result showed that there was significant energy distribution in the high-frequency band, indicating that the nuclear membrane profile had obvious irregularity. The frequency threshold was set to , and the comprehensive regularity degree was obtained, which was significantly higher than that of the normal cell, indicating that the nuclear membrane was highly irregular.

[0061] As shown in Figure 10 and Figure 11 , through comparative analysis, it can be seen that the method of the present application can effectively distinguish the regularity difference between normal and abnormal cell nuclear membranes, realize accurate quantitative evaluation of the regularity of the cell nuclear membrane, and overcome the limitations of traditional methods, providing a quantitative basis for cell pathology diagnosis.

[0062] In summary, the present application has the following beneficial technical effects: (1) By Fourier analysis, the cell nuclear membrane profile is converted from spatial domain to frequency domain, which can effectively separate the overall shape change and local detail change, overcoming the problem that traditional methods cannot distinguish different types of nuclear membrane regularity; (2) A positive frequency, negative frequency and comprehensive nuclear membrane regularity index system is established, providing a multi-dimensional regularity evaluation scheme, improving the accuracy and reliability of the evaluation results; (3) The Euler circle animation visualizes the Fourier series reconstruction process, allowing researchers to intuitively understand the influence of each frequency component on the nuclear membrane shape, facilitating parameter adjustment and result interpretation; (4) The interactive interface allows real-time adjustment of the frequency threshold, allowing researchers to balance the weight of overall shape and local details in regularity evaluation according to specific needs; (5) This method is not only suitable for cell nuclear membrane analysis, but also can be applied to other biological morphology analysis fields, and has wide application prospects.

[0063] Figure 12A structural block diagram of a nuclear membrane regularity calculation device provided by an embodiment of the present application is shown, and the device comprises: A sampling processing module 1210 is configured to sample and pre-process the original nuclear image data to obtain a plurality of contour sequences. A spectrum analysis module 1220 is configured to perform fast Fourier transform and spectrum analysis on the contour data of the contour sequences to calculate frequency components of the complex number contours. A regularity calculation module 1230 is configured to calculate nuclear membrane regularity index data according to a set Fourier transform frequency threshold, wherein the nuclear membrane regularity index data comprises positive frequency nuclear membrane regularity, negative frequency nuclear membrane regularity and comprehensive nuclear membrane regularity. A quantitative evaluation module 1240 is configured to quantitatively evaluate the regularity of the nuclear membrane based on the nuclear membrane regularity index data and output an evaluation result.

[0064] It should be noted that the nuclear membrane regularity calculation device provided in the embodiments of the present application is only exemplified by the division of the above functional modules / functional units, and in actual application, the above functions can be completed by different functional modules / functional units according to needs, that is, the internal structure of the nuclear membrane regularity calculation device is divided into different functional modules / functional units to complete all or part of the functions described above. In addition, the implementation of the nuclear membrane regularity calculation method provided by the method embodiments belongs to the same concept as the implementation of the nuclear membrane regularity calculation device provided by the present embodiment, and the specific implementation process of the nuclear membrane regularity calculation device provided by the present embodiment is described above in the method embodiments, which will not be repeated here.

[0065] The specific embodiments are only an explanation of the present application, and are not a limitation of the present application. Those skilled in the art can make non-creative modifications to the embodiments according to needs after reading the present specification, but as long as the modifications are within the scope of the claims of the present application, they are protected by the patent law.

Claims

1. A method for calculating the regularity of a cell nuclear membrane, characterized in that, The method includes: The raw cell nucleus image data is sampled and preprocessed to obtain complex contour sequences; The contour data of the complex contour sequence are subjected to fast Fourier transform and spectral analysis to calculate the frequency components of the complex contour. Nuclear membrane regularity index data are calculated based on the set Fourier transform frequency threshold; the nuclear membrane regularity index data includes positive frequency nuclear membrane regularity, negative frequency nuclear membrane regularity, and comprehensive nuclear membrane regularity; The regularity of the cell nuclear membrane is quantitatively evaluated based on the nuclear membrane regularity index data, and the evaluation results are output.

2. The method according to claim 1, characterized in that, The process of sampling and preprocessing the raw cell nuclear image data includes: The original cell nucleus image data is sampled for contours according to a sampling frequency, and the sampled contour points are two-dimensional coordinate points. ; For sampling index, This represents the maximum number of sampling points. Calculate the centroid coordinates of the nuclear membrane profile The calculation formula is as follows: Convert all contour points to polar coordinates. The polar coordinate transformation is as follows: ， Among them and These represent the angle and length of the distance between the contour point and the centroid of the contour, respectively. All contour points according to angle The sizes are sorted and resampled at equal intervals to generate N standardized contour points; The standardized contour points are normalized to the centroid of the nuclear membrane contour and converted into the complex contour sequence. It is expressed as follows: Among them It represents the imaginary unit.

3. The method according to claim 2, characterized in that, Sampling interval of N standardized contour points After resampling Angle of a standardized contour point .

4. The method according to claim 1, characterized in that, The Fast Fourier Transform (FFT) is used to calculate the spectral components, as shown below: Among them Indicates the first One frequency component, Indicates the amplitude component, This is a frequency index, with a value range of [value range missing]. , The number of terms with the highest frequency. The maximum number of sampling points, Indicates the first A complex number of contour points, It represents the imaginary unit.

5. The method according to claim 4, characterized in that, The data used to calculate the nuclear membrane regularity index includes: Set frequency threshold The nuclear membrane regularity index is calculated based on the frequency threshold. Positive frequency nuclear membrane regularity: ; Negative frequency nuclear membrane regularity: ; Overall nuclear membrane regularity: ; Among them is and These are the weighting coefficients.

6. The method according to claim 1, characterized in that, A higher regularity result indicates a more irregular nuclear membrane, while a lower result indicates a more regular nuclear membrane. The quantitative assessment of the regularity of the cell nuclear membrane based on the nuclear membrane regularity index data includes: Set a sensitivity threshold; When the calculated nuclear membrane regularity reaches the set sensitivity threshold, the quantitative evaluation result is directly output; When the calculated nuclear membrane regularity does not reach the set sensitivity threshold, the cell nuclear membrane contour is reconstructed based on Fourier series.

7. The method according to claim 6, characterized in that, The steps for reconstructing the nuclear membrane outline are as follows: Frequency components of the cell nuclear membrane profile and Multiplication, forming a product over time Euler circle in rotation ; indicates the following: Rewrite the frequency components in polar coordinate vector form in the complex plane. ;in and For Euler circle The radius and the initial phase of rotation; By time parameters The changes will By synthesizing Euler circles, the reconstructed nuclear membrane contour trajectory is obtained.

8. A device for calculating the regularity of a cell nuclear membrane, characterized in that, The device includes: The sampling and processing module is used to sample and preprocess the raw cell nucleus image data to obtain complex contour sequences. The spectrum analysis module is used to perform fast Fourier transform and spectrum analysis on the contour data of the complex contour sequence to calculate the frequency components of the complex contour. The regularity calculation module is used to calculate nuclear membrane regularity index data based on the set Fourier transform frequency threshold; the nuclear membrane regularity index data includes positive frequency nuclear membrane regularity, negative frequency nuclear membrane regularity, and comprehensive nuclear membrane regularity; The quantitative evaluation module is used to quantitatively evaluate the regularity of the cell nuclear membrane based on the nuclear membrane regularity index data and output the evaluation results.

9. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one instruction, at least one program, code set, or instruction set, the at least one instruction, the at least one program, the code set, or the instruction set being loaded and executed by the processor to implement the cell nuclear membrane regularity calculation method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The readable storage medium stores at least one instruction, at least one program, code set, or instruction set, wherein the at least one instruction, the at least one program, the code set, or instruction set is loaded and executed by a processor to implement the cell nuclear membrane regularity calculation method as described in any one of claims 1 to 7.

Images