Hip joint prosthesis position selection method and device based on mechanical simulation

Abstract

The application provides a hip joint prosthesis position selection method and device based on mechanical simulation, and belongs to the technical field of prostheses. The method comprises the following steps: obtaining CT images to obtain CT images of each hip joint sub-region; obtaining CT image analysis results of each hip joint sub-region by processing, and performing corresponding image enhancement processing; obtaining a whole hip joint region CT image, analyzing the whole hip joint region CT image to obtain imaging quality indicators of the whole hip joint region CT image, and judging the whole imaging verification result; based on the whole imaging verification result, if the whole imaging verification result is qualified, the whole hip joint region CT image is input into a three-dimensional reconstruction module in a mechanical simulation platform to perform hip joint prosthesis implantation simulation testing; and the hip joint prosthesis implantation simulation testing data is counted and processed to obtain a corrected region image set, and the prosthesis position simulation analysis result of the whole hip joint region CT image is obtained after processing the corrected region image set.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of prostheses, in particular to a hip joint prosthesis position selection method and device based on mechanical simulation. BACKGROUND

[0002] With the continuous progress of science and technology, hip joint prosthesis technology has made remarkable development and is widely used in clinical practice. The placement position and angle of the hip joint prosthesis are crucial for the postoperative mechanical balance and motor function recovery of patients. Therefore, targeted position selection is needed before the placement of the hip joint prosthesis.

[0003] For example, the control method of the joint replacement surgery robot disclosed in the patent CN114305697B includes the following steps: capturing the patient's bone image and transmitting the image to the computer system, the computer system obtains the image data of the bone; selecting a suitable prosthesis model in the computer system and placing it on the diseased joint for matching. After the matching is completed, the overlapping part of the joint prosthesis model and the bone is the part that needs to be removed and replaced. A signal source is installed as a signal sending end. The patient's bone image is captured again, and the range of the image includes the signal source and the joint part. In the computer system, the position of the signal source is taken as the origin to establish a coordinate system. The computer system automatically calculates and memorizes the interface between the joint prosthesis model and the bone overlapping part and the coordinates in the coordinate system. The surgical robot receives the coordinates and combines the received signal source to emit a signal. The position coordinates of the diseased joint that needs to be removed in the coordinate system are obtained, and the surgical robot performs surgical operation.

[0004] In the prior art, the patient's bone image is captured and recognized to match the lesion joint, and the position is selected directly using the bone image. However, in actual application, the quality of the bone image is often uneven, and the quality of the bone image can seriously affect the accuracy of positioning. Therefore, the existing technology has the problem of unreasonable prosthesis position selection due to unqualified image quality processing. SUMMARY

[0005] The present application provides a hip joint prosthesis position selection method and device based on mechanical simulation, which solves the problem of unreasonable prosthesis position selection due to unqualified image quality processing in the prior art, and realizes the accuracy and controllability of prosthesis position selection.

[0006] To solve the above-mentioned application purposes, the technical solutions provided by the present application are as follows:

[0007] The embodiment of the application provides a hip joint prosthesis position selection method based on mechanical simulation, which comprises the following steps: S1, acquiring a CT image when a CT upload signal is received by a mechanical simulation platform, segmenting the CT image, obtaining a hip joint region image, and segmenting the hip joint region image to obtain CT images of each hip joint sub-region; S2, acquiring CT image performance parameters of each hip joint sub-region, analyzing the CT image performance parameters to obtain CT image performance indexes of each hip joint sub-region, processing the CT image performance indexes to obtain CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; S3, acquiring the CT images of each hip joint sub-region after the image enhancement processing, performing splicing processing on the CT images, obtaining a whole hip joint region CT image, analyzing the whole hip joint region CT image to obtain a whole hip joint region CT image quality index, and judging a whole imaging verification result according to the whole hip joint region CT image quality index, wherein the whole imaging verification result comprises a whole imaging verification pass and a whole imaging verification fail; S4, discriminating based on the whole imaging verification result, if the whole imaging verification result is the whole imaging verification pass, inputting the whole hip joint region CT image into a three-dimensional reconstruction module in the mechanical simulation platform to perform hip joint prosthesis implantation simulation inspection; and S5, processing hip joint prosthesis implantation simulation inspection data to obtain a corrected region image set, and processing the corrected region image set to obtain a prosthesis position simulation analysis result of the whole hip joint region CT image.

[0008] Optionally, the CT image performance parameters of each hip joint sub-region are acquired, and the CT image performance parameters of each hip joint sub-region comprise average gray values, average densities, noise values, resolutions and contrasts of the hip joint sub-regions; a preset CT image reference set in a database is acquired, and the CT image reference set is compared with the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indexes of each hip joint sub-region; and the CT image performance indexes of each hip joint sub-region are used to represent a CT image quality conformity degree of each hip joint sub-region.

[0009] Optionally, the CT image performance indexes of each hip joint sub-region are obtained by the following method:

[0010]

[0011] In the formula, BX i represents the CT image performance index of the i th hip joint sub-region, i represents the number of the hip joint sub-region, i = 1, 2,..., i max , i max represents the total number of the hip joint sub-regions, e represents a natural constant, HD irepresents the average gray value of the i-th hip joint sub-region, and ΔHD represents the average gray value reference, and ρ i represents the average density of the i-th hip joint sub-region, and Δρ represents the average density reference, and ZS i represents the noise value of the i-th hip joint sub-region, and ΔZS represents the noise reference, and FB i represents the resolution of the i-th hip joint sub-region, and ΔFB represents the resolution reference, and DB i represents the contrast of the i-th hip joint sub-region, and ΔDB represents the contrast reference.

