Method and device for calculating joint rate of acetabular cup bone, computer equipment and storage medium

Abstract

The invention discloses a method and a device for calculating an acetabular cup bone fitting rate, computer equipment and a storage medium. The method for calculating the acetabular cup bone fitting rate comprises the following steps: acquiring a pelvis model and an acetabular cup model; extracting an effective bone binding surface of the acetabular cup model; performing Boolean operation on the effective bone binding surface and the pelvis model to obtain an intersection set and a difference set of the effective bone binding surface relative to the pelvis model; and based on the intersection and the difference set, calculating to obtain the bone fitting rate of the effective bone fitting surface relative to the pelvis model. According to the method for calculating the joint rate of the acetabular cup bone, an objective, accurate and efficient method is provided for evaluation of the joint rate of the acetabular cup bone through geometric model reconstruction, effective bone joint surface extraction, Boolean operation and quantitative calculation, quantitative calculation of the joint rate of the acetabular cup bone is achieved, and the joint rate of the acetabular cup bone is evaluated. The problems that in a traditional method, doctor experience is relied on, and subjective judgment accuracy is insufficient are solved, and scientificity and accuracy of operation planning are remarkably improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of information processing, and in particular to a method and device for calculating the bone contact rate of a hip cup, a computer device and a storage medium. BACKGROUND

[0002] The hip cup is one of the key components of hip replacement surgery, used to replace the damaged acetabulum of the human body, cooperate with the femoral prosthesis (femoral stem + ball head), and reconstruct the function of the hip joint. The bone contact rate of the hip cup is an important indicator for measuring the contact degree of the artificial hip cup with the bone bed surface, reflecting its initial mechanical stability and postoperative bone ingrowth potential. High bone contact rate helps to achieve firm biological fixation, reduces the risk of postoperative loosening and revision, and is a key factor affecting the long-term effect of hip replacement surgery. Therefore, during preoperative planning, special attention should be paid to the size and implantation position of the hip cup to improve the bone contact rate of the hip cup after surgery.

[0003] The preoperative planning of traditional hip replacement is based on the CT image of the patient, and the doctor adjusts the hip cup model in the image space, places it in the appropriate position of the patient's acetabular CT reconstruction, and observes the contact degree of the acetabular cup and the bone surface by switching different viewing angles. The above method depends on the experience of the doctor, and the accuracy of the contact degree of the hip cup is not high. SUMMARY

[0004] The purpose of the present application is to overcome the shortcomings of the prior art, provide a method and device for calculating the bone contact rate of a hip cup, a computer device and a storage medium, to solve the technical problem that the contact degree of the traditional hip cup and the bone surface depends on the experience of the doctor, and the accuracy of the contact degree of the hip cup is not high.

[0005] To achieve the above-mentioned purpose, the present application adopts the following technical solutions:

[0006] In a first aspect, the present application provides a method for calculating the bone contact rate of a hip cup, comprising:

[0007] obtaining a pelvis model and a hip cup model;

[0008] extracting the effective bone contact surface of the hip cup model;

[0009] performing a Boolean operation on the effective bone contact surface and the pelvis model to obtain the intersection and difference set of the effective bone contact surface relative to the pelvis model;

[0010] based on the intersection and the difference set, calculating the bone contact rate of the effective bone contact surface relative to the pelvis model.

[0011] In a second aspect, the present application provides a device for calculating the bone contact rate of a hip cup, comprising:

[0012] An acquisition module is configured to acquire a pelvis model and an acetabular cup model;

[0013] An extraction module is configured to extract an effective bone-fitting surface of the acetabular cup model;

[0014] A Boolean operation module is configured to perform a Boolean operation on the effective bone-fitting surface and the pelvis model to obtain an intersection and a difference set of the effective bone-fitting surface relative to the pelvis model;

[0015] A bone-fitting rate calculation module is configured to calculate a bone-fitting rate of the effective bone-fitting surface relative to the pelvis model based on the intersection and the difference set.

[0016] In a third aspect, the present application provides a computer device, which comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method for calculating the acetabular cup bone-fitting rate according to the first aspect when executing the computer program.

[0017] In a fourth aspect, the present application provides a storage medium, which stores a computer program, and the computer program comprises program instructions executable by a processor to implement the method for calculating the acetabular cup bone-fitting rate according to the first aspect.

