A registration system and method for femoral original model

Abstract

The present invention relates to a registration system and method for a femoral original model. The method comprises at least the following steps: reconstructing two-dimensional tomographic image data of a femoral head portion to establish a three-dimensional mesh model of the femoral head, wherein the three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system, and the three-dimensional mesh model at least includes an original necrosis model A and an original bone model B. The method further comprises: performing at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm based on the acetabulum lunate surface or the original necrosis model A, thereby sequentially completing a spatial coarse registration process and a spatial fine registration process, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

Description

[0001] The original basis of this divisional application is the patent application with application number CN201910878678.2, application date September 17, 2019, and invention name “A femoral integrity analysis system and a femoral integrity analysis model construction method”. Technical Field

[0002] The present invention relates to the field of medical imaging technology, and in particular to a registration system and method for a femur original model. Background Art

[0003] It is an indisputable fact that femoral head necrosis seriously endangers the health and quality of life of patients. Its treatment has also become a global problem. The difficulties are nothing more than: patients with femoral head necrosis are relatively young and it is not easy to undergo joint replacement. If joint replacement is performed, they will face the risk of multiple revisions; once the femoral head collapses, no matter what treatment is used, the original shape of the femoral head cannot be restored, and secondary osteoarthritis is likely to occur; due to the lack of precise preoperative design, the current success rate of hip preservation treatment is not completely satisfactory. If hip preservation is successful in younger patients, joint replacement can be avoided or significantly delayed, greatly reducing the risk of multiple revisions; if the femoral head collapse can be prevented, irreversible damage caused by changes in the femoral head morphology and the risk of secondary osteoarthritis can be avoided. Therefore, preventing and treating femoral head collapse and maintaining the original shape of the hip joint are important principles for maintaining femoral head health.

[0004] Treatment options for osteonecrosis of the femoral head are based on clinical data. However, conventional MRI can only estimate the area of ​​necrosis from slices taken at different planes, failing to determine the volume of necrosis. Furthermore, the scanned planes before and after treatment are inconsistent, making it difficult to accurately assess treatment efficacy. An increasing number of researchers are using finite element analysis software to verify the feasibility of surgical options. With the increasing popularity of minimally invasive surgery, more people are choosing novel treatment options. For example, core decompression tantalum rods or titanium alloy stents are implanted into the femoral head via decompression channels within the femoral head for minimally invasive surgery. However, the support area provided by these devices is relatively small, resulting in certain drawbacks. The memory metal nickel-titanium alloy tennis ball designed by Wang Yan et al. is made of woven wire. While it provides a relatively large support area, its support force is correspondingly low. Furthermore, this metal ball requires opening the femoral head, which is relatively invasive. Therefore, while there are numerous treatment options for avascular necrosis of the femoral head, no single, most reliable one exists. Finite element analysis has revealed numerous issues.

[0005] The finite element analysis method is to obtain real data through clinical experiments, use a computer to construct a three-dimensional solid model of the femoral head, and use finite element analysis software for calculation and analysis. With the help of the finite element model, in vivo tissues can be non-invasively examined, assisting in the formulation of surgical diagnosis and treatment plans and the simulation of quantitative surgery. Compared with the experimental model, it greatly reduces the influence of multiple factors, making quantitative experiments more illustrative. In order to simulate the repair effect of the support on patients with ischemic femoral head necrosis more realistically and more accurately, and to complete research that is difficult to complete with clinical experiments, it is necessary to establish a finite element model obtained by applying finite element analysis to the ischemic necrotic femoral head.

[0006] A Chinese patent (publication number CN107802378A) discloses a partial talar prosthesis with a porous structure and its design and manufacturing method. The method includes the following steps: extracting a model of the talar region based on patient CT data, reconstructing the 3D model, and exporting it as an STL file; importing the STL file into the 3D modeling software Rhinoceros for reverse modeling and design; extracting and fitting the locally necrotic surface to obtain a surface that matches the curvature of the surrounding healthy surfaces; using the fitted surface to design the prosthesis, determine the position of the fixation pins, and complete the modeling of the porous structure using the Grasshopper plug-in; fabricating a personalized partial talar prosthesis using laser selective melting technology; polishing the upper surface of the prosthesis and acid-etching the lower porous structure. This patent enables rapid prosthetic repairs, shortening patient waiting times and improving the responsiveness of prosthetic supply.

[0007] A Chinese patent (publication number CN104462636B) discloses a method for modeling a necrotic femoral head repair model based on an umbrella-shaped femoral head support. The method comprises the following steps: 1. Acquiring a 3D model of the femoral head to be repaired: Obtaining a NURBS surface model of the femoral head to be repaired, wherein the femoral head to be repaired is a femoral head with a necrotic area of ​​femoral tissue and previously repaired with a femoral head support; the femoral head support comprises an umbrella-shaped support and a support sleeve; 2. Determining the necrotic area to be separated based on the shape of the umbrella-shaped support, and establishing a necrotic femoral head model; 3. Establishing the femoral head support model; 4. Establishing a necrotic femoral head implant model: Establishing a necrotic femoral head implant model with an implantation channel and a 3D model of the implanted bone; 5. Establishing a necrotic femoral head repair model. This patented method features simple steps, a rational design, easy implementation, and excellent performance. It can easily and quickly establish a high-quality repair model of a necrotic femoral head implanted with a femoral head support.

[0008] The above patents in the prior art can all establish a necrotic femoral head model through CT / MR scanning image data of femoral head necrosis, or obtain the corresponding repair surface of the lesion site, but the modeling method provided by it can only establish a corresponding repair model of the necrotic femoral head with a femoral head support device implanted. However, for the subsequent surgical treatment of femoral head necrosis, the modeling method provided by this patent cannot accurately segment the femoral head model, and the spatial three-dimensional relationship of femoral head necrosis is not judged on the femoral head model obtained after segmentation during the model construction process. The model extraction accuracy is low, which will have a great impact on the subsequent surgical treatment of femoral head necrosis. It not only cannot meet the femoral head necrosis surgical simulation training of orthopedic surgeons, but also cannot fully adapt to the clinical needs of femoral head necrosis surgical treatment. In addition, since whether femoral head necrosis will progress to collapse is closely related to its necrotic surface area, it cannot further achieve early prediction and prevention of collapse of femoral head avascular necrosis.

[0009] In addition, on the one hand, the understanding of the technical personnel in this field by the applicant is bound to be different from that of the examination department; on the other hand, the inventor studied a large number of documents and patents when making this invention, but due to space limitations, not all details and contents are listed in detail. However, this does not mean that the present invention does not have the characteristics of these prior arts. On the contrary, the present invention already has all the characteristics of the prior art, and the applicant reserves the right to add relevant prior art to the background technology at any time in accordance with the relevant provisions of the examination guidelines. Summary of the Invention

[0010] In view of the shortcomings of the existing technology, the present invention provides a femoral integrity analysis system, which includes at least a three-dimensional model creation module, a three-dimensional model extraction module and a data analysis module. The three-dimensional model creation module is configured to: reconstruct the two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system. The three-dimensional mesh model includes at least an original necrosis model A and an original bone model B. The three-dimensional model extraction module is used to use a distance threshold about the acetabulum lunate surface or about the original necrosis model A. The algorithm performs at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, respectively, thereby completing the spatial coarse alignment process and the spatial fine alignment process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area. The data analysis module is used to cut the femoral head weight-bearing precise area and the necrosis weight-bearing precise area according to the cutting surface it establishes, so as to obtain the corresponding weight-bearing stress concentration areas on the femoral head weight-bearing precise area and the necrosis weight-bearing precise area, and calculate the completeness rate of at least one of the weight-bearing stress concentration areas.