[0012] Optionally, the CT image analysis result of each hip joint sub-region is obtained by the processing, and corresponding image enhancement processing is performed based on the CT image analysis result, and the specific steps include: obtaining a preset CT image performance threshold in a database, and comparing the CT image performance threshold with the CT image performance index of each hip joint sub-region, if the CT image performance index of a certain hip joint sub-region is less than the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis unqualified, if the CT image performance index of a certain hip joint sub-region is above the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis qualified; if the CT image analysis result is CT image analysis qualified, the image enhancement processing is not performed, and if the CT image analysis result is CT image analysis unqualified, the image enhancement processing is performed.

[0013] Optionally, the imaging quality index of the overall CT image of the hip joint region is obtained by the analysis, and specifically includes: obtaining the CT image of each hip joint sub-region after the image enhancement processing, and splicing the CT image according to the order before the segmentation to obtain the overall CT image of the hip joint region; obtaining the imaging quality parameters of the overall CT image of the hip joint region, the imaging quality parameters of the overall CT image of the hip joint region including the artifact area ratio, the gray value range, the edge region gradient, the image noise and the image granularity; obtaining a preset imaging quality reference set in a database, and comparing the imaging quality reference set with the imaging quality parameters of the overall CT image of the hip joint region to obtain the imaging quality index of the overall CT image of the hip joint region; the imaging quality index of the overall CT image of the hip joint region is used to represent the imaging quality coincidence degree of the overall CT image of the hip joint region.

[0014] Optionally, the overall imaging verification result is determined according to the imaging quality index of the overall CT image of the hip joint region, and the specific steps include: obtaining a preset imaging quality threshold of the overall CT image of the hip joint region in a database, and comparing the imaging quality threshold with the imaging quality index of the overall CT image of the hip joint region; if the imaging quality index of the overall CT image of the hip joint region is above the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification qualified; and if the imaging quality index of the overall CT image of the hip joint region is less than the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification unqualified.

[0015] Optionally, the overall imaging verification result is determined according to the imaging quality index of the overall CT image of the hip joint region, and the specific steps include: obtaining a preset imaging quality threshold of the overall CT image of the hip joint region in a database, and comparing the imaging quality threshold with the imaging quality index of the overall CT image of the hip joint region; if the imaging quality index of the overall CT image of the hip joint region is above the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification qualified; and if the imaging quality index of the overall CT image of the hip joint region is less than the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification unqualified.

[0016] Optionally, the overall imaging verification result is determined according to the imaging quality index of the overall CT image of the hip joint region, and the specific steps include: obtaining a preset imaging quality threshold of the overall CT image of the hip joint region in a database, and comparing the imaging quality threshold with the imaging quality index of the overall CT image of the hip joint region; if the imaging quality index of the overall CT image of the hip joint region is above the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification qualified; and if the imaging quality index of the overall CT image of the hip joint region is less than the imaging quality threshold of the overall CT image of the hip joint region, the overall imaging verification result is overall imaging verification unqualified.

[0017] Optionally, the method further comprises obtaining the force parameters of the hip joint prosthesis in each simulation, and performing analysis to obtain a force index of the hip joint prosthesis, thereby outputting and displaying effective simulation information, and the specific steps include: obtaining the force parameters of the hip joint prosthesis in each simulation; obtaining a set of force demand factors of the hip joint prosthesis preset in the database, and comprehensively analyzing the force parameters of the hip joint prosthesis in each simulation to obtain a force index of the hip joint prosthesis in each simulation; obtaining a standard force index interval of the hip joint prosthesis preset in the database, and comparing the force index of the hip joint prosthesis in each simulation with the standard force index interval of the hip joint prosthesis, if the force index of the hip joint prosthesis in each simulation is within the standard force index interval of the hip joint prosthesis, the simulation of the hip joint prosthesis in each simulation is defined as effective simulation, and the effective simulation information is output and displayed.

[0018] Optionally, a hip joint prosthesis position selection device based on mechanical simulation, characterized in that it comprises: a CT image acquisition module, a CT image result analysis module, a whole imaging verification module, a hip joint prosthesis implantation simulation test execution judgment module and a result analysis module; the CT image acquisition module is used to acquire CT images after the mechanical simulation platform receives a CT upload signal, and segment the CT images to obtain hip joint region images, and segment the hip joint region images to obtain CT images of each hip joint sub-region; the CT image result analysis module is used to obtain CT image performance parameters of each hip joint sub-region for analysis to obtain CT image performance indexes of each hip joint sub-region, thereby processing to obtain CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; the whole imaging verification module is used to splice the CT images of each hip joint sub-region after image enhancement processing to obtain a whole hip joint region CT image, and analyze to obtain an imaging quality index of the whole hip joint region CT image, and judge a whole imaging verification result according to the imaging quality index of the whole hip joint region CT image, the whole imaging verification result including whole imaging verification qualification and whole imaging verification disqualification; the hip joint prosthesis implantation simulation test execution judgment module is used to judge based on the whole imaging verification result, if the whole imaging verification result is whole imaging verification qualification, the whole hip joint region CT image is input to a three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implantation simulation test; the result analysis module is used to process the hip joint prosthesis implantation simulation test data to obtain a corrected region image set, and process the corrected region image set to obtain a prosthesis position simulation analysis result of the whole hip joint region CT image.

[0019] The above technical solution has at least the following beneficial effects compared with the prior art:

[0020] 1. The hip joint prosthesis position selection method based on mechanical simulation provided by the application, through CT image segmentation based on the hip joint region image, and extraction of the CT image performance parameters of each hip joint sub-region for analysis and processing, the imaging quality index of the overall CT image of the hip joint region is obtained, and then the overall imaging verification result is obtained, effectively solving the problem of unreasonable prosthesis position selection caused by unqualified image quality processing in the prior art.

