Method and system for analyzing vertebral compression fracture state, electronic equipment and medium

Abstract

The invention provides a centrum compression fracture state analysis method and system, electronic equipment and a medium. The analysis method comprises the steps that a target centrum image is acquired; segmenting and extracting each target vertebral body in the target vertebral body image; obtaining a target position and a target theoretical area corresponding to each target vertebral body; calculating to obtain an area difference value of theoretical areas of a target vertebral body at the current target position and an adjacent vertebral body; and based on the area difference value, determining whether the current target vertebral body is subjected to integral compression or not. On the basis of segmentation and extraction of each cone in a cone image, area associated data of the cones are considered, and the area difference condition of adjacent cones is calculated to determine whether the current target cone is subjected to integral compression or not. The precision, efficiency and reliability of a determination result of whether the centrum in any target centrum image is subjected to integral compression or not output by the centrum compression fracture analysis model are effectively ensured.

Description

Technical Field

[0001] The present disclosure relates to the field of image processing technologies, and particularly to a method, a system, an electronic device, and a medium for analyzing the state of vertebral compression fractures. Background Art

[0002] Currently, for the diagnosis of vertebral compression fractures, X-ray films are generally used, with low detection costs and faster detection speeds. Clinically, when identifying compression fractures, it is necessary to evaluate and determine the location and degree of the fractures. However, currently, the treatment plans for vertebral (such as thoracolumbar) compression fractures generally have problems such as complex treatment processes and low accuracy. Summary of the Invention

[0003] The technical problem to be solved by the present disclosure is to overcome the above-mentioned defects in the prior art and provide a method, a system, an electronic device, and a medium for analyzing the state of vertebral compression fractures.

[0004] The present disclosure solves the above technical problems through the following technical solutions:

[0005] The present disclosure provides a method for analyzing the state of vertebral compression fractures, and the analysis method includes:

[0006] Obtain a target vertebral image;

[0007] Segment and extract each target vertebra in the target vertebral image;

[0008] Obtain the target position and the target theoretical area corresponding to each target vertebra;

[0009] Calculate the area difference between the target vertebra at the current target position and the theoretical areas of adjacent vertebrae;

[0010] Based on the area difference, determine whether the current target vertebra has undergone overall compression.

[0011] Optionally, the step of determining whether the vertebra in the target vertebral image has undergone overall compression based on the area difference includes:

[0012] In response to the area difference not falling within the preset range representing the area difference between normal vertebrae, determine that the current target vertebra has undergone overall compression;

[0013] In response to the area difference falling within the preset range, determine that the target vertebra has not undergone overall compression.

[0014] Optionally, the analysis method further includes:

[0015] Determine an optimal area difference threshold from the preset range for determining whether the vertebral bodies at different positions have undergone overall compression;

[0016] The step of determining whether the vertebral body in the target vertebral body image has undergone overall compression based on the area difference includes:

[0017] In response to the area difference being greater than the optimal area difference threshold, determine that the current target vertebral body has undergone overall compression;

[0018] In response to the area difference being less than or equal to the optimal area difference threshold, determine that the target vertebral body has not undergone overall compression.

[0019] Optionally, the analysis method further includes:

[0020] In response to the current target vertebral body not having undergone overall compression, obtain the target true area corresponding to the current target vertebral body;

[0021] Calculate a first ratio of the target true area and the target theoretical area corresponding to each target vertebral body;

[0022] Determine the fracture classification of the current target vertebral body based on the first ratio;

[0023] Wherein, different fracture classifications correspond to different preset ratio ranges.

[0024] Optionally, the analysis method further includes:

[0025] In response to the current target vertebral body having undergone overall compression, use a preset correction scheme to determine the fracture classification of the target vertebral body.

[0026] Optionally, the step of, in response to the current target vertebral body having undergone overall compression, using a preset correction scheme to determine the fracture classification of the target vertebral body includes:

[0027] In response to the current target vertebral body having undergone overall compression, obtain a first theoretical area of the previous vertebral body of the current target vertebral body and a second theoretical area of the next vertebral body;

[0028] Use the weighted average of the first theoretical area and the second theoretical area as the new theoretical area of the current target vertebral body;

[0029] Obtain the target true area of the current target vertebral body;

[0030] Calculate a second ratio of the target true area and the new theoretical area of the current target vertebral body;

[0031] Determine the fracture degree of the current target vertebral body based on the second ratio;

[0032] Among them, different fracture degrees correspond to different preset ratio ranges;

[0033] Or,

[0034] The step of determining the fracture degree of the target vertebral body by adopting a preset correction scheme in response to the overall compression of the current target vertebral body includes:

[0035] In response to the overall compression of the current target vertebral body, obtain the target true area and the target theoretical area of the current target vertebral body;

[0036] Calculate the third ratio of the target true area and the new theoretical area of the current target vertebral body, and multiply the third ratio by a preset coefficient to obtain a fourth ratio;

[0037] Determine the fracture degree of the current target vertebral body based on the fourth ratio;

[0038] Among them, different fracture degrees correspond to different preset ratio ranges.

[0039] Optionally, the step of calculating the target true area of the target vertebral body includes:

[0040] Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain a true edge line corresponding to each edge;

[0041] Based on each true edge line, generate a true edge graph corresponding to the vertebral body;

[0042] Obtain the area corresponding to the true edge graph and use it as the target true area of the target vertebral body;

[0043] And / or,

[0044] The step of calculating the target theoretical area of the target vertebral body includes:

[0045] Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain a true edge line corresponding to each edge;

[0046] Among them, the true edge line includes a true leading edge line, a true trailing edge line, a true upper edge line, and a true lower edge line;

[0047] Based on each true edge line, form a closed graph and use it as the true edge graph corresponding to the vertebral body;

[0048] Generate a circumscribed rectangle of the true edge pattern based on the true leading edge line, the true trailing edge line, and the true edge pattern, and use the circumscribed rectangle as the theoretical shape of the target vertebral body;

[0049] Obtain the area corresponding to the theoretical shape and use it as the target theoretical area of the target vertebral body.

[0050] The present disclosure also provides an analysis system for the state of vertebral compression fractures, the analysis system including:

[0051] A target image acquisition module for acquiring a target vertebral body image;

[0052] A target vertebral body segmentation module for segmenting and extracting each target vertebral body in the target vertebral body image;