[0011] The femoral integrity analysis system includes at least a three-dimensional model creation module. The three-dimensional model creation module is configured to reconstruct the two-dimensional tomographic image data of the femoral head to establish a three-dimensional grid model of the femoral head. The three-dimensional grid model of the femoral head is composed of a series of triangular facet grid sets in a three-dimensional coordinate system. The three-dimensional grid model includes at least the original necrosis model A and the original bone model B. The femoral integrity analysis system includes at least a three-dimensional model extraction module. The three-dimensional model extraction module is used to perform at least one distance threshold iterative calculation on the triangular facet grid sets of the original necrosis model A and the original bone model B, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence. The three-dimensional model extraction module is used to extract the obtained femoral head weight-bearing precise area and the necrosis weight-bearing precise area after completing the spatial coarse registration process and the spatial fine registration process in sequence. The distance threshold algorithm about the acetabulum lunate surface or about the original necrosis model A is used to perform at least one distance threshold iterative calculation on the triangular facet grid sets of the original necrosis model A and the original bone model B. The three-dimensional model extraction module is used to perform at least one distance threshold iterative calculation on the triangular facet mesh set of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area. The three-dimensional model extraction module is used to perform at least one distance threshold iterative calculation on the triangular facet mesh set of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm about the acetabulum lunate surface or about the original necrosis model A, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area. The femoral integrity analysis system includes at least a data analysis module. The data analysis module is used to cut the femoral head weight-bearing precise area and the necrosis weight-bearing precise area respectively according to the cutting surface it establishes, so as to obtain the corresponding load-bearing stress concentration area on the femoral head weight-bearing precise area and the necrosis weight-bearing precise area, and calculate the integrity rate of at least one of the load-bearing stress concentration areas. The data analysis module is configured to cut the femoral head weight-bearing precise region and the necrosis weight-bearing precise region according to the cutting plane it establishes, thereby obtaining corresponding weight-bearing stress concentration regions on the femoral head weight-bearing precise region and the necrosis weight-bearing precise region. The data analysis module is configured to calculate the integrity rate of at least one of the weight-bearing stress concentration regions.

[0012] According to a preferred embodiment, the three-dimensional model extraction module uses a distance threshold algorithm about the lunate surface of the acetabulum to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, respectively, to complete the spatial coarse alignment process and extract the femoral head weight-bearing area and the necrotic weight-bearing area; the three-dimensional model extraction module uses a distance threshold algorithm about the original necrosis model A to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the femoral head weight-bearing area and the necrotic weight-bearing area, respectively, to complete the spatial fine alignment process and extract the further optimized femoral head weight-bearing precise area and the necrotic weight-bearing precise area.

[0013] The three-dimensional model extraction module uses the distance threshold algorithm about the lunate surface of the acetabulum to complete the spatial coarse registration process. The three-dimensional model extraction module uses the distance threshold algorithm about the lunate surface of the acetabulum to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, thereby completing the spatial coarse registration process and extracting the femoral head weight-bearing area and the necrosis weight-bearing area. The three-dimensional model extraction module uses the distance threshold algorithm about the original necrosis model A to complete the spatial fine registration process. The three-dimensional model extraction module uses the distance threshold algorithm about the original necrosis model A to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the femoral head weight-bearing area and the necrosis weight-bearing area, thereby completing the spatial fine registration process and extracting the further optimized femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

[0014] According to a preferred embodiment, the three-dimensional model extraction module is further configured to perform necrosis segmentation preprocessing on the original necrosis model A before executing the spatial coarse registration process and the spatial fine registration process, and use the triangular facet mesh set corresponding to the necrotic surface model A1 obtained after the preprocessing as the data source of the spatial coarse registration process, wherein the necrosis segmentation preprocessing is that the three-dimensional model extraction module uses a distance threshold algorithm to traverse each triangular facet mesh of the original necrotic model A to generate a triangular facet mesh set S 11 ; The three-dimensional model extraction module extracts the triangular face mesh set S 11 The connected region with the largest area is determined to be the necrotic surface region, and a necrotic surface model A1 is established in the necrotic surface region.

[0015] The three-dimensional model extraction module is further configured to perform necrosis segmentation preprocessing on the original necrosis model A before executing the spatial coarse registration process and the spatial fine registration process. The three-dimensional model extraction module uses the triangular facet mesh set corresponding to the necrosis surface model A1 obtained after preprocessing as the data source for the spatial coarse registration process. The necrosis segmentation preprocessing at least includes the three-dimensional model extraction module using a distance threshold algorithm to traverse each triangular facet mesh of the original necrosis model A to generate a triangular facet mesh set S 11 The necrosis segmentation preprocessing at least includes the three-dimensional model extraction module extracting the triangular face mesh set S 11 The connected region with the largest area is determined to be the necrotic surface region, and a necrotic surface model A1 is established in the necrotic surface region.

[0016] According to a preferred embodiment, the spatial coarse registration process of the original necrosis model A is as follows: based on the necrosis surface model A1 obtained by performing necrosis segmentation preprocessing on the original necrosis model A by the three-dimensional model extraction module, the three-dimensional model extraction module uses a distance threshold algorithm to traverse each triangular facet mesh of the necrosis surface model A1 to generate a triangular facet mesh set S 12 ; The three-dimensional model extraction module extracts the triangular face mesh set S 12 The connected region with the largest area is determined to be the necrotic weight-bearing region, and a necrotic weight-bearing region model A2 is established in the necrotic weight-bearing region.

[0017] The spatial coarse registration process of the original necrosis model A comprises at least one or more of the following steps: performing necrosis segmentation preprocessing on the original necrosis model A based on the three-dimensional model extraction module to obtain a necrosis surface model A1; the three-dimensional model extraction module uses a distance threshold algorithm to traverse each triangular facet mesh of the necrosis surface model A1; generating a triangular facet mesh set S 12 ; The three-dimensional model extraction module extracts the triangular face mesh set S 12 determining the extracted connected region as the necrotic weight-bearing region; and establishing a necrotic weight-bearing region model A2 in the necrotic weight-bearing region.

[0018] According to a preferred embodiment, the spatial coarse registration process of the original bone model B is as follows: the three-dimensional model extraction module uses a distance threshold algorithm to traverse each triangular face mesh of the original bone model B to generate a triangular face mesh set S 21 ; The three-dimensional model extraction module extracts the triangular face mesh set S 21The connected area with the largest area is determined to be the femoral head weight-bearing area, and a femoral head weight-bearing area B1 is established in the femoral head weight-bearing area.

[0019] According to a preferred embodiment, the execution order of the spatial fine registration process of the original bone model B takes precedence over the execution order of the spatial fine registration process of the original necrosis model A, wherein the three-dimensional model extraction module uses at least a portion of the triangular facet mesh obtained after executing the spatial fine registration process on the original bone model B as a partial data source for the spatial fine registration process of the original necrosis model A. The execution order of the spatial fine registration process of the original bone model B takes precedence over the execution order of the spatial fine registration process of the original necrosis model A. The three-dimensional model extraction module uses at least a portion of the triangular facet mesh obtained after executing the spatial fine registration process on the original bone model B as a partial data source for the spatial fine registration process of the original necrosis model A.