[0018] The method for calculating the acetabular cup bone-fitting rate according to the present application provides an objective, accurate and efficient method for evaluating the acetabular cup bone-fitting rate through geometric model reconstruction, effective bone-fitting surface extraction, Boolean operation and quantitative calculation, realizes quantitative calculation of the bone-fitting rate, overcomes the problem of insufficient accuracy of subjective judgment in the traditional method which depends on the experience of doctors, and significantly improves the scientificity and accuracy of surgical planning.

[0019] The above description is only a summary of the technical solutions of the present application, in order to more clearly understand the technical means of the present application, the content of the specification can be implemented, and in order to make the above and other purposes, characteristics and advantages of the present application more obvious and easy to understand, the following preferred embodiments are described in detail as follows. BRIEF DESCRIPTION OF DRAWINGS

[0020] Figure 1 A flowchart of the method for calculating the acetabular cup bone-fitting rate according to the embodiments of the present application;

[0021] Figure 2 A first sub-flowchart of the method for calculating the acetabular cup bone-fitting rate according to the embodiments of the present application;

[0022] Figure 3 A second sub-flowchart of the method for calculating the acetabular cup bone-fitting rate according to the embodiments of the present application;

[0023] Figure 4A third sub-flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application;

[0024] Figure 5 A fourth sub-flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application;

[0025] Figure 6 A fifth sub-flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application;

[0026] Figure 7 A sixth sub-flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application;

[0027] Figure 8 A schematic diagram of the device for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application;

[0028] Figure 9 A schematic block diagram of a computer device provided by the embodiment of the present application. DETAILED DESCRIPTION

[0029] In order to make the objects, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0030] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings of the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the present application.

[0031] In the description of the present application, it should be understood that the terms “center”, “longitudinal”, “transverse”, “length”, “width”, “thickness”, “upper”, “lower”, “front”, “rear”, “left”, “right”, “resin”, “horizontal”, “top”, “bottom”, “inner”, “outer”, “clockwise”, “counterclockwise” and the like indicate the orientation or positional relationship described based on the drawings, and are only for the convenience of describing the present application and simplifying the description, and therefore cannot be understood as indicating or implying that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as limiting the present application.

[0032] In addition, the terms “first” and “second” are only for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the technical features referred to. Therefore, the features defined with “first” and “second” can explicitly or implicitly include one or more of the features. In the description of the present application, the meaning of “multiple” is two or more, unless otherwise specifically limited.

[0033] In the present application, unless specifically defined and limited otherwise, the terms "mount", "connect", "connection", "fixed", and the like, should be construed broadly and can include direct connection, or indirect connection through an intermediate medium, mechanical connection, or electrical connection, direct connection, or indirect connection through an intermediate medium, or the internal communication of two elements or the interaction relationship between two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0034] In the present application, unless specifically defined and limited otherwise, the first feature "on" or "under" the second feature can include that the first and second features are in direct contact, or that the first and second features are not in direct contact but are in contact through another feature between them. Moreover, the first feature "on", "above" and "above" the second feature includes that the first feature is directly above and obliquely above the second feature, or only indicates that the horizontal height of the first feature is higher than that of the second feature. The first feature "under", "below" and "below" the second feature includes that the first feature is directly below and obliquely below the second feature, or only indicates that the horizontal height of the first feature is less than that of the second feature.

[0035] In the description of the present application, the description of the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present application, the illustrative description of the above terms should not be understood as necessarily referring to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0036] Please refer to Figure 1 , Figure 1 The flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application. The embodiment of the present application provides a method for calculating the bone fitting rate of the acetabular cup, which comprises S100-S400:

[0037] S100, obtaining a pelvic model and an acetabular cup model.

[0038] As shown in Figure 2 , Figure 2 The first sub-flowchart of the method for calculating the bone fitting rate of the acetabular cup of the embodiment of the present application, the step S100 comprises S110-S120:

[0039] S110, obtaining the pelvic image data of the patient by CT scanning.

[0040] CT scan is a medical imaging technology based on X-ray imaging principle, which can perform tomographic scanning on internal structures of human body and generate a large amount of two-dimensional tomographic image data, which contains information such as the shape and density of the pelvis.

[0041] S120, reconstructing the pelvis model based on the pelvis image data.