[0053] A vertebral body information acquisition module for acquiring the target position and the target theoretical area corresponding to each target vertebral body;

[0054] An area difference calculation module for calculating the area difference between the theoretical areas of the target vertebral body and adjacent vertebral bodies at the current target position;

[0055] A state analysis module for determining whether the current target vertebral body has undergone overall compression based on the area difference.

[0056] The present disclosure also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and used to run on the processor. When the processor executes the computer program, the analysis method for the state of vertebral compression fractures as described above is implemented.

[0057] The present disclosure also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the analysis method for the state of vertebral compression fractures as described above is implemented.

[0058] The present disclosure also provides a computer program product, including a computer program. When the computer program is executed by a processor, the analysis method for the state of vertebral compression fractures as described above is implemented.

[0059] On the basis of conforming to the common knowledge in the art, the above preferred conditions can be combined arbitrarily to obtain various preferred examples of the present disclosure.

[0060] The positive and progressive effects of the present disclosure are as follows:

[0061] In the present disclosure, based on the segmentation and extraction of each vertebral body in the vertebral body image, considering the area correlation data of the vertebral body, the area difference between adjacent vertebral bodies is calculated to determine whether the current target vertebral body has undergone overall compression, effectively ensuring the accuracy, efficiency, and reliability of the determination result of whether the vertebral body in any target vertebral body image has undergone overall compression by the vertebral body compression fracture analysis model, and having the advantages of convenient analysis process and small calculation amount. Description of the Drawings

[0062] Figure 1 It is a flowchart of the method for analyzing the vertebral body compression fracture state in Embodiment 1 of the present disclosure;

[0063] Figure 2 It is the first flowchart of the method for analyzing the vertebral body compression fracture state in Embodiment 2 of the present disclosure;

[0064] Figure 3 It is the second flowchart of the method for analyzing the vertebral body compression fracture state in Embodiment 2 of the present disclosure;

[0065] Figure 4 It is a flowchart of obtaining the target true area of the target vertebral body in Embodiment 2 of the present disclosure;

[0066] Figure 5 It is the original schematic diagram of the vertebral body image in Embodiment 2 of the present disclosure;

[0067] Figure 6 It is a schematic diagram of the actual graph of the vertebral body image in Embodiment 2 of the present disclosure;

[0068] Figure 7 It is a schematic diagram of the posterior edge of the vertebral body image in Embodiment 2 of the present disclosure;

[0069] Figure 8 It is a schematic diagram of the presence of bilateral edges at the upper and lower edges of the vertebral body image in Embodiment 2 of the present disclosure;

[0070] Figure 9 It is a schematic diagram of the midlines at the upper and lower edges of the vertebral body image in Embodiment 2 of the present disclosure;

[0071] Figure 10 It is a schematic diagram of the anterior edge of the vertebral body image in Embodiment 2 of the present disclosure;

[0072] Figure 11 It is a flowchart of obtaining the target theoretical area of the target vertebral body in Embodiment 2 of the present disclosure;

[0073] Figure 12 It is a module schematic diagram of the vertebral body compression fracture state analysis system in Embodiment 3 of the present disclosure;

[0074] Figure 13Schematic diagram of the modules of the vertebral compression fracture status analysis system according to Embodiment 4 of the present disclosure;

[0075] Figure 14 Schematic diagram of the structure of the electronic device according to Embodiment 5 of the present disclosure. Specific embodiments

[0076] The present disclosure will be further described below by way of embodiments, but the present disclosure is not limited to the scope of the described embodiments.

[0077] In the embodiments of the present disclosure, prefix words such as "first" and "second" are only used to distinguish different described objects, and have no limiting effect on the position, order, priority, quantity or content of the described objects. The use of ordinal numbers and other prefix words for distinguishing described objects in the embodiments of the present disclosure does not constitute a limitation on the described objects. The statement of the described objects refers to the description in the claims or the context of the embodiments, and should not constitute an unnecessary limitation because of the use of such prefix words. In addition, in the description of this embodiment, unless otherwise specified, the meaning of "a plurality" is two or more.

[0078] Embodiment 1

[0079] As Figure 1 shown, the method for analyzing the vertebral compression fracture status in this embodiment includes:

[0080] S101. Obtain a target vertebral image;

[0081] Among them, the target vertebral image can be a lateral radiograph of the thoracic and lumbar vertebrae generated by an X-ray machine.

[0082] S102. Segment and extract each target vertebra in the target vertebral image;

[0083] Among them, the target vertebra is other vertebrae except the basic vertebra (that is, the vertebra other than the first vertebra), and the first vertebra is a non-integrally compressed vertebra;

[0084] Specifically, image segmentation technology (such as the yolo-xray image segmentation algorithm, etc.) is used to segment the target vertebral image to obtain each vertebra (that is, each vertebra) in the target vertebral image, so as to ensure the accurate segmentation of the vertebra, and then ensure the accuracy of the determination of parameters such as the theoretical area and the actual area of each vertebra, so as to ensure the accuracy and reliability of the analysis results of the vertebral compression fracture status.

[0085] S103. Obtain the target position and the target theoretical area corresponding to each target vertebra;

[0086] Among them, the target position is the position of each vertebral body in the spine, such as which vertebral body; for example, different vertebral bodies from top to bottom are: T4 - T12 - L1 - L4.

[0087] S104. Calculate the area difference between the target vertebral body at the current target position and the adjacent vertebral body; preferably, the adjacent vertebral body is the previous adjacent vertebral body of the current target vertebral body.

[0088] S105. Based on the area difference, determine whether the current target vertebral body has undergone overall compression.

[0089] In the present disclosure, based on the segmentation and extraction of each vertebral body in the vertebral body image, considering the area correlation data of the vertebral body, calculating the area difference situation between adjacent vertebral bodies to determine whether the current target vertebral body has undergone overall compression effectively ensures the accuracy, efficiency, and reliability of the determination result of whether the vertebral body in any target vertebral body image has undergone overall compression in the vertebral body compression fracture analysis model, and has the advantages of convenient analysis process and small calculation amount.

[0090] Embodiment 2

[0091] The analysis method for the state of vertebral body compression fracture in this embodiment is a further improvement of Embodiment 1. Specifically:

[0092] In an implementable solution, as Figure 2 shown, step S105 includes:

[0093] S1051. In response to the area difference not falling within the preset range representing the area difference between normal vertebral bodies, determine that the current target vertebral body has undergone overall compression;