[0020] According to a preferred embodiment, the spatial precise registration process of the original bone model B is as follows: the three-dimensional model extraction module uses a distance threshold algorithm to traverse each triangular facet mesh of the femoral head weight-bearing area B1 to generate a triangular facet mesh set S 22 ; Use the distance threshold algorithm to traverse the triangle mesh set S 22 Each triangular face mesh in generates a triangular face mesh set S 23 and the triangle mesh set S 24 ; For the triangular mesh set S 23 Perform denoising to remove isolated noise points in the set and establish a mesh with the triangle mesh set S 23 The corresponding load-bearing area model B2 and the triangle mesh set S 24 The corresponding necrotic weight-bearing area model A3, the spatial precise registration process of the original necrotic model A is: perform a Boolean union operation on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3, and obtain the final necrotic weight-bearing area model A4 after the operation. The spatial precise registration process of the original bone model B is: the three-dimensional model extraction module uses the distance threshold algorithm to traverse each triangular facet mesh of the femoral head weight-bearing area B1, and generate a triangular facet mesh set S 22 ; Use the distance threshold algorithm to traverse the triangle mesh set S 22 Each triangular face mesh in generates a triangular face mesh set S 23 and the triangle mesh set S 24 ; For the triangular mesh set S 23Perform denoising to remove isolated noise points in the set and establish a mesh with the triangle mesh set S 23 The corresponding load-bearing area model B2 and the triangle mesh set S 24 The corresponding necrotic weight-bearing area model A3. The spatial precise registration process of the original necrotic model A is: performing a Boolean union operation on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3 to obtain the final necrotic weight-bearing area model A4.

[0021] According to a preferred embodiment, the data analysis module establishes a cutting surface of the necrotic weight-bearing area model A4 according to the direction of the gravity line corresponding to the precise weight-bearing area of ​​the femoral head, and based on the fact that the two-thirds of the necrotic front and outside of the precise weight-bearing area is its weight-bearing stress concentration area, the necrotic weight-bearing area model A4 is cut along the cutting surface to extract the model A5 of the necrotic front and outside of the precise weight-bearing area of ​​the femoral head; the data analysis module establishes a cutting surface of the weight-bearing area model B2 according to the direction of the gravity line corresponding to the precise weight-bearing area of ​​the femoral head, and based on the fact that the two-thirds of the femur front and outside of the precise weight-bearing area of ​​the femoral head is its weight-bearing stress concentration area, the weight-bearing area model B2 is cut along the cutting surface to extract the model B3 of the two-thirds of the femur front and outside of the precise weight-bearing area of ​​the femoral head. The data analysis module establishes a cutting surface of the necrotic weight-bearing area model A4 according to the direction of the gravity line corresponding to the precise weight-bearing area of ​​the femoral head. The two-thirds of the necrotic front and outside of the precise weight-bearing area of ​​the necrotic weight-bearing area is its weight-bearing stress concentration area. The data analysis module cuts the necrotic weight-bearing area model A4 along the cutting surface to extract the model A5 at the anterolateral two-thirds of the necrotic weight-bearing area model A4. The data analysis module establishes the cutting surface of the weight-bearing area model B2 according to the direction of the gravity line corresponding to the precise weight-bearing area of ​​the femoral head. The anterolateral two-thirds of the femur on the precise weight-bearing area of ​​the femoral head is its weight-bearing stress concentration area. The data analysis module cuts the weight-bearing area model B2 along the cutting surface to extract the model B3 at the anterolateral two-thirds of the femur on the weight-bearing area model B2.

[0022] A method for constructing a femoral integrity analysis model comprises at least the following steps:

[0023] S1: reconstructing the two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system. The three-dimensional mesh model includes at least an original necrosis model A and an original bone model B.

[0024] S2: using a distance threshold algorithm on the acetabulum lunate surface or on the original necrosis model A, performing at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area;

[0025] S3: The precise weight-bearing area of ​​the femoral head and the precise weight-bearing area of ​​the necrosis are cut respectively according to the established cutting surface, so as to obtain the corresponding weight-bearing stress concentration areas on the precise weight-bearing area of ​​the femoral head and the precise weight-bearing area of ​​the necrosis, and calculate the integrity rate of at least one of the weight-bearing stress concentration areas.

[0026] According to a preferred embodiment, the construction method also includes the following steps: using the distance threshold algorithm about the lunate surface of the acetabulum to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, respectively, to complete the spatial coarse alignment process and extract the femoral head weight-bearing area and the necrotic weight-bearing area; using the distance threshold algorithm about the original necrosis model A to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the femoral head weight-bearing area and the necrotic weight-bearing area, respectively, to complete the spatial fine alignment process and extract the further optimized femoral head weight-bearing precise area and the necrotic weight-bearing precise area.

[0027] The present invention also relates to a registration method for an original femoral model, comprising at least the following steps: reconstructing two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head, wherein the three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system, and the three-dimensional mesh model includes at least an original necrotic model A and an original bone model B. The method further comprises traversing the original necrotic model A using a distance threshold algorithm to extract regions where the original necrotic model A does not intersect or overlap with the femoral head.

[0028] According to a preferred embodiment, the method further includes: prior to executing the coarse spatial registration process and the fine spatial registration process, the original necrosis model A is subjected to necrosis segmentation preprocessing, and the triangular facet mesh set corresponding to the necrotic surface model A1 obtained after the preprocessing is used as the data source for the coarse spatial registration process. The distance threshold algorithm based on the acetabulum lunate surface or the original necrosis model A is used to perform at least one distance threshold iteration calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, respectively, thereby completing the coarse spatial registration process and the fine spatial registration process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

[0029] According to a preferred embodiment, the method further comprises: a spatial coarse registration process of the original necrosis model A, so as to extract the weight-bearing area of ​​the necrotic surface area in the acetabulum; wherein the spatial coarse registration process of the original necrotic model A comprises: traversing the necrotic surface model A1 using a distance threshold algorithm; extracting a triangular facet mesh set S 12 The connected region with the largest area is determined to be the necrotic weight-bearing region, and a necrotic weight-bearing region model A2 is established in the necrotic weight-bearing region.

[0030] According to a preferred embodiment, the method further comprises: a spatial coarse registration process of the original bone model B, for extracting the femoral head weight-bearing area of ​​the femoral head surface area located in the acetabulum area; wherein the spatial coarse registration process of the original bone model B comprises: traversing the necrotic surface model A1 using a distance threshold algorithm; extracting a triangular facet mesh set S 21 The connected area with the largest area is determined to be the femoral head weight-bearing area, and the femoral head weight-bearing area B1 is established within the femoral head weight-bearing area. The femoral head weight-bearing area refers to the area corresponding to the surface of the femoral head projected onto the surface of the femoral head along the direction of the gravity line of the hip bone located in the acetabulum.

[0031] According to a preferred embodiment, the method further comprises: a spatial precise registration process of the original bone model B and a spatial precise registration process of the original necrosis model A; wherein the spatial precise registration process of the original bone model B is used to perform denoising on the extracted femoral head weight-bearing area and eliminate noise data; wherein the spatial precise registration process of the original bone model B comprises: using a distance threshold algorithm to traverse the femoral head weight-bearing area B1; using a distance threshold algorithm to traverse the triangular facet mesh set S22; and traversing the triangular facet mesh set S23. 23 Perform denoising to remove isolated noise points in the set and establish a mesh with the triangle mesh set S 23 The corresponding load-bearing area model B2 and the triangle mesh set S 24 The corresponding necrotic weight-bearing area model A3.