[0042] The three-dimensional reconstruction based on the pelvis image data can adopt a manual segmentation technique or a deep learning method. The manual segmentation technique requires a doctor or a professional to manually outline the contour of the pelvis on each tomographic image of the CT image, and then connects the contours to form a three-dimensional model through a computer algorithm. The deep learning method (such as U-Net) trains the neural network using a large amount of labeled CT image data, so that the model can automatically identify and segment the pelvis region, and then reconstruct a three-dimensional model. Alternatively, the three-dimensional reconstruction based on the pelvis image data can also be based on manual segmentation, combined with semi-automatic segmentation technology, such as using image processing algorithms to automatically identify the initial contour, and then corrected by the doctor, to improve the segmentation efficiency and accuracy.

[0043] In some embodiments, as shown in Figure 3 , Figure 3 is a second sub-flow chart of the method for calculating the acetabular cup bone fitting rate of the embodiment of the application, and the step S100 further includes S130-S150:

[0044] S130, selecting the acetabular cavity region in the pelvis model and fitting to obtain the inner diameter of the acetabular cavity region.

[0045] Specifically, this step selects the acetabular cavity region in the reconstructed pelvis model and fits the inner diameter of the acetabular cavity region through a geometric fitting algorithm to obtain the inner diameter size. It should be explained that the operation of selecting the acetabular cavity region can be manually selected by a doctor.

[0046] S140, according to the inner diameter of the acetabular cavity region, screening the acetabular cup meeting the preset condition from the medical consumable database.

[0047] The medical consumable database is provided by a medical institution or an acetabular cup supplier, which should contain detailed acetabular cup parameter information, including but not limited to brand, material (such as metal, ceramic, polyethylene, etc.), available fixation method (such as biological fixation, cement fixation, etc.), inventory quantity, etc., and key size parameters of the acetabular cup, especially the outer diameter range, so as to match the inner diameter of the acetabular cavity.

[0048] S150, displaying the parameter information of the acetabular cup meeting the preset condition, the parameter information including brand, material, available fixation method, inventory quantity and the difference between the outer diameter of the acetabular cup and the inner diameter of the acetabular cavity.

[0049] Among them, in the step of screening out an acetabular cup with an outer diameter matching the inner diameter of the acetabular fossa obtained by fitting in the medical consumables database, a preset condition can be set that the difference between the outer diameter of the acetabular cup and the inner diameter of the acetabular fossa is within 1mm-3mm as a screening condition. Then, acetabular cups with a difference between the outer diameter of the acetabular cup and the inner diameter of the acetabular fossa within 1mm-3mm are all acetabular cups that meet the preset conditions, so as to ensure that the acetabular cup can be successfully implanted and maintain a certain stability.

[0050] Furthermore, in this embodiment, the parameter information of the acetabular cup that meets the preset conditions is displayed to facilitate the doctor to select the acetabular cup that meets the actual situation of the patient.

[0051] S200: Extracting the effective bone-fitting surface of the acetabular cup model.

[0052] It is understandable that not all parts of the acetabular cup need to fit with the bone. Its effective bone fitting surface is only the outermost hemispherical shell part. In order to accurately calculate the bone fitting rate of the acetabular cup, it is necessary to extract its effective bone fitting surface, that is, the outermost hemispherical shell part of the acetabular cup.

[0053] like Figure 4 As shown, Figure 4 This is a first sub-flowchart of the method for calculating the acetabular cup bone fitting rate according to an embodiment of the present invention, wherein step S200 includes steps S210 to S230:

[0054] S210: Calculate the normal vectors of all triangular facets constituting the acetabular cup model.

[0055] It's understandable that determining the directional characteristics of each triangular facet of the acetabular cup model provides foundational data for subsequently distinguishing different geometric feature regions, thereby accurately identifying the effective bone-fitting surface. Because the effective bone-fitting surface of the acetabular cup model has specific directional characteristics, calculating the normal vector can capture these characteristic differences, making it a key prerequisite for extracting the effective bone-fitting surface.