[0094] In response to the area difference falling within the equal preset range, determine that the area target vertebral body has not undergone overall compression.

[0095] Furthermore, the analysis method further includes:

[0096] Determine the optimal area difference threshold for determining whether vertebral bodies at different positions have undergone overall compression from the preset range;

[0097] The step of determining whether the vertebral body in the target vertebral body image has undergone overall compression based on the area difference includes:

[0098] In response to the area difference being greater than the optimal area difference threshold, determine that the current target vertebral body has undergone overall compression;

[0099] In response to the area difference being less than or equal to the optimal area difference threshold, determine that the area target vertebral body has not undergone overall compression.

[0100] Specifically, in this solution, by pre - constructing a data set, the data structure S1 is {T4, area1, flag1}, where flag is 0 or 1. 0 represents no overall compression, and 1 represents overall compression; the data structure S2 is {{T5, T4}, area difference, type}, where the area difference is the theoretical area of two adjacent vertebrae T5 in the current image - the theoretical area of T4, and type is no overall compression and overall compression. According to the position information of {T5, T4} in S2, all data are divided into two groups by type category. One group corresponds to overall compression, and the other group corresponds to no overall compression. The data and distribution of the area differences of normal vertebrae, as well as the data and distribution of the area differences of compressed vertebrae, are respectively counted.

[0101] Based on the above - mentioned pre - constructed data set, it is determined that the range of the area difference between normal vertebrae at the position of each vertebra (such as T4, T12, L1, L4) is [a1, b1], and the range of the area difference between the vertebra with overall compression and adjacent normal vertebrae is [a2, b2]; then, the area difference s between any vertebra and its adjacent vertebra is calculated. If the area difference s falls within the range [a2, b2], it is determined that the current target vertebra has overall compression; otherwise, it is determined that the current target vertebra has no overall compression. Among them, the ranges of a1, b1, a2, and b2 are [0, 1], and a1 < b1, a2 < b2;

[0102] In addition, for the data set at each vertebra position, the corresponding optimal area - difference threshold c for identifying whether overall compression occurs at each vertebra position (such as T4, T12, L1, L4) is obtained through the grid - search algorithm, and this optimal area - difference threshold c is also used to distinguish whether the vertebra has overall compression;

[0103] Among them, the process of the grid - search algorithm for determining the corresponding area - difference threshold for identifying whether overall compression occurs at each vertebra position is as follows:

[0104] The grid - search algorithm adjusts the area difference to judge whether all vertebrae have overall compression in the sample space, and finally selects the optimal area - difference threshold c through the AUC (area under the curve) index. Specifically:

[0105] Define the search space:

[0106] Set the interpolation coefficient The search range of, that is, the number ∈(0, 1];

[0107] Determine the step size step of the grid search, such as step = 0.01;

[0108] Among them, the step size determines the fineness of the search. A smaller step size will increase the accuracy of the search, but at the same time will significantly increase the computational amount. For example, in this experiment, there are 100 results.

[0109] Implementation logic for obtaining the area difference threshold of each vertebral body:

[0110] for current area difference ∈ (0, 1);

[0111] for each current vertebral body in list (all vertebral bodies);

[0112] Calculate whether the current vertebral body undergoes overall compression through the current area difference -> label1

[0113] If (current vertebral body.label == label1)

[0114] TrueCount++;

[0115] Else

[0116] FalseCount++;

[0117] Calculate the current AUC_i (i ≤ 100 and i is a positive integer)

[0118] current area difference = current area difference + step;

[0119] Through traversal search , obtain the optimal AUC_i, which means that the current area difference corresponding to the optimal AUC_i is the optimal area difference threshold c at the current vertebral body position; and so on, to obtain the optimal area difference threshold c at each vertebral body position.

[0120] In addition, based on experience, the optimal area difference threshold c can be determined as a preset fixed value (such as 0.2) for direct judgment; among them, the formula for calculating the optimal area difference threshold c is: c = |theoretical area of the current vertebral body - theoretical area of the previous vertebral body| / theoretical area of the previous vertebral body = 0.2;

[0121] After obtaining the optimal area difference threshold c at each vertebral body position, when actually judging whether a vertebral body undergoes overall compression, it is only necessary to judge whether the theoretical area difference between the current vertebral body and the previous vertebral body is greater than the optimal area difference threshold c. If it is greater than the optimal area difference threshold c, it is determined that the vertebral body undergoes overall compression; otherwise, it is determined that the vertebral body does not undergo overall compression. The value of the optimal area difference threshold c can be determined according to experience, and of course, it can also be adaptively adjusted according to actual needs; in a preferred solution, the value of the optimal area difference threshold c is 0.2,

[0122] In this solution, the difference between the theoretical area of the current vertebral body and the theoretical area of the adjacent vertebral body is calculated, and based on the situation of this difference, it is determined whether the vertebral body has undergone overall compression or not, ensuring the accuracy, efficiency, and reliability of the determination result.

[0123] In an implementable solution, this analysis method further includes:

[0124] In response to the current target vertebral body not having undergone overall compression, obtain the target true area corresponding to the current target vertebral body;

[0125] Calculate the first ratio of the target true area and the target theoretical area corresponding to each target vertebral body;

[0126] Determine the fracture grade of the current target vertebral body based on the first ratio;

[0127] Among them, different fracture grades correspond to different preset ratio ranges.

[0128] The first ratio = the target true area of the vertebral body / the target theoretical area of the vertebral body; different fracture grades correspond to different preset ratio ranges: for r = 0 degree: [0, 10]; for r = 1 degree: (10, 20]; for r = 2 degree: (20, 40]; for r = 3 degree: (40, 100];

[0129] In this solution, the vertebral compression fracture analysis model is used to predict whether overall compression has occurred in the vertebral body. Once it is determined that the vertebral body has not undergone overall compression, the default scheme for determining the fracture grade of the vertebral body is adopted for calculation and processing to ensure the efficiency and accuracy of timely determination of the fracture grade of the target vertebral body.

[0130] In an implementable solution, as Figure 3 shown, this analysis method further includes:

[0131] S106. In response to the current target vertebral body having undergone overall compression, adopt a preset correction scheme to determine the fracture grade of the target vertebral body.