[0032] According to a preferred embodiment, the spatial precise registration process of the original necrosis model A is used to merge the extracted necrotic weight-bearing area model A3 with the necrotic weight-bearing area model A2 to achieve accurate extraction of the necrotic weight-bearing area; the spatial precise registration process of the original necrotic model A includes: performing a Boolean union operation on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3, and obtaining the final necrotic weight-bearing area model A4 after the operation.

[0033] According to a preferred embodiment, the method further includes: cutting the necrotic weight-bearing precise area according to the established cutting surface 12, thereby obtaining the corresponding weight-bearing stress concentration area on the necrotic weight-bearing precise area, and the weight-bearing stress concentration area on the necrotic weight-bearing area model A4 is the anterolateral two-thirds of the necrosis; cutting the femoral head weight-bearing precise area according to the established cutting surface 12, thereby obtaining the corresponding weight-bearing stress concentration area on the femoral head weight-bearing precise area, and the weight-bearing stress concentration area on the femoral head weight-bearing precise area is the anterolateral two-thirds of the femur. Step S302 is used to extract the femoral weight-bearing area in the anterior two-thirds of the lateral side of the acetabulum; calculating the completeness rate of at least one of the weight-bearing stress concentration areas to calculate the completeness rate of the weight-bearing area in the anterolateral two-thirds of the femur.

[0034] According to a preferred embodiment, the method further comprises: obtaining the surface area S of the necrotic anterior and lateral two-thirds model A5. A , based on the model B3 at the anterolateral two-thirds of the femur, its surface area S is obtained B , given the surface area S A and surface area S B The total surface area S of the anterolateral two-thirds weight-bearing area and the integrity rate of the lateral femoral weight-bearing area are formed by the following: , , where the surface area S is obtained based on the model B3 at the anterolateral two-thirds of the femur B The surface area S is the surface area of ​​the intact part of the femur in the anterior two-thirds of the lateral femur, and the surface area S is obtained based on the model A5 of the anterolateral two-thirds of the necrosis. A The surface area of ​​the necrotic area in the anterior two-thirds of the lateral femur is obtained based on the corresponding load-bearing stress concentration area on the necrotic load-bearing precise area. A , based on the corresponding load-bearing stress concentration area on the femoral head load-bearing area, its surface area S is obtained. B , given the surface area S A and surface area S B The total surface area S of the anterolateral two-thirds of the weight-bearing area and the integrity rate of the lateral femoral weight-bearing area, where S=S A +S B ,Rate=S B / S, where the integrity rate of the lateral femoral weight-bearing area is calculated by calculating the area of ​​the intact surface of the femur in the anterolateral two-thirds of the weight-bearing area as a percentage of the total area of ​​the weight-bearing area in that area.

[0035] The present invention also relates to a registration system for an original femoral model, comprising at least a three-dimensional model creation module, a three-dimensional model extraction module, and a data analysis module. The three-dimensional model creation module is configured to reconstruct two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system, and the three-dimensional mesh model includes at least an original necrotic model A and an original bone model B. The three-dimensional model extraction module uses a distance threshold algorithm to traverse the original necrotic model A to extract the non-intersecting and overlapping areas of the original necrotic model A and the femoral head.

[0036] According to a preferred embodiment, the three-dimensional model extraction module preferentially performs necrosis segmentation preprocessing on the original necrosis model A before executing the spatial coarse registration process and the spatial fine registration process, and uses the triangular facet mesh set corresponding to the necrotic surface model A1 obtained after the preprocessing as the data source for the spatial coarse registration process. The three-dimensional model extraction module is used to perform at least one distance threshold iterative calculation on the triangular facet mesh set of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm about the acetabulum lunate surface or about the original necrosis model A, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

[0037] The femoral integrity analysis system and femoral integrity analysis model construction method provided by the present invention have at least the following beneficial technical effects:

[0038] The femoral integrity analysis model construction method provided by the present invention uses a distance threshold algorithm to limit the positional distance between the extended line segment on the normal line of the triangular facet vertex and the acetabular lunate surface or the original necrotic model, and performs a spatial coarse registration process and a spatial fine registration process respectively. This greatly reduces the target area and can effectively ensure the convergence of the calculation results, thereby achieving accurate segmentation of the three-dimensional mesh model of the femoral head area.

[0039] In addition, the femoral integrity analysis model construction method provided by the present invention, during the model construction process, interactively analyzes and calculates the weight-bearing area model extracted by the spatial coarse alignment process and the weight-bearing area model extracted after the spatial fine alignment process segmentation, and can calculate the necrosis rate of the femoral head weight-bearing area to the maximum extent in three-dimensional space, greatly improving the accuracy of the subsequent calculated femoral head weight-bearing area necrosis rate, thereby providing more reliable data support for accurately judging the spatial stereoscopic relationship of femoral head necrosis, adapting to the clinical needs of surgical treatment of femoral head necrosis and planning of femoral head necrosis osteotomy surgery, solving the problems in the existing technology that the spatial stereoscopic relationship of femoral head necrosis cannot be accurately judged and the extracted femoral head model has low accuracy, and can further achieve early prediction and prevention of collapse of avascular necrosis of the femoral head. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 It is a simplified flowchart of the method for constructing a femoral integrity analysis model provided by the present invention;

[0041] Figure 2 1 is a simplified flow chart of a preferred method for constructing a femoral integrity analysis model provided by the present invention;

[0042] Figure 3 This is a simplified schematic diagram of the module connection relationship of the femoral integrity analysis system provided by the present invention;

[0043] Figure 4 is a simplified structural diagram of a preferred femoral integrity analysis system provided by the present invention; and

[0044] Figure 5 It is a simplified schematic diagram of the preferred cutting surface provided by the present invention.

[0045] Reference Signs List

[0046] 1: 3D model creation module 2: 3D model extraction module 3: Data analysis module

[0047] 4: Server 5: Storage device 6: Processor

[0048] 7: Bus 8: Network adapter 9: I / O interface

[0049] 10: Display 11: External device 12: Cutting surface

[0050] 13: Direction of gravity DETAILED DESCRIPTION

[0051] The present invention will be described in detail below with reference to the accompanying drawings. Example 1

[0052] like Figure 1As shown, the method for constructing a femoral integrity analysis model includes at least the following steps S1 to S3:

[0053] Step S1: Reconstruct the two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system. The three-dimensional mesh model includes at least the original necrosis model A and the original bone model B.

[0054] Preferably, adopt the three-dimensional medical imaging surface reconstruction algorithm based on segmentation to reconstruct the three-dimensional grid model to the two-dimensional tomographic image data, obtain the three-dimensional grid model of the femoral head by threshold screening, interactive segmentation and three-dimensional reconstruction process.Wherein the two-dimensional tomographic image data can be obtained by imaging equipment CT and / or imaging equipment MRI.The three-dimensional grid model of the femoral head is gridded, and the continuous geometric model is discretized to obtain the finite element model corresponding thereto.Wherein, preferably, the three-dimensional grid model of the femoral head is carried out face grid division.This three-dimensional grid model at least comprises original necrosis model A and original bone model B.