[0056] Specifically, each triangular facet in the three-dimensional model can be regarded as a plane, and the normal vector thereof can be calculated by the three vertex coordinates of the triangular facet, specifically, by performing a cross multiplication operation on vectors formed by any two sides of the triangular facet to obtain a vector perpendicular to the plane of the triangular facet, i.e., the normal vector. For example, for three vertexes A(x1, y1, z1), B(x2, y2, z2), and C(x3, y3, z3) on a triangular facet, vector AB=(x2-x1, y2-y1, z2-z1), and vector AC=(x3-x1, y3-y1, z3-z1), the normal vector Q=ABxAC=((y2-y1)(z3-z1)-(z2-z1)(y3-y1), (z2-z1)(x3-x1)-(x2-x1)(z3-z1), (x2-x1)(y3-y1)-(y2-y1)(x3-x1)).

[0057] S220, the triangular facets adjacent to each other and having an included angle of normal vectors less than or equal to a preset threshold are set as the same grid data, and the triangular facets adjacent to each other and having an included angle of normal vectors greater than the preset threshold are set as mutually independent grid data.

[0058] It can be understood that, according to the geometric characteristics of the acetabular cup model, the adjacent triangular facets having similar directions are grouped into the same region, and the adjacent triangular facets having a large direction difference are separated, so that the acetabular cup model is split into multiple grid data parts having different geometric characteristics, facilitating screening of a region where an effective bone fitting surface is located and effectively removing interference of irrelevant structures.

[0059] In this embodiment, the included angle of normal vectors of adjacent triangular facets can be calculated by a dot product formula. For normal vectors Q1 and Q2, the included angle θ satisfies cosθ=(Q1·Q2) / (|Q1||Q2|). The preset threshold can be selected as 5°. When the included angle of normal vectors of adjacent triangular facets is less than or equal to the preset threshold, it is considered that they belong to the same geometric characteristic region and are merged into the same grid data. Otherwise, if the included angle is greater than the threshold, they are set as mutually independent grid data. By traversing all pairs of adjacent triangular facets and performing merging or independent processing according to the above rule, the acetabular cup model is finally split into multiple grid data parts.

[0060] S230, screening the grid data with the largest spatial volume from the multiple grid data, and determining that the grid data is an effective bone fitting surface of the acetabular cup model.

[0061] In this embodiment, the effective bone fitting surface of the acetabular cup is the outermost hemispherical shell part, which usually occupies a large spatial volume in the model. By screening the grid data with the largest spatial volume, the effective bone fitting surface of the acetabular cup can be accurately determined, and subsequent bone fitting rate calculation can be based on the correct region.

[0062] S300: Perform a Boolean operation on the effective bone fitting surface and the pelvic model to obtain an intersection and a difference between the effective bone fitting surface and the pelvic model.

[0063] It can be understood that the Boolean operation can determine the spatial relationship between the effective bone fit surface of the acetabular cup and the pelvic model, specifically finding the intersection (i.e., the area of ​​actual fit) and the difference (i.e., the area of ​​non-fit) between the two, providing a direct geometric basis for the subsequent accurate calculation of the bone fit rate. In particular, as the doctor manipulates the fit position of the acetabular cup model and the pelvic model, the intersection and difference results of the Boolean operation can be updated in real time, thereby also reflecting the bone fit rate of the acetabular cup model and the pelvic model in real time, which helps guide the doctor in planning treatment plans.

[0064] In some embodiments, specialized mesh data processing libraries, such as Meshlib, CGAL, and libigl, can be used to perform precise Boolean operations on complex three-dimensional mesh models. Using the extracted acetabular cup effective bone contact surface and the pelvic model as input, the intersection (regions belonging to both the acetabular cup effective bone contact surface and the pelvic model) and difference (regions belonging to the acetabular cup effective bone contact surface but not the pelvic model) between the two can be obtained.

[0065] S400: Calculate and obtain a bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference.

[0066] Among them, through the quantitative calculation of the intersection and difference areas, the numerical index of the acetabular cup bone fitting rate can be obtained, thereby realizing the quantitative evaluation of the degree of fit between the acetabular cup and the pelvis, providing doctors with an objective and accurate quantitative basis for surgical planning, and improving the scientificity and rationality of surgical planning.

[0067] In the first embodiment, as Figure 5 As shown, Figure 5 This is a fourth sub-flowchart of the method for calculating the acetabular cup bone fitting rate according to an embodiment of the present invention, wherein step S400 includes steps S401-S403:

[0068] S401 : Optimize the intersection and the difference to make the mesh vertices of the intersection and the difference uniformly distributed.