[0132] In this solution, the vertebral compression fracture analysis model is used to predict whether overall compression has occurred in the vertebral body. Once it is determined that a certain vertebral body has undergone overall compression, the default scheme for determining the fracture grade of the vertebral body is no longer adopted, and it is necessary to promptly call the correction scheme to perform timely deviation correction on this vertebral body to ensure the accuracy of timely determination of the fracture grade of the target vertebral body.

[0133] In an implementable solution, step S106 includes:

[0134] In response to the current target vertebral body having undergone overall compression, obtain the first theoretical area of the vertebral body before the current target vertebral body and the second theoretical area of the vertebral body after it;

[0135] Use the weighted average of the first theoretical area and the second theoretical area as the new theoretical area of the current target vertebral body;

[0136] Obtain the target true area of the current target vertebral body;

[0137] Calculate the second ratio of the target true area and the new theoretical area of the current target vertebral body;

[0138] Determine the fracture classification of the current target vertebral body based on the second ratio;

[0139] Among them, different fracture classifications correspond to different preset ratio ranges;

[0140] Second ratio = target true area of the vertebral body / new theoretical area of the vertebral body; different fracture classifications correspond to different preset ratio ranges: r is 0 degree: [0, 10]; r is 1 degree: (10, 20]; r is 2 degree: (20, 40]; r is 3 degree: (40, 100];

[0141] In this solution, the value obtained by weighted averaging the theoretical areas of the adjacent vertebral bodies before and after is used as the theoretical area of the current target vertebral body, rather than directly calculating the theoretical area of the current target vertebral body for subsequent fracture classification calculations, so as to ensure the accuracy of determining the fracture classification of the vertebral body with overall compression.

[0142] In an implementable solution, step S106 includes:

[0143] In response to overall compression of the current target vertebral body, obtain the target true area and the target theoretical area of the current target vertebral body;

[0144] Calculate the third ratio of the target true area and the new theoretical area of the current target vertebral body, and multiply the third ratio by a preset coefficient to obtain a fourth ratio;

[0145] Determine the fracture classification of the current target vertebral body based on the fourth ratio;

[0146] Among them, different fracture classifications correspond to different preset ratio ranges.

[0147] Fourth ratio = true area of the vertebral body / theoretical area of the vertebral body; different fracture classifications correspond to different preset ratio ranges: r is 0 degree: [0, 10]; r is 1 degree: (10, 20]; r is 2 degree: (20, 40]; r is 3 degree: (40, 100];

[0148] Specifically, the fourth ratio = the true area of the vertebral body / (a * the theoretical area of the vertebral body), where a is a preset coefficient; the range of a is [1.19 - 1.3]. The fracture grading after correction can highly coincide with the gold standard grading result x of the vertebra with overall compression marked by experts.

[0149] In this solution, the theoretical area of the vertebral body with overall compression is directly corrected by a preset system, rather than using the theoretical area directly calculated for the current target vertebral body for subsequent fracture grading calculation, which ensures the accuracy of determining the fracture grading of the vertebral body with overall compression, simplifies the intermediate data processing process, and further improves the overall data processing efficiency.

[0150] In an implementable solution, as Figure 4 shown, the steps of calculating the target true area of the target vertebral body include:

[0151] S201. Adopt a preset edge processing rule to identify each edge in the target vertebral body to obtain the true edge line corresponding to each edge;

[0152] Among them, obtain the edge association information of each edge of the vertebral body in the vertebral body image;

[0153] Among them, the edges of the vertebral body in the vertebral body image include the front edge, the rear edge, the upper edge, and the lower edge; the edge association information includes all information that can characterize and explain the corresponding edges, including but not limited to whether there are multiple curves and the curve position information.

[0154] Judge whether the edge association information meets the corresponding preset conditions. If it meets, execute the processing of the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge;

[0155] Among them, different edges correspond to different preset conditions and preset edge processing rules; that is, for the front edge, the rear edge, the upper edge, and the lower edge, due to their different presented states, they each have their own characteristics, and the matching preset conditions are pre - constructed to automatically identify the actual situation of the corresponding edge, and then the matching preset edge processing rule is used for targeted processing to extract the true edge line corresponding to each edge.

[0156] S202. Generate the true edge graph corresponding to the vertebral body based on each true edge line;

[0157] According to the four true edge lines of the front edge, the rear edge, the upper edge, and the lower edge, automatically form a closed quadrilateral to obtain the true edge graph corresponding to the vertebral body. As Figure 5 shown, it is the vertebral body image without processing; as Figure 6As shown, where the quadrilateral (black frame) is the true edge graph of the vertebral body.

[0158] S203. Obtain the area corresponding to the true edge graph and use it as the target true area of the target vertebral body;

[0159] Among them, different edges of the vertebral body are processed using corresponding preset edge processing rules to automatically, quickly and accurately extract each actual edge line of each vertebral body, ensuring the processing accuracy and efficiency of segmenting the actual shape of the vertebral body from the vertebral body image, and further ensuring the accuracy and efficiency of determining the true area of the vertebral body.

[0160] Furthermore, for the posterior edge of the vertebral body, the edge association information is the posterior edge information; among them, the posterior edge is the edge at the connection with the pedicle in the vertebral body.

[0161] The steps of judging whether the edge association information meets the corresponding preset conditions include:

[0162] Judge whether there is a ghosting area in the vertebral body image according to the posterior edge information. If so, execute the processing of the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge;

[0163] Among them, whether there is a ghosting area is distinguished by judging whether there are multiple initial posterior edge curves. That is, once there are ≥2 initial posterior edge curves, there will be a ghosting area, that is, the actual position of the posterior edge cannot be directly distinguished.

[0164] The steps of processing the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge include:

[0165] Obtain multiple initial posterior edge curves;

[0166] Identify the pedicle area in the vertebral body image;

[0167] Use image recognition technology to analyze the vertebral body image to obtain the pedicle area and the area where the vertebral body is located.

[0168] Obtain the initial posterior edge curve that is farthest from the pedicle area as the target posterior edge curve;

[0169] Specifically, each irregular initial posterior edge curve can be unified into a line segment shape, as Figure 7 shown, calculate the distances between different line segments and the pedicle area to obtain the initial posterior edge curve that is farthest from the pedicle area.

[0170] Obtain the first upper endpoint and the first lower endpoint corresponding to the target posterior edge curve;

[0171] Among them, the intersection point of the target trailing edge curve and the curve closest to the inner side at the upper edge of the vertebral body is the first upper endpoint; the intersection point of the target trailing edge curve and the curve closest to the inner side at the lower edge of the vertebral body is the first lower endpoint; that is, the target trailing edge curve intersects with other curves only at the endpoints, and there is no intersection with other curves at other places.