[0055] Step S2: The three-dimensional model extraction module 2 is used to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm about the acetabulum lunate surface or about the original necrosis model A, thereby completing the spatial coarse alignment process and the spatial fine alignment process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area. The femoral integrity analysis model construction method provided by the present invention uses a distance threshold algorithm to limit the positional distance between the extended line segment on the triangular facet vertex normal and the acetabulum lunate surface or the original necrosis model, and respectively performs the spatial coarse alignment process and the spatial fine alignment process, which greatly reduces the target area and can effectively ensure the convergence of the calculation results, thereby achieving accurate segmentation of the three-dimensional mesh model of the femoral head.

[0056] Among them, Figure 1 As shown, step S2 at least includes the following steps S201 to S203:

[0057] Step S201: Before executing the spatial coarse registration process and the spatial fine registration process, the original necrosis model A is preferentially preprocessed for necrosis segmentation, and the triangular facet mesh set corresponding to the necrosis surface model A1 obtained after the preprocessing is used as the data source for the spatial coarse registration process.

[0058] Preferably, step S201: using a distance threshold algorithm to traverse each triangular face mesh of the original necrotic model A to generate a triangular face mesh set S 11 ; Extract the triangular mesh set S 11The connected region with the largest area is determined as the necrotic surface region, and a necrotic surface model A1 is established in the necrotic surface region.

[0059] Further preferably, the specific steps of traversing the original necrosis model A using the distance threshold algorithm are as follows: traverse each triangular facet mesh of the original necrosis model A, extend the three vertices of each triangular facet mesh by 0.013m along the corresponding normal direction to obtain the corresponding line segments, and if the line segments corresponding to the three vertices of the triangular facet mesh do not intersect with the original bone model B, then take out the facet and put it in the set S 11 Step S201 is used to extract the non-intersecting overlapping area of ​​the original necrosis model A and the femoral head.

[0060] Step S202: Use the distance threshold algorithm about the lunate surface of the acetabulum to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, so as to complete the spatial coarse alignment process and extract the femoral head weight-bearing area and the necrotic weight-bearing area.

[0061] More specifically, step S202 includes at least a rough spatial registration process of the original necrosis model A, which is used to extract the weight-bearing area of ​​the necrotic surface area within the acetabulum region, such as Figure 2 As shown, the process is steps S2021~S2022.

[0062] Step S2021: Use the distance threshold algorithm to traverse each triangular face mesh of the necrotic surface model A1 to generate a triangular face mesh set S 12 ;

[0063] The specific steps of using the distance threshold algorithm to traverse the necrotic surface model A1 are as follows: traverse each triangular facet mesh of the original necrotic model A, extend the three vertices of each triangular facet mesh by 0.013m along the corresponding normal direction to obtain the corresponding line segments, and if the line segments corresponding to the three vertices of the triangular facet mesh intersect with the acetabulum lunate surface, then take out the facet and put it in the set S 12 middle.

[0064] Step S2022: Extracting a triangular mesh set S 12 The connected region with the largest area is determined to be the necrotic weight-bearing region, and a necrotic weight-bearing region model A2 is established in the necrotic weight-bearing region.

[0065] The necrotic weight-bearing area refers to an area of ​​the hip bone located in the acetabulum that can be projected onto the necrotic surface model A1 along the gravity line direction 13 and corresponds to the surface of the necrotic surface model A1.

[0066] More specifically, step S202 includes at least a rough spatial registration process of the original bone model B, which is used to extract the femoral head weight-bearing area of ​​the femoral head surface area located in the acetabulum area, such as Figure 2 As shown, the process is steps S2023~S2024.

[0067] Step S2023: Use the distance threshold algorithm to traverse each triangular face mesh of the original bone model B to generate a triangular face mesh set S 21 .

[0068] The specific steps of using the distance threshold algorithm to traverse the necrotic surface model A1 are as follows: traverse each triangular mesh of the original bone model B, extend the three vertices of each triangular mesh by 0.013m along the corresponding normal direction to obtain the corresponding line segments, and if the line segments corresponding to the three vertices of the triangular mesh intersect with the acetabulum lunate surface, then take out the facet and put it in the triangular mesh set S 21 middle.

[0069] Step S2024: Extracting a triangular mesh set S 21 The connected area with the largest area is determined to be the femoral head weight-bearing area, and a femoral head weight-bearing area B1 is established in the femoral head weight-bearing area.

[0070] The weight-bearing area of ​​the femoral head refers to the area on the surface of the femoral head that can be projected onto the surface of the femoral head along the direction of the gravity line 13 of the hip bone located in the acetabulum. Preferably, the weight-bearing area of ​​the femoral head refers to the area on the original bone model B that can be projected onto the original bone model B along the direction of the gravity line 13 of the hip bone located in the acetabulum.

[0071] Step S203: Use the distance threshold algorithm for the original necrosis model A to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the femoral head weight-bearing area and the necrosis weight-bearing area, respectively, so as to complete the spatial precise alignment process and extract the further optimized femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

[0072] More specifically, step S203 includes at least a spatial fine registration process for the original bone model B and a spatial fine registration process for the original necrotic model A. The spatial fine registration process for the original bone model B is executed in a sequence that precedes the sequence that precedes the spatial fine registration process for the original necrotic model A. The three-dimensional model extraction module 2 uses at least a portion of the triangular facet mesh obtained after executing the spatial fine registration process on the original bone model B as a partial data source for the spatial fine registration process of the original necrotic model A.

[0073] The femoral integrity analysis model construction method provided by the present invention, during the model construction process, interactively analyzes and calculates the weight-bearing area model extracted by the spatial coarse alignment process and the weight-bearing area model extracted after the spatial fine alignment process segmentation, and can calculate the necrosis rate of the femoral head weight-bearing area to the maximum extent in three-dimensional space, greatly improving the accuracy of the subsequent calculated femoral head weight-bearing area necrosis rate, thereby providing more reliable data support for accurately judging the spatial stereoscopic relationship of femoral head necrosis, adapting to the clinical needs of surgical treatment of femoral head necrosis and planning of femoral head necrosis osteotomy surgery, solving the problems in the existing technology that the spatial stereoscopic relationship of femoral head necrosis cannot be accurately judged and the extracted femoral head model has low accuracy, and can further achieve early prediction and prevention of collapse of avascular necrosis of the femoral head.

[0074] Among them, the spatial precise registration process of the original bone model B is used to perform denoising on the weight-bearing area of ​​the femoral head extracted in the previous step and remove noise data, such as Figure 2 As shown, the process is steps S2031 to S2033.

[0075] Step S2031: Use the distance threshold algorithm to traverse each triangular facet mesh of the femoral head weight-bearing area B1 to generate a triangular facet mesh set S 22 .

[0076] The specific steps of using the distance threshold algorithm to traverse the femoral head weight-bearing area B1 are as follows: traverse each triangular facet mesh of the femoral head weight-bearing area B1, extend the three vertices of each triangular facet mesh by 0.013m along the corresponding normal direction to obtain the corresponding line segments, and if the line segments corresponding to the three vertices of the triangular facet mesh do not intersect with the original necrosis model A, then take out the facet and put it in the triangular facet mesh set S 22 middle.

[0077] Step S2032: Use the distance threshold algorithm to traverse the triangle mesh set S 22 Each triangular face mesh in generates a triangular face mesh set S 23 and the triangle mesh set S 24 .

[0078] Among them, the distance threshold algorithm is used to traverse the triangular mesh set S 22 The specific steps are: traverse the triangle mesh set S 22 For each triangular mesh, extend the three vertices of each triangular mesh by 0.002m in the opposite direction of their corresponding normals to obtain their corresponding line segments. If the line segments corresponding to the three vertices of the triangular mesh do not intersect with the original necrotic model A, then take out the face and put it in the triangular mesh set S. 23Otherwise, place it in the triangle mesh set S 24 middle.