[0069] It is understandable that the intersection and difference after Boolean operations may contain uneven distribution of mesh vertices, such as vertices in some areas are too dense and vertices in some areas are too sparse. This unevenness may cause the vertex count statistics to not accurately reflect the actual fit. By optimizing the vertex distribution to make it uniform, it can ensure that the vertex count calculation is more representative, thereby improving the accuracy of subsequent bone fit rate calculations.

[0070] In some embodiments, mesh optimization methods such as Delaunay triangulation and Voronoi diagrams can be used to achieve uniform vertex distribution. Delaunay triangulation reconnects existing vertices so that the circumcircle of each triangle contains no other vertices, thereby generating a relatively uniform triangular mesh. Voronoi diagrams divide the space into regions so that the area surrounding each vertex is relatively uniform in size.

[0071] S402 : Obtain the number N of mesh vertices in the intersection set after optimization processing, and the number M of mesh vertices in the difference set after optimization processing.

[0072] After the vertex distribution is uniform, by counting the number of mesh vertices in the intersection and difference sets, a simple and effective quantitative basis for calculating the bone fitting rate can be provided. Specifically, the number of vertices can reflect the size ratio of the fitting area and the non-fitting area to a certain extent. Especially when the vertex distribution is uniform, the number of vertices has a certain correlation with the area of ​​the region. Therefore, it can be used as an approximate indicator for the calculation of the bone fitting rate, providing doctors with fast and intuitive evaluation results. By traversing the mesh data of the intersection and difference sets, the number of vertices contained therein can be counted respectively to obtain the values ​​of N and M. In some embodiments, the vertex list in the mesh data structure or the index information of the vertex can be used for counting. This is a relatively direct and computationally low-cost operation that can quickly obtain the required quantitative data.

[0073] S403: Calculate the bone fitting rate of the effective bone fitting surface relative to the pelvic model as N / (N+M).

[0074] It can be understood that, under the premise of uniform vertex distribution, the number of vertices N in the intersection represents the relative size of the fitted area, and the number of vertices M in the difference represents the relative size of the non-fitted area. Adding the two together gives the total number of vertices N+M, which represents the total size of the contact area between the entire effective bone fitting surface and the pelvic bone. By calculating the ratio of N to the total number of vertices N / (N+M), a value between 0 and 1 can be obtained. This value is the bone fitting rate, and a larger value indicates a higher degree of fit.

[0075] In this first embodiment, intuitive evaluation metrics are provided to the surgeon, enabling them to quickly understand the fit between the acetabular cup and the pelvis, thereby assisting in surgical planning decisions. Furthermore, due to the high computational efficiency of the vertex-based calculation method, calculations can be completed and results provided in a relatively short time, improving the efficiency and convenience of the entire surgical planning process.

[0076] In the second embodiment, as Figure 6 As shown,Figure 6 The fifth sub-flow chart for calculating the bone fit rate of the acetabular cup in the embodiment of the present application is shown. The steps S400 include S411-S412.

[0077] S411, calculate the area S1 of the intersection region and the area S2 of the difference set region.

[0078] It can be understood that the area of each triangular facet in the three-dimensional model can be calculated by the geometric formula of the triangle, that is, half of the length of the vector obtained by the cross product of any two sides of the triangle.

[0079] In some embodiments, for the intersection and difference set regions, the total areas S1 and S2 can be obtained by accumulating the areas of all triangular facets in them, respectively. Specifically, for each triangular facet, two vectors are determined by the coordinates of the three vertices, the length of the cross product of the two vectors is calculated, and then divided by 2 to obtain the area of the triangular facet. Then, all the triangular facets in the intersection or difference set are traversed, and their areas are accumulated to obtain the total area of the corresponding region.

[0080] S412, calculate the bone fit rate of the effective bone fit surface relative to the pelvic model as S1 / (S1+S2).

[0081] It can be understood that the area S1 of the intersection region represents the part actually fitted by the acetabular cup and the pelvic, and the area S2 of the difference set region represents the part not fitted. Adding them together to obtain the total area S1+S2 represents the total size of the entire effective bone fit surface and the pelvic contact area. By calculating the ratio of S1 to the total area S1 / (S1+S2), a value between 0 and 1 is obtained, which is the bone fit rate, which can directly reflect the fitting degree of the acetabular cup and the pelvic.