[0172] Form a first line segment by connecting the first upper endpoint and the first lower endpoint.

[0173] Based on the first line segment, obtain the true trailing edge line of the vertebral body.

[0174] Furthermore, for the upper edge of the vertebral body, the edge association information is the upper edge information.

[0175] The steps of determining whether the edge association information meets the corresponding preset conditions include:

[0176] Judge whether there are multiple initial upper edge curves according to the upper edge information. If so, execute the step of processing the edge association information with the matching preset edge processing rule to generate the true edge line corresponding to each edge.

[0177] The steps of processing the edge association information with the matching preset edge processing rule to generate the true edge line corresponding to each edge include:

[0178] Obtain a first target curve within the range of the first preset position between the two initial upper edge curves located on the outermost and innermost sides, and obtain the true upper edge line based on the obtained first target curve.

[0179] Among them, the first target curve can be obtained by direct selection or automatic generation, etc.

[0180] Furthermore, the steps of determining whether the edge association information meets the corresponding preset conditions include:

[0181] Judge whether there are multiple initial upper edge curves according to the upper edge information. If not, execute the step of processing the edge association information with the matching preset edge processing rule to generate the true edge line corresponding to each edge.

[0182] The steps of processing the edge association information with the matching preset edge processing rule to generate the true edge line corresponding to each edge include:

[0183] If there is only one initial upper edge curve, obtain the true upper edge line based on the initial upper edge curve.

[0184] Furthermore, for the lower edge of the vertebral body, the edge association information is the lower edge information.

[0185] The steps of determining whether the edge association information meets the corresponding preset conditions include:

[0186] Judge whether there are multiple initial lower edge curves according to the lower edge information. If so, perform processing on the edge association information using the matching preset edge processing rules to generate the true edge line corresponding to each edge;

[0187] The steps of performing processing on the edge association information using the matching preset edge processing rules to generate the true edge line corresponding to each edge include:

[0188] Obtain a second target curve within the range of the second preset position between the two initial lower edge curves located at the innermost and outermost sides, and obtain the true lower edge line based on the obtained second target curve;

[0189] Among them, the second target curve can be obtained by direct selection or automatic generation, etc.;

[0190] Furthermore, the steps of judging whether the edge association information meets the corresponding preset conditions include:

[0191] When it is determined according to the lower edge information that there are no multiple initial lower edge curves, perform processing on the edge association information using the matching preset edge processing rules to generate the true edge line corresponding to each edge;

[0192] The steps of performing processing on the edge association information using the matching preset edge processing rules to generate the true edge line corresponding to each edge include: if there is only one initial lower edge curve, obtain the true lower edge line based on the initial lower edge curve.

[0193] Furthermore, if the vertebral body image contains two initial upper edge curves, form a first target curve within the position range from one-third to two-thirds between the two initial upper edge curves, and obtain the true upper edge line based on the formed first target curve;

[0194] Among them, the trend of the first target curve is between the trends of the two initial upper edge curves.

[0195] As Figure 8 shown, judge whether there is a bilateral edge at the upper edge. When there are 2 initial upper edge curves (curves 1 and 2 marked in the figure), arbitrarily generate a curve between one-third and two-thirds between the 2 initial lower edge curves as the first target curve, and use this first target curve as the basis to finally determine the true upper edge line of the vertebral body, avoiding determining the true upper edge line based on the curve located more outside or the curve located more inside, so as to avoid the situation where the recognition accuracy of the true upper edge line cannot be guaranteed, and effectively ensure the accuracy of obtaining the true upper edge line.

[0196] Further, when two initial lower edge curves are included in the vertebral body image, a second target curve is formed within the range of the position area from one-third to two-thirds between the two initial lower edge curves, and the true lower edge line is obtained based on the formed second target curve;

[0197] Among them, the trend of the second target curve is between the trends of the two initial lower edge curves.

[0198] As Figure 8 shown, it is judged whether there is a bilateral edge at the lower edge. When there are 2 initial lower edge curves (curves 3 and 4 marked in the figure), an arbitrary curve is generated between 1 / 3 and 2 / 3 between the 2 initial lower edge curves as the second target curve, and the true lower edge line of the vertebral body is finally determined based on this second target curve, avoiding determining the true lower edge line based on the curve located more outside or more inside, so as to prevent the situation where the recognition accuracy of the true lower edge line cannot be guaranteed, and effectively ensuring the accuracy of obtaining the true lower edge line.

[0199] Further, the first target curve is the midline between the two initial upper edge curves;

[0200] As Figure 9 shown, the midline between the 2 initial upper edge curves (the curve in the middle of curves 1 and 2 marked in the figure) is used as the first target curve, and the overall undulating trend of the first target curve is between the 2 initial upper edge curves.

[0201] Further, the second target curve is the midline between the two initial lower edge curves.

[0202] As Figure 9 shown, the midline between the 2 initial lower edge curves (the curve in the middle of curves 3 and 4 marked in the figure) is used as the second target curve, and the overall undulating trend of the second target curve is between the 2 initial lower edge curves.

[0203] Further, the first intersection point of the straight line where the first line segment is located and the first target curve, and the second intersection point with the second target curve are obtained;

[0204] The second line segment between the first intersection point and the second intersection point is used as the true posterior edge line of the vertebral body;

[0205] Among them, the first intersection point is the first left end point of the true upper edge line, and the second intersection point is the second left end point of the true lower edge line.

[0206] Among them, based on the first line segment, the curve where the true upper edge line is located, and the curve where the true lower edge line is located, the intersection points existing between the straight line where the first line segment is located and the first target curve and the second target curve area are automatically obtained to finally determine the two end points of the true posterior edge line, thereby ensuring the rationality and accuracy of obtaining the true posterior edge line.

[0207] Further, for the anterior edge of the vertebral body, the edge association information is the anterior edge information;

[0208] The steps of determining whether the edge association information meets the corresponding preset conditions include:

[0209] Determine whether there is a bend on the initial anterior edge curve with an inward concavity greater than the first set threshold. If so, perform processing on the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge;

[0210] The steps of performing processing on the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge include:

[0211] Determine that there is osteophyte formation on the anterior edge of the vertebral body;