[0021] 2. The application discriminates based on the overall imaging verification result, if the overall imaging verification result is overall imaging verification qualified, the overall CT image of the hip joint region is input to the three-dimensional reconstruction module for hip joint prosthesis implant simulation test, so as to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region, and then the effective simulation information is output and displayed.

[0022] 3. The application obtains the deviation adjustment index of each implant sub-region by statistically analyzing the hip joint prosthesis implant simulation test data, combining the CT image performance index of each implant sub-region and the implant performance influence factor, and then obtains the corrected region image set and processes it. BRIEF DESCRIPTION OF DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the application, the following will briefly introduce the drawings needed to be used in the embodiment description. Obviously, the drawings in the following description are only some embodiments of the application, and other drawings can be obtained by those skilled in the art without creative labor.

[0024] Figure 1 The flowchart of the hip joint prosthesis position selection method based on mechanical simulation of the application;

[0025] Figure 2 The structural diagram of the hip joint prosthesis position selection device based on mechanical simulation of the application. DETAILED DESCRIPTION

[0026] In order to make the purpose, technical scheme and advantages of the embodiments of the application more clear, the technical scheme of the embodiments of the application will be described clearly and completely in combination with the drawings of the embodiments of the application. Obviously, the described embodiments are part of the embodiments of the application, not all the embodiments. Based on the described embodiments of the application, all other embodiments obtained by those skilled in the art without creative labor belong to the protection scope of the application.

[0027] As Figure 1As shown, a flow chart of a hip joint prosthesis position selection method provided by an embodiment of the application is provided, the method is applied to a hip joint prosthesis position selection device based on mechanical simulation, and the method comprises the following steps: S1, after the mechanical simulation platform receives a CT upload signal, CT images are acquired, the CT images are segmented to obtain a hip joint region image, and the hip joint region image is segmented to obtain CT images of each hip joint sub-region; S2, CT image performance parameters of each hip joint sub-region are acquired for analysis to obtain CT image performance indicators of each hip joint sub-region, and thus CT image analysis results of each hip joint sub-region are obtained, and corresponding image enhancement processing is performed based on the CT image analysis results; S3, the CT images of each hip joint sub-region after the image enhancement processing are acquired for splicing processing to obtain a whole hip joint region CT image, and imaging quality indicators of the whole hip joint region CT image are analyzed, and a whole imaging verification result is determined according to the imaging quality indicators of the whole hip joint region CT image, and the whole imaging verification result comprises whole imaging verification qualification and whole imaging verification disqualification;

[0028] S4, based on the whole imaging verification result, if the whole imaging verification result is whole imaging verification qualification, the whole hip joint region CT image is input to a three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implant simulation testing; S5, hip joint prosthesis implant simulation testing data are counted and processed to obtain a corrected region image set, and the corrected region image set is processed to obtain a prosthesis position simulation analysis result of the whole hip joint region CT image.

[0029] In this embodiment, it should be noted that the three-dimensional reconstruction module in the mechanical simulation platform is a component of the mechanical simulation platform, for example: VGSTUDIO MAX (Chinese name: VGSTUDIO MAX three-dimensional reconstruction and analysis software), Abaqus and CT2Model 3D plug-in (Chinese name: Abaqus finite element analysis software and CT2Model 3D three-dimensional reconstruction plug-in), Amira (Chinese name: Amira visualization and analysis software), etc.

[0030] Through segmentation and processing of the hip joint region image, key CT image performance parameters can be extracted, splicing of the whole CT image and determination of the imaging quality indicators provide a basis for accurate selection of the prosthesis position, the mechanical simulation platform can effectively analyze the stress state of the prosthesis under different positions, and through comprehensive analysis with a pre-set stress demand factor set of the hip joint prosthesis in the database, each simulation stress indicator of the hip joint prosthesis is obtained, and effective simulation information is further analyzed.

[0031] The CT image performance parameters of each hip joint sub-region are obtained, and the CT image performance indicators of each hip joint sub-region are obtained by analyzing the CT image performance parameters. The specific steps include: obtaining the CT image performance parameters of each hip joint sub-region, the CT image performance parameters of each hip joint sub-region including the average gray value, the average density, the noise value, the resolution and the contrast of each hip joint sub-region; obtaining the preset CT image reference set in the database, and comparing the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indicators of each hip joint sub-region; and the CT image performance indicators of each hip joint sub-region are used to represent the CT image quality compliance degree of each hip joint sub-region.

[0032] In the embodiment, the CT image performance parameters (average gray value, average density, noise value, resolution and contrast) of each hip joint sub-region can be obtained by statistical tools or programming libraries (such as MATLAB, Python and ImageJ / Fiji) in image processing software.

[0033] The CT image performance indicators of each hip joint sub-region are obtained, and the specific method is:

[0034]

[0035] In the formula, BX i represents the CT image performance indicators of the i-th hip joint sub-region, i represents the number of the hip joint sub-region, i = 1, 2,..., i max , i max represents the total number of hip joint sub-regions, e represents a natural constant, HD i represents the average gray value of the i-th hip joint sub-region, and ΔHD represents the average gray reference value, ρ i represents the average density of the i-th hip joint sub-region, and Δρ represents the average density reference value, ZS i represents the noise value of the i-th hip joint sub-region, and ΔZS represents the noise reference value, FB i represents the resolution of the i-th hip joint sub-region, and ΔFB represents the resolution reference value, DB i represents the contrast of the i-th hip joint sub-region, and ΔDB represents the contrast reference value.