[0082] In the second embodiment, the area is taken as the quantitative indicator to provide an accurate and intuitive geometric basis for the calculation of the bone fit rate. The area calculation can directly reflect the size relationship between the fitted region and the unfitted region, and the calculation result has a clear geometric meaning, which is easy for doctors to understand and apply, so that the quantitative result of the acetabular cup bone fit rate is obtained in a intuitive and accurate way, providing a reliable evaluation basis for doctors, so that they can quickly and accurately understand the fitting condition of the acetabular cup and the pelvic, thereby assisting them in making surgical planning decisions, and providing strong support for improving the quality and success rate of surgery.

[0083] In the third embodiment, as shown in Figure 7 , Figure 7 The sixth sub-flow chart for calculating the bone fit rate of the acetabular cup in the embodiment of the present application is shown. The steps S400 include S421-S426:

[0084] S421、optimizing the intersection and the difference set, so that the grid vertex distribution of the intersection and the difference set is homogenized;

[0085] S422、obtaining the number N of grid vertices in the optimized intersection and the number M of grid vertices in the optimized difference set;

[0086] S423、calculating the first fitting rate of the effective bone fitting surface relative to the pelvis model as N / (N+M);

[0087] S424、calculating the area S1 of the intersection region and the area S2 of the difference set region;

[0088] S425、calculating the second fitting rate of the effective bone fitting surface relative to the pelvis model as S1 / (S1+S2);

[0089] S426、calculating the average of the first fitting rate and the second fitting rate to obtain the bone fitting rate of the effective bone fitting surface relative to the pelvis model.

[0090] It can be understood that the method of obtaining the first fitting rate is the same as the method of obtaining the bone fitting rate in the first embodiment, and the method of obtaining the second fitting rate is the same as the method of obtaining the bone fitting rate in the second embodiment. By comprehensively calculating the first fitting rate and the second fitting rate based on the number of vertices and the area, and taking the average as the final bone fitting rate, the advantages of the two methods are combined to improve the comprehensiveness and accuracy of the calculation result, and a more comprehensive and accurate bone fitting rate evaluation is obtained.

[0091] As shown in Figure 8 , the embodiment of the application also provides a device 500 for calculating the bone fitting rate of the acetabular cup, which comprises an acquisition module 501, an extraction module 502, a Boolean operation module 503, and a bone fitting rate calculation module 504.

[0092] The acquisition module 501 is used to acquire the pelvis model and the acetabular cup model.

[0093] The extraction module 502 is used to extract the effective bone fitting surface of the acetabular cup model.

[0094] The Boolean operation module 503 is used to perform Boolean operation on the effective bone fitting surface and the pelvis model to obtain the intersection and the difference set of the effective bone fitting surface relative to the pelvis model.

[0095] The bone fitting rate calculation module 504 is used to calculate the bone fitting rate of the effective bone fitting surface relative to the pelvis model based on the intersection and the difference set.

[0096] In some embodiments, the extraction module 502 is further configured to calculate normal vectors of all triangular facets constituting the acetabular cup model; set the triangular facets adjacent to each other and having an included angle of normal vectors less than or equal to a preset threshold as the same mesh data, and set the triangular facets adjacent to each other and having an included angle of normal vectors greater than the preset threshold as mesh data independent of each other; and select the mesh data with the largest spatial volume from the plurality of mesh data, and determine the mesh data as the effective bone-fitting surface of the acetabular cup model.

[0097] It should be noted that the specific implementation process of the above device and each unit can be clearly understood by those skilled in the art, and the corresponding description in the foregoing method embodiments can be referred to. For the convenience and brevity of description, it will not be repeated here.

[0098] The above device for calculating the acetabular cup bone-fitting rate can be implemented in the form of a computer program, which can run on a computer device as shown in the accompanying drawings. Figure 9 The computer device 600 includes a processor 602, a memory, and a network interface 605 connected through a system bus 601, wherein the memory can include a non-volatile storage medium 603 and an internal memory 604.

[0099] Referring to Figure 9 , the computer device 600 includes a processor 602, a memory, and a network interface 605 connected through a system bus 601, wherein the memory can include a non-volatile storage medium 603 and an internal memory 604.