[0212] For example, if there is a bend on the curve with an inward concavity greater than 15°, it is determined that there is osteophyte formation at this bend.

[0213] Based on the vertex of the bend, perform tangent processing using the preset slope to obtain the initial straight line;

[0214] Among them, as Figure 10 shown in the initial straight line, the preset slope is determined based on the first slope of the first line segment;

[0215] The initial straight line obtained thereby is used to remove the part with osteophyte formation at the anterior edge of the vertebral body to avoid the influence of osteophyte formation on the determination of the edge of the vertebral body.

[0216] Obtain the third intersection point of the initial straight line and the first target curve, and the fourth intersection point with the second target curve;

[0217] Take the third line segment between the third intersection point and the fourth intersection point as the true anterior edge line of the vertebral body;

[0218] Among them, the third intersection point is the first right endpoint of the true upper edge line, and the fourth intersection point is the second right endpoint of the true lower edge line.

[0219] Further, the difference between the preset slope and the first slope is less than the second set threshold; where the second set threshold is a non-zero number;

[0220] Or, the preset slope is equal to the first slope.

[0221] Further, the steps of determining whether the edge association information meets the corresponding preset conditions include:

[0222] When there is no bend on the initial leading edge curve with a concave curvature greater than the first set threshold, the edge association information is processed using the matching preset edge processing rule to generate the true edge line corresponding to each edge;

[0223] The step of processing the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge includes:

[0224] Determine that there is no bone hyperplasia on the leading edge of the vertebral body, and use the initial leading edge curve as the target leading edge curve;

[0225] Obtain the fifth intersection point of the target leading edge curve and the true upper edge line, and the sixth intersection point with the true lower edge line;

[0226] Use the connecting curve between the fifth intersection point and the sixth intersection point as the true leading edge line of the vertebral body;

[0227] Among them, the fifth intersection point is the first right endpoint of the true upper edge line, and the sixth intersection point is the second right endpoint of the true lower edge line;

[0228] Furthermore, the step of processing the edge association information using the matching preset edge processing rule to generate the true edge line corresponding to each edge includes:

[0229] Obtain the second upper endpoint and the second lower endpoint corresponding to the target leading edge curve;

[0230] Use the fourth line segment formed by connecting the second upper endpoint and the second lower endpoint;

[0231] Obtain the seventh intersection point of the straight line where the fourth line segment is located and the true upper edge line, and the eighth intersection point with the true lower edge line;

[0232] Use the fifth line segment between the seventh intersection point and the eighth intersection point as the true leading edge line of the vertebral body;

[0233] Among them, the seventh intersection point is the first right endpoint of the true upper edge line, and the eighth intersection point is the second right endpoint of the true lower edge line.

[0234] The step of obtaining the true upper edge line based on the acquired first target curve includes:

[0235] Use the curve between the first left endpoint and the first right endpoint as the true upper edge line;

[0236] Furthermore, the step of obtaining the true lower edge line based on the acquired second target curve includes:

[0237] Use the curve between the second left endpoint and the second right endpoint as the true lower edge line.

[0238] Further, based on the true trailing edge line, true leading edge line, true upper edge line, and true lower edge line, a closed quadrilateral is formed, and the quadrilateral is used as the true edge figure of the vertebral body.

[0239] In an implementable solution, as Figure 11 shown, the steps for calculating the target theoretical area of the target vertebral body include:

[0240] S301. Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain the true edge line corresponding to each edge;

[0241] Among them, the true edge line includes a true leading edge line, a true trailing edge line, a true upper edge line, and a true lower edge line;

[0242] S302. Based on each true edge line, form a closed figure and use it as the true edge figure corresponding to the vertebral body;

[0243] S303. Generate a circumscribed rectangle of the true edge figure according to the true leading edge line, true trailing edge line, and true edge figure, and use the circumscribed rectangle as the theoretical shape of the target vertebral body;

[0244] S304. Obtain the area corresponding to the theoretical shape and use it as the target theoretical area of the target vertebral body.

[0245] Among them, different edges of the vertebral body are processed using corresponding preset edge processing rules to automatically, quickly, and accurately extract each actual edge line of each vertebral body, ensuring the processing accuracy and efficiency of segmenting the actual shape of the vertebral body from the vertebral body image; furthermore, a circumscribed rectangle is formed with the true leading edge line, true trailing edge line, and the quadrilateral closed by the four true edge lines as the corresponding theoretical shape of the vertebral body in the vertebral body image, thereby ensuring the accuracy and rationality of determining the theoretical shape of the vertebral body, and further ensuring the accuracy of determining the theoretical area.

[0246] Further, step S303 includes:

[0247] Select the longer line segment among the true leading edge line and the true trailing edge line as the width side length of the circumscribed rectangle;

[0248] Among them, the true leading edge line and the true trailing edge line are parallel to each other;

[0249] Use the distance between the true leading edge line and the true trailing edge line as the length side length of the circumscribed rectangle;

[0250] Generate a circumscribed rectangle of the true edge figure according to the width side length and the length side length, and use the circumscribed rectangle as the theoretical shape of the vertebral body.

[0251] The specific implementation processes in the other steps S301 and S302 are the same as those for calculating the target true area of the vertebral body above, and thus will not be elaborated herein.

[0252] Embodiment 3

[0253] As Figure 12 shown, the analysis system for the state of vertebral body compression fractures in this embodiment includes:

[0254] A target image acquisition module 1 for acquiring a target vertebral body image;

[0255] A target vertebral body segmentation module 2 for segmenting and extracting each target vertebral body in the target vertebral body image;

[0256] A vertebral body information acquisition module 3 for acquiring the target position and target theoretical area corresponding to each target vertebral body;

[0257] An area difference calculation module 4 for calculating the area difference between the theoretical areas of the target vertebral body and the adjacent vertebral bodies at the current target position;

[0258] A state analysis module 5 for determining whether the current target vertebral body has undergone overall compression based on the area difference.

[0259] In the present disclosure, based on the segmentation and extraction of each vertebral body in the vertebral body image, considering the area correlation data of the vertebral bodies, calculating the area difference situation between adjacent vertebral bodies to determine whether the current target vertebral body has undergone overall compression effectively ensures the accuracy, efficiency, and reliability of the determination result of whether the vertebral body in any target vertebral body image has undergone overall compression output by the vertebral body compression fracture analysis model, and has the advantages of convenient analysis process and small calculation amount.