[0079] Step S2033: Triangle mesh set S 23 Perform denoising to remove isolated noise points in the set and establish a mesh with the triangle mesh set S 23 The corresponding load-bearing area model B2 and the triangle mesh set S 24 The corresponding necrotic weight-bearing area model A3.

[0080] More specifically, step S203 includes at least a precise spatial registration process for the original necrosis model A. This process is used to merge the necrosis weight-bearing region model A3 extracted in the previous step with the necrosis weight-bearing region model A2 to achieve precise extraction of the necrosis weight-bearing region. This process includes at least step S2034.

[0081] Step S2034: performing a Boolean union operation on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3 to obtain a final necrotic weight-bearing area model A4.

[0082] Preferably, a Boolean union operation is performed on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3 to obtain a final necrotic weight-bearing area model A4. The final necrotic weight-bearing area model A4 is obtained by adding and merging the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3 into one object. Specifically, image processing software, such as the Boolean Command tool in Mimics, can be used to perform a Boolean union operation on the necrotic weight-bearing area model A2 and the necrotic weight-bearing area model A3 to obtain the final necrotic weight-bearing area model A4. Boolean operation is defined as a digital symbolized logical deduction method, including union, intersection, and subtraction. This logical operation method is referenced in graphics processing operations to combine simple basic graphics to produce new shapes, and develop from two-dimensional Boolean operations to Boolean operations of three-dimensional graphics. Boolean union operation is defined as a method used to merge two models, the intersecting parts will be deleted, and after the operation is completed, the two objects will become one object.

[0083] Step S3: Cut the femoral head weight-bearing precise area and the necrosis weight-bearing precise area respectively according to the established cutting surface 12, so as to obtain the corresponding weight-bearing stress concentration areas on the femoral head weight-bearing precise area and the necrosis weight-bearing precise area, and calculate the integrity rate of at least one of the weight-bearing stress concentration areas.

[0084] Wherein, step S3 at least includes the following steps S301 to S303:

[0085] Step S301: cutting the necrotic load-bearing precise area according to the established cutting surface 12, so as to obtain the load-bearing stress concentration area corresponding to the necrotic load-bearing precise area.

[0086] Preferably, the load-bearing stress concentration area on the necrotic weight-bearing area model A4 is the anterior and lateral two-thirds of the necrosis. Step S301 is used to extract the necrotic weight-bearing area in the anterior and lateral two-thirds of the acetabulum. Step S301 can be: according to the gravity line direction 13 corresponding to the precise weight-bearing area of ​​the femoral head, a cutting surface 12 of the necrotic weight-bearing area model A4 is established, based on the fact that the anterior and lateral two-thirds of the necrotic weight-bearing area on the precise weight-bearing area is its load-bearing stress concentration area, the necrotic weight-bearing area model A4 is cut along the cutting surface 12 to extract the necrotic anterior and lateral two-thirds model A5 on the necrotic weight-bearing area model A4.

[0087] Among them, Figure 5 As shown, the cutting surface 12 is the reference surface for distinguishing between the inside and outside of the femoral head. The reference surface is defined as the plane of the hip bone in the acetabulum along the gravity line 13, which can be used to distinguish the inside and outside of the femoral head weight-bearing area. The user can determine the inside and outside of the weight-bearing area by adjusting the position of the reference surface.

[0088] Step S302: cutting the precise weight-bearing area of ​​the femoral head according to the established cutting surface 12, so as to obtain the corresponding weight-bearing stress concentration area on the precise weight-bearing area of ​​the femoral head.

[0089] Preferably, the weight-bearing stress concentration area on the femoral head weight-bearing precise area is the anterolateral two-thirds of the femur. Step S302 is used to extract the femoral weight-bearing area in the lateral front two-thirds of the acetabulum. Step S302 can be: according to the gravity line direction 13 corresponding to the femoral head weight-bearing precise area, a cutting surface 12 of the weight-bearing area model B2 is established, based on the anterolateral two-thirds of the femur on the femoral head weight-bearing precise area being its weight-bearing stress concentration area, the weight-bearing area model B2 is cut along the cutting surface 12 to extract the anterolateral two-thirds of the femur model B3 on the weight-bearing area model B2.

[0090] Step S303: Calculating the integrity rate of at least one of the load-bearing stress concentration areas.

[0091] Preferably, step S303 is used to calculate the integrity rate of the weight-bearing area at the anterolateral two-thirds of the femur. Step S303 specifically comprises: obtaining the surface area S of the necrotic anterolateral two-thirds model A5. A , based on the model B3 at the anterolateral two-thirds of the femur, its surface area S is obtained B , given the surface area S A and surface area S BThe total surface area S of the anterolateral two-thirds weight-bearing area and the integrity rate of the lateral femoral weight-bearing area are formed by the following: , Among them, the surface area S is obtained based on the model B3 at the anterolateral two-thirds of the femur. B The surface area S is the surface area of ​​the intact part of the femur in the anterior two-thirds of the lateral femur, and the surface area S is obtained based on the model A5 of the anterolateral two-thirds of the necrosis. A It is the surface area of ​​necrosis in the anterior two-thirds of the lateral femur.

[0092] Preferably, step S303: obtaining the surface area S of the load stress concentration area corresponding to the necrotic load precise area. A , based on the corresponding load-bearing stress concentration area on the femoral head load-bearing area, its surface area S is obtained. B , given the surface area S A and surface area S B The total surface area S of the anterolateral two-thirds of the weight-bearing area and the integrity rate of the lateral femoral weight-bearing area, where S=S A +S B ,Rate=S B / S. The integrity rate of the lateral weight-bearing area of ​​the femur, Rate, is calculated as the area of ​​the intact surface portion of the femur located in the anterolateral two-thirds of the weight-bearing area relative to the total area of ​​the weight-bearing area in the area. Further preferably, step S303 further includes displaying the surface projection area and the necrotic projection area of ​​the current intact part of the femur through a display 10 or an external device 11, such as a VR display device, that is, projecting the model A5 of the anterolateral two-thirds of the necrotic part and the model B3 of the anterolateral two-thirds of the femur onto the original necrotic model A and the original bone model B, respectively.

[0093] Preferably, the femoral integrity analysis model construction method further includes an extraction process of the acetabulum lunate surface, and the extraction process includes at least the following steps:

[0094] Reconstructing the two-dimensional tomographic image data of the hip joint to establish a three-dimensional mesh model of the hip joint. The three-dimensional mesh model of the hip joint is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system.

[0095] Traverse each vertex of the three-dimensional mesh model of the hip joint, extend each vertex along its corresponding normal direction by 0.013m to obtain its corresponding line segment. If the corresponding line segments on the vertex intersect with the femur and the original necrosis model A, extract the vertex number of the intersection and put it in the vertex number set S. 31 middle;

[0096] Traverse each triangular mesh of the three-dimensional mesh model of the hip joint. If the three vertices of the triangular mesh all belong to the vertex number set S 31 The triangular mesh is extracted and forms a triangular mesh set S 32 ;

[0097] like Figure 5 As shown, for the triangular mesh set S 32 Perform denoising to remove isolated noise points in the set and extract the triangular mesh set S 32 The connected area S33 with the largest number of vertices is determined to be the extracted connected area S 33 It is the lunate surface of the acetabulum on the corresponding side.