[0036] In this embodiment, the CT image performance indicators of each hip joint sub-region are obtained by analyzing the CT image performance parameters of each hip joint sub-region. It is considered that there is a mutual influence relationship between these parameters. For example, a higher average gray value indicates a higher image brightness, and a lower contrast will make the gray value difference between different tissues not obvious, affecting the resolution of image content. However, if the gray value is too high or too low, it may lead to insufficient contrast and difficulty in distinguishing each region. If the average density is higher, it is helpful to distinguish bone and soft tissue and improve image clarity. When the noise value is high, the random interference in the image will cover or blur the boundaries between tissues, resulting in a decrease in contrast and making the gray difference between tissues not obvious. Higher image resolution can capture more details, but if the noise is large, the image details will be covered by the noise. Conversely, low resolution cannot clearly display all tissue details, and higher resolution can more clearly display tissue boundaries, thereby enhancing contrast. If the contrast is too low, even if the resolution is high, it is still difficult to distinguish different tissues in the image.

[0037] By obtaining and analyzing the CT image performance parameters of each hip joint sub-region, including average gray value, average density, noise value, resolution and contrast, and comparing with the preset CT image reference set in the database, the image quality of each sub-region is evaluated. The image quality of each sub-region is quantified to ensure the CT image compliance accuracy and avoid the problem of inaccurate subsequent hip joint prosthesis analysis caused by unqualified image quality. By discovering and processing image problems in time before hip joint prosthesis implantation, the accuracy and reliability of prosthesis position selection are improved.

[0038] The CT image analysis results of each hip joint sub-region are obtained by the processing, and corresponding image enhancement processing is performed based on the CT image analysis results. The specific steps include: obtaining the preset CT image performance threshold in the database, and comparing with the CT image performance indicators of each hip joint sub-region. If the CT image performance indicator of a certain hip joint sub-region is less than the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis unqualified. If the CT image performance indicator of a certain hip joint sub-region is above the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis qualified. If the CT image analysis result is CT image analysis qualified, no image enhancement processing is performed. If the CT image analysis result is CT image analysis unqualified, image enhancement processing is performed.

[0039] In this embodiment, it should be noted that if the CT image analysis result is CT image analysis unqualified, image enhancement processing is performed, and the CT image is processed by adaptive histogram equalization for image enhancement processing. The processing can be realized by image processing software (such as Photoshop or professional image processing library such as OpenCV).

[0040] By performing image enhancement processing based on the CT image analysis result in the application, the accuracy of hip prosthesis position selection and the image quality can be effectively improved. By comparing the CT image performance indicators of each hip joint sub-region with the preset threshold, it can be quickly identified whether the CT image of each hip joint sub-region meets the quality standard, avoiding the waste of computing resources caused by using unqualified images for prosthesis position selection. For unqualified images, the image enhancement processing improves the clarity, contrast and details of the image, effectively reducing the influence of image quality problems on the accuracy of prosthesis positioning. By automatically determining whether image enhancement processing is needed, the complexity and error of manual intervention are reduced, the operation efficiency is improved, and the high precision and high reliability of hip prosthesis position selection are ensured.

[0041] In the embodiment, it is necessary to explain that the imaging quality reference set includes the artifact area ratio reference value, the gray value range reference value, the edge region gradient reference value, the image noise reference value and the image granularity reference value.

[0042] In the embodiment, it is necessary to explain that the imaging quality reference set includes the artifact area ratio reference value, the gray value range reference value, the edge region gradient reference value, the image noise reference value and the image granularity reference value.

[0043] The imaging quality index of the overall CT image of the hip joint region is obtained, and the specific method is as follows:

[0044]

[0045] In the formula, ZL represents the imaging quality index of the overall CT image of the hip joint region, WY represents the artifact area ratio, ΔWY represents the artifact area ratio reference value, HC represents the gray value range, ΔHC represents the gray value range reference value, TD represents the edge region gradient, ΔTD represents the edge region gradient reference value, TZ represents the image noise, ΔTZ represents the image noise reference value, KL represents the image granularity, ΔKL represents the image granularity reference value, and e represents the natural constant.

[0046] It should be noted that the artifact area ratio refers to the proportion of the artifact area to the total area, and the artifact is usually caused by device error or irregular noise in the imaging process. By analyzing the imaging quality parameters (artifact area ratio, gray value range, edge region gradient, image noise and image granularity) of the overall CT image of the hip joint region, the imaging quality index of the overall CT image of the hip joint region is obtained, considering the mutual influence relationship between these parameters, for example: the greater the image noise, the more the appearance of artifacts will increase, resulting in an increase in the artifact area ratio, the gray value range is the gray difference between the brightest and darkest areas in the image, the greater the gray difference, the higher the gradient value of the edge region, which is conducive to improving the clarity and tissue differentiation of the image. The edge region gradient represents the distinctness of the structure in the image, the higher the image granularity, the rougher the image, resulting in a decrease in the gradient of the edge region, affecting the selection of the prosthesis position and the accuracy of the image. Lower image granularity helps to clearly present the edge details, and larger image granularity is prone to noise and artifacts, thereby reducing image quality, and an increase in image noise will affect the stability of the gray value, resulting in uneven changes in the gray value.

[0047] In the embodiment, by comparing the imaging quality index of the overall CT image of the hip joint region with the preset quality threshold, it can be effectively judged whether the image quality meets the standard, if the image quality meets the standard, the further mechanical simulation can be carried out on the basis of the high-quality image, and the accuracy and reliability of the simulation result are improved. If the image quality does not meet the standard, the corresponding image processing can be identified and made in time, which avoids the influence of image quality problems on the accuracy of the selection of the prosthesis position, effectively improves the automation level of image quality control, reduces manual intervention, and ensures the efficiency and accuracy of subsequent processing.

[0048] In the embodiment, by comparing the imaging quality index of the overall CT image of the hip joint region with the preset quality threshold, it can be effectively judged whether the image quality meets the standard, if the image quality meets the standard, the further mechanical simulation can be carried out on the basis of the high-quality image, and the accuracy and reliability of the simulation result are improved. If the image quality does not meet the standard, the corresponding image processing can be identified and made in time, which avoids the influence of image quality problems on the accuracy of the selection of the prosthesis position, effectively improves the automation level of image quality control, reduces manual intervention, and ensures the efficiency and accuracy of subsequent processing.