[0260] For the system embodiment, since it basically corresponds to the method embodiment, the relevant parts can refer to the partial description of the method embodiment. The system embodiments described above are merely illustrative, where the units described as separate components may or may not be physically separated, and the components as units may or may not be physical units, that is, they may be located in one place or may be distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the present disclosure solution.

[0261] Embodiment 4

[0262] As Figure 13 shown, the analysis system for the state of vertebral body compression fractures in this embodiment is a further improvement on Embodiment 3. Specifically:

[0263] In an implementable solution, the state analysis module 5 is configured to:

[0264] In response to the area difference not falling within the preset range representing the area difference between normal vertebral bodies, it is determined that the current target vertebral body has undergone overall compression;

[0265] In response to the area difference falling within the preset range, it is determined that the area target vertebral body has not undergone overall compression.

[0266] In an implementable solution, the analysis system further includes:

[0267] A threshold determination module 6, configured to determine, from the preset range, an optimal area difference threshold for determining whether vertebral bodies at different positions have undergone overall compression;

[0268] The state analysis module 5 is configured to determine that the current target vertebral body has undergone overall compression in response to the area difference being greater than the optimal area difference threshold;

[0269] In response to the area difference being less than or equal to the optimal area difference threshold, it is determined that the area target vertebral body has not undergone overall compression.

[0270] In an implementable solution, the analysis system further includes a fracture grading determination module 7;

[0271] The fracture grading determination module 7 is configured to obtain the target true area corresponding to the current target vertebral body in response to the current target vertebral body not having undergone overall compression;

[0272] Calculate the first ratio of the target true area and the target theoretical area corresponding to each target vertebral body;

[0273] Determine the fracture grade of the current target vertebral body based on the first ratio;

[0274] Among them, different fracture grades correspond to different preset ratio ranges.

[0275] In an implementable solution, the fracture grading determination module 7 is configured to determine the fracture grade of the target vertebral body by using a preset correction scheme in response to the current target vertebral body having undergone overall compression.

[0276] In an implementable solution, the fracture grading determination module 7 is configured to:

[0277] In response to the current target vertebral body having undergone overall compression, obtain the first theoretical area of the vertebral body before the current target vertebral body and the second theoretical area of the vertebral body after the current target vertebral body;

[0278] Use the weighted average of the first theoretical area and the second theoretical area as the new theoretical area of the current target vertebral body;

[0279] Obtain the target true area of the current target vertebral body;

[0280] Calculate the second ratio of the target true area and the new theoretical area of the current target vertebral body;

[0281] Determine the fracture classification of the current target vertebral body based on the second ratio;

[0282] Among them, different fracture classifications correspond to different preset ratio ranges;

[0283] In an implementable solution, the fracture classification determination module 7 is used for:

[0284] In response to the overall compression of the current target vertebral body, obtain the target true area and the target theoretical area of the current target vertebral body;

[0285] Calculate the third ratio of the target true area and the new theoretical area of the current target vertebral body, and multiply the third ratio by a preset coefficient to obtain the fourth ratio;

[0286] Determine the fracture classification of the current target vertebral body based on the fourth ratio;

[0287] Among them, different fracture classifications correspond to different preset ratio ranges.

[0288] In an implementable solution, the steps of calculating the target true area of the target vertebral body include:

[0289] Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain the true edge line corresponding to each edge;

[0290] Generate the true edge graph corresponding to the vertebral body based on each true edge line;

[0291] Obtain the area corresponding to the true edge graph and use it as the target true area of the target vertebral body;

[0292] In an implementable solution, the steps of calculating the target theoretical area of the target vertebral body include:

[0293] Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain the true edge line corresponding to each edge;

[0294] Among them, the true edge line includes the true front edge line, the true rear edge line, the true upper edge line and the true lower edge line;

[0295] Based on each true edge line, form a closed graph and use it as the true edge graph corresponding to the vertebral body;

[0296] Generate the circumscribed rectangle of the true edge graph according to the true front edge line, the true rear edge line and the true edge graph, and use the circumscribed rectangle as the theoretical shape of the target vertebral body;

[0297] Obtain the area corresponding to the theoretical shape and use it as the target theoretical area of the target vertebral body.

[0298] For the system embodiment, since it basically corresponds to the method embodiment, the relevant parts can be referred to the partial description of the method embodiment. The system embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated. The components as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the present disclosure solution.

[0299] Embodiment 5

[0300] Figure 14 FIG. is a schematic structural diagram of an electronic device shown in an exemplary embodiment of the present disclosure. The electronic device includes a memory, a processor, and a computer program stored in the memory and configured to run on the processor. When the processor executes the computer program, the method described in any of the above embodiments is implemented. Figure 14 The displayed electronic device 90 is merely an example and should not impose any limitations on the functions and usage scopes of the embodiments of the present disclosure.

[0301] As Figure 14 shown, the electronic device 90 may be presented in the form of a general computing device. For example, it may be a server device. The components of the electronic device 90 may include, but are not limited to: at least one of the above processors 91, at least one of the above memories 92, and a bus 93 connecting different system components (including the memory 92 and the processor 91).

[0302] The bus 93 includes a data bus, an address bus, and a control bus.

[0303] The memory 92 may include volatile memory, such as a random access memory (RAM) 921 and / or a cache memory 922, and may further include a read-only memory (ROM) 923.

[0304] The memory 92 may further include a program tool 925 (or utility) having a set of (at least one) program modules 924. Such program modules 924 include, but are not limited to: an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include the implementation of a network environment.

[0305] The processor 91 executes various functional applications and data processing by running the computer program stored in the memory 92, such as the method provided in any of the above embodiments.

[0306] The electronic device 90 may also communicate with one or more external devices 94 (such as a keyboard, a pointing device, etc.). Such communication may be carried out through the input / output (I / O) interface 95. Moreover, the electronic device 90 may also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN) and / or a public network, such as the Internet) through the network adapter 96. As shown in the figure, the network adapter 96 communicates with other modules of the electronic device 90 through the bus 93. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in combination with the electronic device 90, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID (redundant array of independent disks) systems, tape drives, and data backup storage systems, etc.

[0307] It should be noted that, although several units / modules or sub-units / modules of the electronic device are mentioned in the above detailed description, such division is merely exemplary and not mandatory. In fact, according to the embodiments of the present disclosure, the features and functions of two or more units / modules described above may be embodied in one unit / module. Conversely, the features and functions of one unit / module described above may be further divided and embodied by multiple units / modules.