[0098] Preferably, a three-dimensional medical imaging surface reconstruction algorithm based on segmentation is used to reconstruct a three-dimensional grid model of the two-dimensional tomographic image data, and a three-dimensional grid model of the femoral head is obtained through threshold screening, interactive segmentation and three-dimensional reconstruction process. The two-dimensional tomographic image data can be obtained by imaging equipment CT and / or imaging equipment MRI. The three-dimensional grid model of the hip joint part is gridded, and the continuous geometric model is discretized to obtain a finite element model corresponding thereto. Preferably, the three-dimensional grid model of the hip joint part is subjected to surface grid division.

[0099] Preferably, the normal corresponding to the vertex mentioned in each embodiment of the present invention is the vertex normal, and the vertex normal is defined as the normal of the three-dimensional vertex P when the coordinates of the vertex P are x, y, z. Preferably, the triangular facet mesh mentioned in each embodiment of the present invention is defined as including at least one fixed edge and two elastically deformable edges, and the fixed edge and the two elastically deformable edges together constitute a rigid body elastic system. The two elastically deformable edges can be used for finite element calculations, for example, by differentiating the triangular facet to an infinitely close unit element, thereby integrating a corresponding model that is infinitely close to the actual hip bone and femur. Example 2

[0100] This embodiment provides a schematic structural diagram of a server 4 applicable to an embodiment of the present invention, as shown in FIG. Figure 4The example shown is only an example and does not limit the functionality and scope of application of the embodiments of the present invention. The server 4 is generally represented by a server 4, which includes at least one processor 6, a storage device 5, and a bus 7 for connecting the processor 6 and the storage device 5. The bus 7 represents one or more of several types of bus 7 structures, including a storage device 5 bus 7 or a storage device 5 controller, a peripheral bus 7, a graphics acceleration port, a processor 6, or a local bus 7 using any of a variety of bus 7 structures. For example, these architectures include, but are not limited to, one or more of an industry standard architecture bus 7, a microchannel architecture bus 7, an enhanced ISA bus 7, and a peripheral component interconnect bus 7. The server 4 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the server 4, including volatile and non-volatile media, removable and non-removable media.

[0101] The storage device 5 may include computer system readable media in the form of volatile memory, such as random access memory and / or cache memory. The server 4 may further include other removable / non-removable, volatile / non-volatile computer system storage media. For example, the storage device 5 may be used to read and write to non-removable, non-volatile magnetic media (hard disk drives). A disk drive for reading and writing to removable non-volatile disks, as well as an optical disk drive for reading and writing to removable non-volatile optical disks, such as read-only disks, digital video disks, or other optical media, may be provided. In these cases, each drive may be connected to the bus 7 via one or more data media interfaces. The storage device 5 may include at least one program product having a set or at least one program module configured to perform the functions of various embodiments of the present invention.

[0102] A program / utility having a set or at least one program module may be stored in, for example, a storage device 5, such program modules including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include the implementation of a network environment. The program modules typically perform the functions and / or methods described in the embodiments of the present invention. The server 4 may also communicate with one or more external devices 11, such as a keyboard, a pointing terminal, a display 10, etc., and may also communicate with one or more terminals that enable a user to interact with the server 4, and / or with any terminal that enables the server 4 to communicate with one or more other computing terminals, such as a network card, a modem, etc. Such communication may be performed through an input / output (I / O) interface. Furthermore, the server 4 may also communicate with one or more networks, such as a local area network, a wide area network, and / or a public network such as the Internet, through a network adapter 8. As Figure 4As shown, the network adapter 8 communicates with other modules of the server 4 via the bus 7. Although not shown in the figure, other hardware and / or software modules may be used in conjunction with the server 4, including at least one or more of the following: microcode, terminal drivers, redundant processors 6, external disk drive arrays, disk array systems, tape drives, and data backup storage systems.

[0103] The processor 6 executes various functional applications and data processing by running the programs stored in the storage device 5, for example, the method for constructing a femoral integrity analysis model provided in an embodiment of the present invention includes at least the following steps:

[0104] S1: reconstructing a CT sequence image of the femoral head to obtain a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system. The three-dimensional mesh model at least includes an original necrosis model A and an original bone model B.

[0105] S2: using a distance threshold algorithm on the acetabulum lunate surface or on the original necrosis model A, performing at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing area and the necrotic weight-bearing area;

[0106] S3: The femoral head weight-bearing area and the necrotic weight-bearing area are cut respectively according to the cutting surface 12 established on the femoral head weight-bearing area to obtain areas on the femoral head weight-bearing area and the necrotic weight-bearing area corresponding to the front two-thirds of the outer side of the acetabulum, and the integrity rate of the weight-bearing area in the outer two-thirds of the femoral head is calculated based on at least one area obtained by cutting. Example 3

[0107] like Figure 3 As shown, the three-dimensional model creation module 1, the three-dimensional model extraction module 2, and the data analysis module 3 provided by the present invention can be respectively a computer-readable storage medium, which stores a computer program. When the program is executed by the processor 6, the femoral integrity analysis model construction method provided in the embodiment of the present invention is implemented, which at least includes the following steps:

[0108] S1: reconstructing a CT sequence image of the femoral head to obtain a three-dimensional mesh model of the femoral head. The three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system. The three-dimensional mesh model at least includes an original necrosis model A and an original bone model B.

[0109] S2: using a distance threshold algorithm on the acetabulum lunate surface or on the original necrosis model A, performing at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B, thereby completing the spatial coarse registration process and the spatial fine registration process in sequence, and extracting the processed femoral head weight-bearing area and the necrotic weight-bearing area;

[0110] S3: The femoral head weight-bearing area and the necrotic weight-bearing area are cut respectively according to the cutting surface 12 established on the femoral head weight-bearing area to obtain areas on the femoral head weight-bearing area and the necrotic weight-bearing area corresponding to the front two-thirds of the outer side of the acetabulum, and the integrity rate of the weight-bearing area in the outer two-thirds of the femoral head is calculated based on at least one area obtained by cutting.

[0111] Computer storage media may be any combination of one or more computer-readable media. Computer-readable media may be computer-readable signal media or computer-readable storage media. Computer-readable storage media may be, for example, but not limited to, electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or components, or any combination thereof. More specific examples of computer-readable storage media may include one or more of the following: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory, a read-only memory, an erasable programmable read-only memory, an optical fiber, a portable compact disk read-only memory, an optical storage device, a magnetic storage device, or any suitable combination thereof. Preferably, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, device or device.

[0112] A computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, which carries computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0113] Program code embodied on a computer readable medium may be transmitted using any appropriate medium, such as wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

[0114] The computer program code for carrying out the operations of the present invention may be written in one or more programming languages, or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or terminal. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network or a wide area network, or may be connected to an external computer.

[0115] It should be noted that the above-described specific embodiments are illustrative only. Those skilled in the art may devise various solutions based on the disclosure of the present invention, and such solutions fall within the scope of the present invention and are intended to be protected by the present invention. Those skilled in the art should understand that the present description and its accompanying drawings are intended to be illustrative only and are not intended to limit the scope of the claims. The scope of protection of the present invention is defined by the claims and their equivalents.

Claims