[0049] The overall imaging verification result is determined based on the overall imaging verification, and specifically includes: if the overall imaging verification result is overall imaging verification qualified, the hip joint region overall CT image is input to the three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implantation simulation test, if the overall imaging verification result is overall imaging verification unqualified, the CT image is reacquired, and S1 and S2 are repeatedly executed within a preset number of times; if the overall imaging verification result after repeatedly executing S1 and S2 within the preset number of times is overall imaging verification qualified, the repeated execution is stopped, and if the overall imaging verification result after repeatedly executing S1 and S2 within the preset number of times is overall imaging verification unqualified, a warning prompt information is sent.

[0050] In the embodiment, when the image quality is determined to be qualified or not through the overall imaging verification result, if the image quality is qualified, the hip joint region overall CT image can be directly input to the three-dimensional reconstruction module in the mechanical simulation platform for simulation test of hip joint prosthesis implantation. Through this step, the rationality of prosthesis implantation can be analyzed in depth, and the rationality of prosthesis implantation is improved. If the image quality fails to pass the overall imaging verification, the step of reacquiring the CT image and performing image processing is automatically triggered. This process will repeatedly execute the image acquisition and analysis steps (S1 and S2) according to a preset number of times, to ensure that qualified image data is obtained after multiple attempts. If the image quality still fails to reach the standard after multiple attempts, a warning prompt information is sent, which includes warning reasons and warning time, etc.

[0051] The prosthesis position simulation analysis result of the hip joint region overall CT image is obtained after processing the correction region image set, and the specific steps include: acquiring the hip joint region overall CT image after hip joint prosthesis implantation, segmenting to obtain the CT image of each hip joint sub-region after hip joint prosthesis implantation, and marking as the CT image of each implantation sub-region, acquiring the CT image performance parameters of each implantation sub-region, and analyzing to obtain the CT image performance indicators of each implantation sub-region; acquiring the preset implantation performance influence factor in the database, and comprehensively analyzing based on the implantation performance influence factor, the CT image performance indicators of each hip joint sub-region and the CT image performance indicators of each implantation sub-region to obtain the deviation adjustment indicators of each implantation sub-region; acquiring the preset deviation adjustment threshold in the database, and comparing with the deviation adjustment indicators of each implantation sub-region, if the deviation adjustment indicators of a certain implantation sub-region are above the deviation adjustment threshold, the implantation sub-region is marked as the correction region image set; after processing the correction region image set, if the deviation adjustment indicators of the correction region image set are above the deviation adjustment threshold, the prosthesis position simulation analysis result of the hip joint region overall CT image is divided into simulation verification unqualified, and if the deviation adjustment indicators of the correction region image set are less than the deviation adjustment threshold, the prosthesis position simulation analysis result of the hip joint region overall CT image is divided into simulation verification qualified.

[0052] In the present embodiment, it is to be noted that the prosthesis position simulation analysis result includes simulation verification pass and simulation verification fail. The revision region image set includes each implant sub-region image corresponding to the deviation adjustment index of the implant sub-region being above the deviation adjustment threshold.

[0053] It is to be noted that if the prosthesis position simulation analysis result of the overall CT image of the hip joint region is divided into simulation verification pass, the next step of stress simulation of the hip joint prosthesis is performed, and if the prosthesis position simulation analysis result of the overall CT image of the hip joint region is divided into simulation verification fail, the stress simulation is not performed and a warning prompt information is issued, the warning prompt information includes issuing a prompt of simulation verification fail and sounding an alarm bell.

[0054] It is to be noted that each hip joint sub-region and each implant sub-region are in a one-to-one correspondence, that is, the hip joint sub-region and the implant sub-region are in the same position, and only the difference before and after the hip joint prosthesis implantation exists.

[0055] It is to be further noted that the CT image performance parameter of each implant sub-region can be processed by the CT image performance index of each hip joint sub-region to obtain the CT image performance index of each implant sub-region.

[0056] The deviation adjustment index of each implant sub-region is obtained, and the specific method is as follows:

[0057]

[0058] In the formula, PT j represents the deviation adjustment index of the jth implant sub-region, j represents the number of implant sub-regions, j = 1, 2,..., j max , j max represents the total number of implant sub-regions, BX i represents the CT image performance index of the ith hip joint sub-region, i represents the number of hip joint sub-regions, i = 1, 2,..., i max , i max represents the total number of hip joint sub-regions, RBX j represents the CT image performance index of the jth implant sub-region, e represents a natural constant, Y z represents the implant performance influence factor.

[0059] The application can effectively ensure the accuracy and stability of the prosthesis position by segmenting, analyzing and correcting the region image set of the overall CT image after the hip joint prosthesis implantation. After obtaining the post-implantation CT image and analyzing each implant sub-region, the implant performance influencing factors preset in the database are combined to comprehensively evaluate the adaptation of the prosthesis and the hip joint region image quality, and the deviation adjustment index is optimized according to the analysis result to ensure that the image quality of each sub-region is optimized. The implant performance influencing factors preset in the database are used to represent the influence degree after the hip joint prosthesis implantation. By comparing the deviation adjustment index with the preset threshold value, the regions that need to be corrected can be automatically identified and marked. The automatic correction and verification mechanism reduces manual intervention and improves the efficiency and accuracy of the operation.