[0100] The non-volatile storage medium 603 can store an operating system 6031 and a computer program 6032. The computer program 6032 includes program instructions, which, when executed, can cause the processor 602 to perform a method for calculating the acetabular cup bone-fitting rate: obtaining a pelvic model and an acetabular cup model; extracting an effective bone-fitting surface of the acetabular cup model; performing a Boolean operation on the effective bone-fitting surface and the pelvic model to obtain an intersection and a difference set of the effective bone-fitting surface relative to the pelvic model; and calculating a bone-fitting rate of the effective bone-fitting surface relative to the pelvic model based on the intersection and the difference set.

[0101] The processor 602 is configured to provide computing and control capabilities to support the operation of the entire computer device 600.

[0102] The internal memory 604 provides an environment for the running of the computer program 6032 in the non-volatile storage medium 603. When the computer program 6032 is executed by the processor 602, the processor 602 can perform a method for calculating the acetabular cup bone-fitting rate.

[0103] The network interface 605 is configured to perform network communication with other devices. Those skilled in the art can understand that Figure 9The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device 600 to which the scheme of the present application is applied. Specifically, the computer device 600 can include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0104] The processor 602 is configured to run the computer program 6032 stored in the memory to implement the following steps S100-S400:

[0105] S100, acquire a pelvis model and an acetabular cup model;

[0106] S200, extract an effective bone fitting surface of the acetabular cup model;

[0107] S300, perform Boolean operation on the effective bone fitting surface and the pelvis model to obtain an intersection and a difference set of the effective bone fitting surface relative to the pelvis model;

[0108] S400, based on the intersection and the difference set, calculate a bone fitting rate of the effective bone fitting surface relative to the pelvis model.

[0109] It should be understood that, in the embodiments of the present application, the processor 602 can be a central processing unit (CPU), and the processor 602 can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0110] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiments of the method can be completed by a computer program instructing related hardware. The computer program includes program instructions, and the computer program can be stored in a storage medium, which is a computer-readable storage medium. The program instructions are executed by at least one processor in the computer system to implement the process steps of the above-mentioned embodiments of the method.

[0111] Therefore, the application further provides a storage medium. The storage medium can be a computer readable storage medium. The storage medium stores a computer program, wherein the computer program includes program instructions, and the program instructions, when executed by a processor, can implement the method for calculating the bone fitting rate of the acetabular cup. The storage medium stores a computer program, and the computer program includes program instructions, and the program instructions, when executed by a processor, can implement the method. The program instructions include the following steps:

[0112] S100, acquiring a pelvic model and an acetabular cup model;

[0113] S200, extracting an effective bone fitting surface of the acetabular cup model;

[0114] S300, performing a Boolean operation on the effective bone fitting surface and the pelvic model to obtain an intersection and a difference set of the effective bone fitting surface relative to the pelvic model;

[0115] S400, calculating a bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference set.

[0116] The storage medium can be a U disk, a mobile hard disk, a read-only memory (ROM), a magnetic disk or an optical disk, and various computer readable storage media that can store program codes.

[0117] Those skilled in the art can realize that the units and algorithm steps of the examples described in combination with the embodiments disclosed herein can be realized by electronic hardware, computer software or a combination of both. In order to clearly illustrate the interchangeability of hardware and software, the components and steps of the examples have been described in general terms in the above description. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. A person skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the application.

[0118] In several embodiments provided by the application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of each unit is only a logical function division, and actual implementation can have another division manner. For example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed.

[0119] The steps in the method of the embodiments of the present application can be adjusted in sequence, combined and deleted according to actual needs. The units in the device of the embodiments of the present application can be combined, divided and deleted according to actual needs. In addition, each functional unit in each embodiment of the present application can be integrated in one processing unit, or each unit can exist physically independently, or two or more units can be integrated in one unit.

[0120] The integrated unit, if realized in the form of a software functional unit and sold or used as an independent product, can be stored in a storage medium. Based on such understanding, the technical solutions of the present application, essentially or in the contribution to the prior art, or all or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes several instructions for causing a computer device (which can be a personal computer, a terminal, or a network device, etc.) to execute all or part of the steps of the methods described in the embodiments of the present application.

[0121] The above is only used to further illustrate the technical content of the present application by way of examples, so that the reader can be easier to understand, but does not represent the embodiments of the present application are limited to this, any technical extension or re-creation made according to the present application, is protected by the present application. The scope of protection of the present application is subject to the claims.