[0308] Embodiment 6

[0309] The embodiments of the present disclosure also provide a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the method provided by any one of the above embodiments is implemented.

[0310] Among them, the more specific computer-readable storage medium that can be adopted may include but not limited to: portable disks, hard disks, random access memories, read-only memories, erasable programmable read-only memories, optical storage devices, magnetic storage devices, or any suitable combination of the above.

[0311] Embodiment 7

[0312] The embodiments of the present disclosure also provide a computer program product, including a computer program, and when the computer program is executed by a processor, the method described in any one of the above is implemented.

[0313] Among them, the program code for executing the computer program product of the present disclosure may be written in any combination of one or more programming languages, and the program code may be executed entirely on the user device, partially on the user device, executed as an independent software package, partially on the user device and partially on a remote device, or entirely on a remote device.

[0314] Although the specific embodiments of the present disclosure have been described above, those skilled in the art should understand that this is only an example, and the protection scope of the present disclosure is defined by the appended claims. Without departing from the principle and essence of the present disclosure, those skilled in the art can make various changes or modifications to these embodiments, but these changes and modifications all fall within the protection scope of the present disclosure.

Claims

1. A method for analyzing the state of vertebral compression fractures, characterized in that, The analysis method includes: Obtain a target vertebral body image; Segment and extract each target vertebral body in the target vertebral body image; Obtain the target position and target theoretical area corresponding to each target vertebral body; Calculate the area difference between the theoretical areas of the target vertebral body at the current target position and the adjacent vertebral bodies; Based on the area difference, determine whether the current target vertebral body has undergone overall compression.

2. The method for analyzing the state of vertebral compression fracture according to claim 1, wherein, The step of determining whether the vertebral body in the target vertebral body image has undergone overall compression based on the area difference includes: In response to the area difference not falling within a preset range representing the area difference between normal vertebral bodies, determine that the current target vertebral body has undergone overall compression; In response to the area difference falling within the preset range, determine that the target vertebral body has not undergone overall compression.

3. The method for analyzing the state of vertebral compression fracture according to claim 2, wherein The analysis method further includes: Determine an optimal area difference threshold from the preset range for determining whether vertebral bodies at different positions have undergone overall compression; The step of determining whether the vertebral body in the target vertebral body image has undergone overall compression based on the area difference includes: In response to the area difference being greater than the optimal area difference threshold, determine that the current target vertebral body has undergone overall compression; In response to the area difference being less than or equal to the optimal area difference threshold, determine that the target vertebral body has not undergone overall compression.

4. The method for analyzing the state of vertebral compression fracture according to claim 2 or 3, characterized in that, The analysis method further includes: In response to the current target vertebral body not having undergone overall compression, obtain the target true area corresponding to the current target vertebral body; Calculate the first ratio of the target true area and the target theoretical area corresponding to each target vertebral body; Based on the first ratio, determine the fracture classification of the current target vertebral body; Wherein, different fracture classifications correspond to different preset ratio ranges.

5. The method for analyzing the state of vertebral compression fractures according to claim 4, characterized in that, The analysis method further includes: In response to the current target vertebral body having undergone overall compression, use a preset correction scheme to determine the fracture classification of the target vertebral body.

6. The method for analyzing the state of vertebral compression fractures according to claim 5, wherein The step of, in response to the current target vertebral body having undergone overall compression, using a preset correction scheme to determine the fracture classification of the target vertebral body includes: In response to the current target vertebral body having undergone overall compression, obtain the first theoretical area of the previous vertebral body and the second theoretical area of the next vertebral body of the current target vertebral body; Use the weighted average of the first theoretical area and the second theoretical area as the new theoretical area of the current target vertebral body; Obtain the target true area of the current target vertebral body; Calculate the second ratio of the target true area and the new theoretical area of the current target vertebral body; Based on the second ratio, determine the fracture classification of the current target vertebral body; Wherein, different fracture classifications correspond to different preset ratio ranges; Or, The step of, in response to the current target vertebral body having undergone overall compression, using a preset correction scheme to determine the fracture classification of the target vertebral body includes: In response to the current target vertebral body having undergone overall compression, obtain the target true area and the target theoretical area of the current target vertebral body; Calculate the third ratio of the target true area and the new theoretical area of the current target vertebral body, and multiply the third ratio by a preset coefficient to obtain a fourth ratio; Determine the fracture classification of the current target vertebral body based on the fourth ratio; Among them, different fracture classifications correspond to different preset ratio ranges.

7. The method for analyzing the state of vertebral body compression fracture according to claim 6, wherein The steps of calculating the target true area of the target vertebral body include: Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain a true edge line corresponding to each edge; Generate a true edge graph corresponding to the vertebral body based on each true edge line; Obtain the area corresponding to the true edge graph and use it as the target true area of the target vertebral body; And / or The steps of calculating the target theoretical area of the target vertebral body include: Adopt a preset edge processing rule to identify and process each edge in the target vertebral body to obtain a true edge line corresponding to each edge; Among them, the true edge line includes a true front edge line, a true rear edge line, a true upper edge line, and a true lower edge line; Form a closed graph based on each true edge line and use it as the true edge graph corresponding to the vertebral body; Generate a circumscribed rectangle of the true edge graph according to the true front edge line, the true rear edge line, and the true edge graph, and use the circumscribed rectangle as the theoretical shape of the target vertebral body; Obtain the area corresponding to the theoretical shape and use it as the target theoretical area of the target vertebral body.

8. An analysis system for vertebral compression fracture status, characterized in that, The analysis system includes: A target image acquisition module for acquiring a target vertebral body image; A target vertebral body segmentation module for segmenting and extracting each target vertebral body in the target vertebral body image; A vertebral body information acquisition module for acquiring the target position and the target theoretical area corresponding to each target vertebral body; An area difference calculation module for calculating the area difference between the target vertebral body at the current target position and the theoretical area of the adjacent vertebral body; A state analysis module for determining whether the current target vertebral body has undergone overall compression based on the area difference.

9. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and configured to run on the processor, characterized in that When the processor executes the computer program, it implements the method for analyzing the state of vertebral body compression fractures as described in any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for analyzing the state of vertebral body compression fractures as described in any one of claims 1 to 7.

11. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for analyzing the state of vertebral body compression fractures as described in any one of claims 1 to 7.

Images