1. A registration method for a femur master model, comprising at least the following steps: Reconstructing the two-dimensional tomographic image data of the femoral head to establish a three-dimensional mesh model of the femoral head, wherein the three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system, and the three-dimensional mesh model includes at least an original necrosis model A and an original bone model B. The method further comprises: The distance threshold algorithm about the acetabulum lunate surface or the original necrosis model A is used to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B. The distance threshold algorithm is used to limit the position distance between the extended line segment on the triangular facet vertex normal and the acetabulum lunate surface or the original necrosis model, thereby completing the spatial coarse alignment process and the spatial fine alignment process in sequence, and extracting the precise weight-bearing area of ​​the processed femoral head and the precise weight-bearing area of ​​the necrosis.

2. The registration method of the femur master model according to claim 1, characterized in that: The method further includes: prior to executing the spatial coarse registration process and the spatial fine registration process, performing necrosis segmentation preprocessing on the original necrosis model A, and using the triangular facet mesh set corresponding to the necrosis surface model A1 obtained after the preprocessing as the data source for the spatial coarse registration process, wherein, The pretreatment process includes: The distance threshold algorithm is used to traverse the original necrosis model A to extract the non-intersecting overlapping area between the original necrosis model A and the femoral head.

3. The registration method of the femur master model according to claim 2, characterized in that: The method further comprises: The spatial coarse registration process of the original necrosis model A is used to extract the weight-bearing area of ​​the necrotic surface area within the acetabulum; wherein, The spatial coarse registration process of the original necrosis model A includes: The distance threshold algorithm is used to traverse each triangular face mesh of the necrotic surface model A1 to generate a triangular face mesh set S 12 ; Extract the triangular mesh set S 12 The connected region with the largest area is determined to be the necrotic weight-bearing region, and a necrotic weight-bearing region model A2 is established in the necrotic weight-bearing region.

4. The registration method of the femoral master model according to any one of claims 1 to 3, characterized in that: The method further comprises: The spatial coarse registration process of the original bone model B is used to extract the femoral head weight-bearing area of ​​the femoral head surface area located in the acetabulum area; wherein, The spatial coarse registration process of the original bone model B includes: Use the distance threshold algorithm to traverse each triangular mesh of the original bone model B to generate a triangular mesh set S 21 ; Extract the triangular mesh set S 21 The connected area with the largest area is determined to be the femoral head weight-bearing area, and a femoral head weight-bearing area B1 is established within the femoral head weight-bearing area, wherein: The weight-bearing area of ​​the femoral head refers to the area corresponding to the surface of the femoral head projected onto the surface of the femoral head along the direction of the gravity line (13) of the hip bone located in the acetabulum.

5. The registration method of the femur master model according to claim 4, characterized in that: The method further comprises: The spatial fine registration process of the original bone model B and the spatial fine registration process of the original necrosis model A; wherein, The spatial precise registration process of the original bone model B is used to perform denoising on the extracted weight-bearing area of ​​the femoral head and remove noise data; The spatial precise registration process of the original bone model B includes: The distance threshold algorithm is used to traverse each triangular mesh of the femoral head weight-bearing area B1 to generate a triangular mesh set S 22 ; Use the distance threshold algorithm to traverse the triangle mesh set S 22 Each triangular face mesh of generates a triangular face mesh set S 23 and the triangle mesh set S 24 ; For the triangle mesh set S 23 Perform denoising to remove isolated noise points in the set and establish a mesh with the triangle mesh set S 23 The corresponding load-bearing area model B2 and the triangle mesh set S 24 The corresponding necrotic weight-bearing area model A3.

6. The registration method of the femur master model according to claim 5, characterized in that: The spatial precise registration process of the original necrosis model A is used to merge the extracted necrosis weight-bearing area model A3 with the necrosis weight-bearing area model A2 to achieve accurate extraction of the necrosis weight-bearing area; The spatial precise registration process of the original necrosis model A includes: performing a Boolean union operation on the necrosis weight-bearing area model A2 and the necrosis weight-bearing area model A3, and obtaining the final necrosis weight-bearing area model A4 after the operation.

7. The registration method of the femur master model according to claim 6, characterized in that: The method further comprises: The necrotic weight-bearing precise area is cut according to the established cutting surface (12), so as to obtain the corresponding weight-bearing stress concentration area on the necrotic weight-bearing precise area. The weight-bearing stress concentration area on the necrotic weight-bearing area model A4 is the two-thirds of the anterior and lateral part of the necrotic area. Step S302: cutting the femoral head weight-bearing precise region according to the established cutting surface (12), thereby obtaining the corresponding weight-bearing stress concentration region on the femoral head weight-bearing precise region, the weight-bearing stress concentration region on the femoral head weight-bearing precise region being the anterior and lateral two-thirds of the femur, and extracting the femoral weight-bearing region in the anterior two-thirds of the lateral side of the acetabulum; The integrity rate of at least one of the weight-bearing stress concentration areas is calculated to calculate the integrity rate of the weight-bearing area in the anterolateral two-thirds of the femur.

8. The registration method of the femur master model according to claim 7, characterized in that: The method further comprises: The surface area S was obtained based on the model A5 at the anterolateral two-thirds of the necrosis. A , based on the model B3 at the anterolateral two-thirds of the femur, its surface area S is obtained B , given the surface area S A and surface area S B The total surface area S of the anterolateral two-thirds of the weight-bearing area and the integrity rate of the lateral femoral weight-bearing area, where S = S A +S B , The surface area S is obtained based on the model B3 at the anterolateral two-thirds of the femur. B The surface area S is the surface area of ​​the intact part of the femur in the anterior two-thirds of the lateral femur, and the surface area S is obtained based on the model A5 of the anterolateral two-thirds of the necrosis. A It is the surface area of ​​necrosis in the anterior two-thirds of the lateral femur; Among them, the integrity rate of the lateral femoral weight-bearing area is calculated by calculating the area of ​​the intact surface of the femur in the anterolateral two-thirds of the weight-bearing area as a percentage of the total area of ​​the weight-bearing area in this area.

9. A registration system for a femoral original model, comprising at least a three-dimensional model creation module (1), a three-dimensional model extraction module (2) and a data analysis module (3), wherein the three-dimensional model creation module (1) is configured to reconstruct two-dimensional tomographic image data of a femoral head portion to establish a three-dimensional mesh model of the femoral head, wherein the three-dimensional mesh model of the femoral head is composed of a series of triangular facet mesh sets in a three-dimensional coordinate system, and the three-dimensional mesh model includes at least an original necrotic model A and an original bone model B, and is characterized in that: The three-dimensional model extraction module (2) is used to perform at least one distance threshold iterative calculation on the triangular facet mesh sets of the original necrosis model A and the original bone model B respectively using a distance threshold algorithm on the acetabulum lunate surface or on the original necrosis model A, and by using the distance threshold algorithm to limit the position distance between the extended line segment on the triangular facet vertex normal and the acetabulum lunate surface or the original necrosis model, thereby completing the spatial coarse alignment process and the spatial fine alignment process in sequence, and extracting the processed femoral head weight-bearing precise area and the necrosis weight-bearing precise area.

10. The registration system for the femur master model according to claim 9, wherein: The three-dimensional model extraction module (2) performs necrosis segmentation preprocessing on the original necrosis model A before executing the spatial coarse registration process and the spatial fine registration process, and uses the triangular facet mesh set corresponding to the necrosis surface model A1 obtained after the preprocessing as the data source of the spatial coarse registration process, wherein: The three-dimensional model extraction module (2) uses a distance threshold algorithm to traverse the original necrosis model A to extract the area of ​​the non-intersecting overlapping part of the original necrosis model A and the femoral head.

Images