[0060] The stress parameters of the hip joint prosthesis in each simulation are obtained, and the stress indicators of the hip joint prosthesis are obtained by analysis, so that the effective simulation information is output and displayed. The specific steps include: obtaining the stress parameters of the hip joint prosthesis in each simulation; obtaining the preset stress demand factor set of the hip joint prosthesis in the database and comprehensively analyzing the stress parameters of the hip joint prosthesis in each simulation to obtain the stress indicators of the hip joint prosthesis in each simulation; obtaining the standard hip joint prosthesis stress indicator interval preset in the database and comparing it with the stress indicators of the hip joint prosthesis in each simulation. If the stress indicator of the hip joint prosthesis in a certain simulation is within the standard hip joint prosthesis stress indicator interval, the simulation of the hip joint prosthesis is defined as effective simulation, and the effective simulation information is output and displayed. If the stress indicator of the hip joint prosthesis in a certain simulation is not within the standard hip joint prosthesis stress indicator interval, the simulation of the hip joint prosthesis is defined as invalid simulation and is not displayed.

[0061] In the embodiment, the stress parameters of the hip joint prosthesis include the prosthesis and acetabulum contact surface friction force reference value, acetabulum and prosthesis contact stress, prosthesis and bone tissue shear stress, prosthesis and acetabulum normal force, hip joint flexion and extension torque, hip joint abduction and adduction torque, and hip joint internal and external rotation torque.

[0062] The stress demand factor set of the hip joint prosthesis includes: unit acetabulum and prosthesis contact stress demand factor, unit prosthesis and bone tissue shear stress demand factor, unit prosthesis and acetabulum normal force demand factor, unit hip joint abduction and adduction torque demand factor, unit hip joint internal and external rotation torque demand factor, prosthesis and acetabulum contact surface friction force reference value, and hip joint flexion and extension torque reference value.

[0063] The contact stress refers to the pressure applied between the prosthesis and the acetabulum contact surface, the shear stress refers to the transverse stress generated at the contact surface between the prosthesis and the bone tissue due to movement and force, the normal force refers to the force perpendicular to the contact surface between the prosthesis and the acetabulum, indicating the pressure distribution between the prosthesis and the acetabulum, the friction force refers to the resistance generated between the contact surface between the prosthesis and the acetabulum due to contact and movement, the flexion-extension torque refers to the torque received by the hip joint during the flexion-extension movement, generated by the contact relationship and movement between the femoral head and the acetabulum, the abduction-adduction torque refers to the torque applied by the hip joint during the abduction or adduction movement, the internal-external rotation torque refers to the torque applied by the hip joint during the internal-external rotation, and is located between the femoral head and the acetabulum.

[0064] The stress parameters of the hip joint prosthesis in each simulation are obtained by analyzing the stress parameters of the hip joint prosthesis in each simulation, considering the mutual influence relationship between these parameters, for example: the greater the normal force, the greater the contact area between the prosthesis and the acetabulum, thereby increasing the contact stress, and when the normal force increases, the influence of the shear stress will be intensified. Excessive shear stress will cause separation between the prosthesis and the bone tissue, affecting the contact stability. When the normal force is large, the effect of the shear stress is increased, thereby causing the displacement of the prosthesis and the risk of fracture. Moderate friction force can provide stable joint movement, but excessive friction force will cause wear, and too small friction force will cause the prosthesis to be unstable, and a larger flexion-extension torque will increase the torsional torque, thereby affecting the stability of the prosthesis, causing the rotation or displacement of the prosthesis. Excessive abduction-adduction torque will produce excessive rotational torque, causing the displacement or dislocation of the prosthesis, and the contact stress and the shear stress are usually in the same direction or interact with each other. Greater contact stress will intensify the influence of shear stress, and when the contact stress increases, the shear stress will cause damage to the bone tissue or dislocation of the prosthesis, affecting the stability of the joint.

[0065] The stress parameters of the hip joint prosthesis in each simulation are obtained by analyzing the stress parameters of the hip joint prosthesis in each simulation, considering the mutual influence relationship between these parameters, for example: the greater the normal force, the greater the contact area between the prosthesis and the acetabulum, thereby increasing the contact stress, and when the normal force increases, the influence of the shear stress will be intensified. Excessive shear stress will cause separation between the prosthesis and the bone tissue, affecting the contact stability. When the normal force is large, the effect of the shear stress is increased, thereby causing the displacement of the prosthesis and the risk of fracture. Moderate friction force can provide stable joint movement, but excessive friction force will cause wear, and too small friction force will cause the prosthesis to be unstable, and a larger flexion-extension torque will increase the torsional torque, thereby affecting the stability of the prosthesis, causing the rotation or displacement of the prosthesis. Excessive abduction-adduction torque will produce excessive rotational torque, causing the displacement or dislocation of the prosthesis, and the contact stress and the shear stress are usually in the same direction or interact with each other. Greater contact stress will intensify the influence of shear stress, and when the contact stress increases, the shear stress will cause damage to the bone tissue or dislocation of the prosthesis, affecting the stability of the joint.

[0066]

[0067] In the formula, L n represents the stress index of the hip joint prosthesis in the nth simulation, n represents the number of simulation times, n = 1, 2,..., n max , n max represents the total number of simulations, e represents the natural constant, KY u represents the contact stress between the acetabulum and the prosthesis of the hip joint prosthesis in the nth simulation, ΔKY represents the unit acetabulum and prosthesis contact stress requirement factor, KQ n represents the shear stress between the prosthesis and the bone tissue of the hip joint prosthesis in the nth simulation, ΔKQ represents the unit shear stress requirement factor between the prosthesis and the bone tissue of the hip joint prosthesis, KF nΔKF represents the unit prosthesis-acetabulum normal force demand factor, KM n ΔKM represents the prosthesis-acetabulum contact surface friction reference value, KN n ΔKN represents the hip flexion and extension torque reference value, KS n ΔKS represents the unit hip abduction and adduction moment demand factor, KX n ΔKX represents the unit hip internal and external rotation moment demand factor.