Claims

1. A method for calculating the acetabular cup bone fit rate, characterized in that: include: Obtain a pelvic model and acetabular cup model; extracting the effective bone-fitting surface of the acetabular cup model; Performing a Boolean operation on the effective bone fitting surface and the pelvic model to obtain an intersection and a difference between the effective bone fitting surface and the pelvic model; Based on the intersection and the difference, a bone fitting rate of the effective bone fitting surface relative to the pelvic model is calculated.

2. The method for calculating the acetabular cup bone fit rate according to claim 1, wherein: The step of extracting the effective bone fitting surface of the acetabular cup model comprises: Calculating normal vectors of all triangular facets constituting the acetabular cup model; The adjacent triangular facets whose normal vector angle is less than or equal to a preset threshold are set as the same mesh data, and the adjacent triangular facets whose normal vector angle is greater than the preset threshold are set as independent mesh data; The grid data with the largest spatial volume is screened out from the plurality of grid data, and is determined to be the effective bone fitting surface of the acetabular cup model.

3. The method for calculating the acetabular cup bone fitting rate according to claim 1, wherein: The step of calculating the bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference includes: Optimizing the intersection and the difference to make the grid vertices of the intersection and the difference uniformly distributed; Obtaining the number N of mesh vertices in the intersection after optimization processing, and the number M of mesh vertices in the difference after optimization processing; The bone fitting rate of the effective bone fitting surface relative to the pelvic model is calculated as N / (N+M).

4. The method for calculating the acetabular cup bone fitting rate according to claim 1, wherein: The step of calculating the bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference includes: Calculating the area S1 of the intersection region and the area S2 of the difference region; The bone fitting rate of the effective bone fitting surface relative to the pelvic model is calculated as S1 / (S1+S2).

5. The method for calculating the acetabular cup bone fitting rate according to claim 1, wherein: The step of calculating the bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference includes: Optimizing the intersection and the difference to make the grid vertices of the intersection and the difference uniformly distributed; Obtaining the number N of mesh vertices in the intersection after optimization processing, and the number M of mesh vertices in the difference after optimization processing; Calculating a first fitting ratio of the effective bone fitting surface relative to the pelvic model as N / (N+M); Calculating the area S1 of the intersection region and the area S2 of the difference region; Calculating a second fitting ratio of the effective bone fitting surface relative to the pelvic model as S1 / (S1+S2); An average value of the first fitting rate and the second fitting rate is calculated to obtain the bone fitting rate of the effective bone fitting surface relative to the pelvic model.

6. The method for calculating the acetabular cup bone fitting rate according to claim 1, wherein: The step of obtaining the pelvic model comprises: Obtain pelvic imaging data of the patient through CT scanning; The pelvic model is reconstructed based on the pelvic image data.

7. A device for calculating acetabular cup bone fit rate, characterized in that: include: An acquisition module, used for acquiring a pelvis model and an acetabular cup model; An extraction module, configured to extract an effective bone fitting surface of the acetabular cup model; a Boolean operation module, configured to perform a Boolean operation on the effective bone fitting surface and the pelvic model to obtain an intersection and a difference between the effective bone fitting surface and the pelvic model; A bone fitting rate calculation module is used to calculate the bone fitting rate of the effective bone fitting surface relative to the pelvic model based on the intersection and the difference.

8. The device for calculating the acetabular cup bone fitting rate according to claim 7, characterized in that: The extraction module is also used to calculate the normal vectors of all triangular facets that constitute the acetabular cup model; set the adjacent triangular facets whose normal vector angle is less than or equal to a preset threshold as the same grid data, and set the adjacent triangular facets whose normal vector angle is greater than the preset threshold as independent grid data; filter out the grid data with the largest spatial volume from the multiple grid data, and determine it as the effective bone bonding surface of the acetabular cup model.

9. A computer device, characterized in that: The method comprises a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for calculating the acetabular cup bone fitting rate according to any one of claims 1 to 6 is implemented.

10. A storage medium, characterized in that: The storage medium stores a computer program, which includes program instructions. When the program instructions are executed by a processor, the method for calculating the acetabular cup bone fitting rate according to any one of claims 1 to 6 can be implemented.

Images