[0068] It should be noted that the effective simulation information includes the force condition of the hip joint prosthesis, the installation position of the hip joint prosthesis, and the installation angle of the hip joint prosthesis, etc. The force index of the hip joint prosthesis in each simulation is used to represent the force degree of the hip joint prosthesis in the force process of each simulation.

[0069] The processing of the corrected region image set can be processed by using an image registration and fusion algorithm, and the specific method is to register and fuse the image of the corrected region and the CT image of any hip joint sub-region whose CT image analysis result is qualified for CT image analysis, and the registration method based on feature points such as SIFT (full name: Scale Invariant Feature Transform), ORB (full name: Oriented FAST and Rotated BRIEF) processing is adopted.

[0070] As shown in Figure 2 Fig. 1 is a structure diagram of a hip joint prosthesis position selection device based on mechanical simulation provided by an embodiment of the present application. The hip joint prosthesis position selection device based on mechanical simulation provided by the embodiment of the present application comprises:

[0071] The CT image acquisition module, the CT image result analysis module, the overall imaging verification module, the hip joint prosthesis implantation simulation test execution judgment module and the result analysis module; the CT image acquisition module is used for acquiring the CT image after the mechanical simulation platform receives the CT upload signal, and segmenting the CT image to obtain the hip joint region image, and segmenting the hip joint region image to obtain the CT image of each hip joint sub-region; the CT image result analysis module is used for acquiring and analyzing the CT image performance parameters of each hip joint sub-region to obtain the CT image performance index of each hip joint sub-region, thereby processing to obtain the CT image analysis result of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis result; the overall imaging verification module is used for acquiring the CT image of each hip joint sub-region after image enhancement processing to perform splicing processing, obtaining the overall CT image of the hip joint region, and analyzing to obtain the imaging quality index of the overall CT image of the hip joint region, judging the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint region, and the overall imaging verification result includes overall imaging verification qualified and overall imaging verification unqualified; the hip joint prosthesis implantation simulation test execution judgment module is used for judging based on the overall imaging verification result, if the overall imaging verification result is overall imaging verification qualified, the overall CT image of the hip joint region is input to the three-dimensional reconstruction module in the mechanical simulation platform to perform hip joint prosthesis implantation simulation test; the result analysis module is used for processing the hip joint prosthesis implantation simulation test data to obtain the corrected region image set, and processing the corrected region image set to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region.

[0072] The above is only a specific embodiment of the present application, but the protection scope of the present application is not limited thereto, and the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for selecting a position of a hip joint prosthesis based on a mechanical simulation, characterized in that, The method comprises the following steps: S1, acquiring a CT image when a mechanical simulation platform receives a CT upload signal, segmenting the CT image to obtain a hip joint region image, and segmenting the hip joint region image to obtain a CT image of each hip joint sub-region; S2, acquiring CT image performance parameters of each hip joint sub-region to analyze and obtain CT image performance indicators of each hip joint sub-region, thereby processing to obtain CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; The CT image performance parameters of each hip joint sub-region are acquired, and the CT image performance parameters of each hip joint sub-region include average gray value, average density, noise value, resolution and contrast of each hip joint sub-region; a preset CT image reference set in a database is acquired and compared with the CT image performance parameters of each hip joint sub-region to obtain CT image performance indicators of each hip joint sub-region; the CT image performance indicators of each hip joint sub-region are used to represent the CT image quality compliance degree of each hip joint sub-region; The CT image performance indicators of each hip joint sub-region are obtained by the following method: S3, acquiring the CT images of each hip joint sub-region after image enhancement processing to perform splicing processing, obtaining a whole hip joint region CT image, analyzing to obtain a rendering quality indicator of the whole hip joint region CT image, and judging a whole rendering verification result according to the rendering quality indicator of the whole hip joint region CT image; the whole rendering verification result includes whole rendering verification qualified and whole rendering verification unqualified; In the formula, BX i represents the CT image performance index of the i-th hip joint sub-region, i represents the number of the hip joint sub-region, i = 1, 2, …, i max , i max represents the total number of hip joint sub-regions, e represents a natural constant, HD i represents the average gray value of the i-th hip joint sub-region, ΔHD represents the average gray reference value, ρ i represents the average density of the i-th hip joint sub-region, Δρ represents the average density reference value, ZS i represents the noise value of the i-th hip joint sub-region, ΔZS represents the noise reference value, FB i represents the resolution of the i-th hip joint sub-region, ΔFB represents the resolution reference value, DB i represents the contrast of the i-th hip joint sub-region, and ΔDB represents the contrast reference value. S4, discriminating based on the whole rendering verification result; if the whole rendering verification result is whole rendering verification qualified, the whole hip joint region CT image is input into a three-dimensional reconstruction module in the mechanical simulation platform to perform hip joint prosthesis implantation simulation inspection; S5, processing the hip joint prosthesis implantation simulation inspection data to obtain a corrected region image set, and processing the corrected region image set to obtain a prosthesis position simulation analysis result of the whole hip joint region CT image. The CT image analysis results of each hip joint sub-region are obtained by the following method:

2. The mechanics simulation based hip prosthesis position selection method according to claim 1, characterized in that, The CT image performance threshold in the database is acquired and compared with the CT image performance indicators of each hip joint sub-region; if the CT image performance indicator of a certain hip joint sub-region is less than the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis unqualified; if the CT image performance indicator of a certain hip joint sub-region is greater than or equal to the CT image performance threshold, the CT image analysis result of the hip joint sub-region is CT image analysis qualified; If the CT image analysis result is CT image analysis qualified, no image enhancement processing is performed; if the CT image analysis result is CT image analysis unqualified, image enhancement processing is performed. The rendering quality indicator of the whole hip joint region CT image is obtained by the following method:

3. The mechanics simulation based hip prosthesis position selection method of claim 1, wherein, ​ The CT images of each hip joint sub-region after image enhancement processing are obtained, and the CT images are spliced in the order before segmentation to obtain a whole hip joint region CT image; An imaging quality parameter of the whole hip joint region CT image is obtained, and the imaging quality parameter of the whole hip joint region CT image includes an artifact area ratio, a gray value range, an edge region gradient, an image noise, and an image granularity; A preset imaging quality reference set in the database is obtained, and compared with the imaging quality parameter of the whole hip joint region CT image to obtain an imaging quality index of the whole hip joint region CT image. The imaging quality index of the whole hip joint region CT image is used to represent the imaging quality of the whole hip joint region CT image.

4. The mechanics simulation based hip prosthesis position selection method of claim 1, wherein, The imaging quality index of the whole hip joint region CT image is used to represent the imaging quality of the whole hip joint region CT image. The imaging quality threshold of the whole hip joint region CT image in the database is obtained, and compared with the imaging quality index of the whole hip joint region CT image. If the imaging quality index of the whole hip joint region CT image is above the imaging quality threshold of the whole hip joint region CT image, the whole imaging verification result is whole imaging verification qualified, and if the imaging quality index of the whole hip joint region CT image is less than the imaging quality threshold of the whole hip joint region CT image, the whole imaging verification result is whole imaging verification unqualified.

5. The mechanics simulation based hip prosthesis position selection method of claim 1, wherein, The whole imaging verification result is used for judgment, and the specific steps include: If the whole imaging verification result is whole imaging verification qualified, the whole hip joint region CT image is input into a three-dimensional reconstruction module in a mechanical simulation platform for hip joint prosthesis implant simulation test. If the whole imaging verification result is whole imaging verification unqualified, the CT image is reacquired, and S1 and S2 are repeatedly executed within a preset number of times. If the whole imaging verification result after repeatedly executing S1 and S2 within a preset number of times is whole imaging verification qualified, the repeated execution is stopped. If the whole imaging verification result after repeatedly executing S1 and S2 within a preset number of times is whole imaging verification unqualified, a warning prompt information is sent.

6. The mechanics simulation based hip prosthesis positioning method according to claim 1, wherein, The specific steps of obtaining the prosthesis position simulation analysis result of the whole hip joint region CT image after processing the correction region image set include: A whole hip joint region CT image after hip joint prosthesis implantation is obtained, segmented to obtain CT images of each hip joint sub-region after hip joint prosthesis implantation, and marked as CT images of each implant sub-region. CT image performance parameters of each implant sub-region are obtained, analyzed to obtain CT image performance indexes of each implant sub-region; Preset implant performance influence factors in the database are obtained, and comprehensive analysis is performed based on the implant performance influence factors, CT image performance indexes of each hip joint sub-region, and CT image performance indexes of each implant sub-region to obtain deviation adjustment indexes of each implant sub-region; Preset deviation adjustment thresholds in the database are obtained, and compared with the deviation adjustment indexes of each implant sub-region. If the deviation adjustment index of a certain implant sub-region is above the deviation adjustment threshold, the implant sub-region is marked as a correction region image set. After processing the correction region image set, if the deviation adjustment index of the correction region image set is above the deviation adjustment threshold, the prosthesis position simulation analysis result of the hip joint region overall CT image is divided into simulation verification unqualified; if the deviation adjustment index of the correction region image set is less than the deviation adjustment threshold, the prosthesis position simulation analysis result of the hip joint region overall CT image is divided into simulation verification qualified.

7. The mechanics simulation based hip prosthesis positioning method according to claim 1, wherein, Further comprising obtaining the stress parameters of the hip joint prosthesis in each simulation, and analyzing to obtain the stress index of the hip joint prosthesis, thereby outputting and displaying the effective simulation information, the specific steps comprising: obtaining the stress parameters of the hip joint prosthesis in each simulation; obtaining the stress demand factor set of the hip joint prosthesis preset in the database, and comprehensively analyzing the stress parameters of the hip joint prosthesis in each simulation to obtain the stress index of the hip joint prosthesis in each simulation; obtaining the standard hip joint prosthesis stress index interval preset in the database, and comparing with the stress index of the hip joint prosthesis in each simulation, if the stress index of the hip joint prosthesis in each simulation is within the standard hip joint prosthesis stress index interval, the simulation of the hip joint prosthesis is defined as effective simulation, and the effective simulation information is output and displayed.

8. A device for mechanical simulation based hip prosthesis position selection, the device being for implementing the method according to any one of claims 1 to 7, characterized in that, The device comprises a CT image acquisition module, a CT image result analysis module, an overall imaging verification module, a hip joint prosthesis implantation simulation inspection execution judgment module and a result analysis module; The CT image acquisition module is used to acquire the CT image after the mechanical simulation platform receives the CT upload signal, and segment the CT image to obtain the hip joint region image, and segment the hip joint region image to obtain the CT image of each hip joint sub-region; The CT image result analysis module is used to obtain the CT image performance parameters of each hip joint sub-region for analysis to obtain the CT image performance index of each hip joint sub-region, thereby processing to obtain the CT image analysis result of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis result; The overall imaging verification module is used to splice the CT images of each hip joint sub-region after image enhancement processing to obtain the hip joint region overall CT image, and analyze to obtain the imaging quality index of the hip joint region overall CT image, and judge the overall imaging verification result according to the imaging quality index of the hip joint region overall CT image, the overall imaging verification result comprising overall imaging verification qualified and overall imaging verification unqualified; The hip joint prosthesis implantation simulation inspection execution judgment module is used to judge based on the overall imaging verification result, if the overall imaging verification result is overall imaging verification qualified, the hip joint region overall CT image is input to the three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implantation simulation inspection; The result analysis module is used to process the hip joint prosthesis implantation simulation inspection data to obtain the correction region image set, and process the correction region image set to obtain the prosthesis position simulation analysis result of the hip joint region overall CT image.

Images