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

Abstract

The invention provides a hip joint prosthesis position selection method and device based on mechanical simulation, and belongs to the technical field of prostheses, and the method comprises the steps: obtaining a CT image, and obtaining the CT image of each hip joint sub-region; the CT image analysis result of each hip joint sub-region is obtained through processing, and corresponding image enhancement processing is carried out; an integral CT image of the hip joint area is obtained, an imaging quality index of the integral CT image of the hip joint area is obtained through analysis, and an integral imaging verification result is judged; judging based on the overall imaging verification result, and if the overall imaging verification result shows that the overall imaging verification is qualified, inputting the overall CT image of the hip joint area into a three-dimensional reconstruction module in a mechanical simulation platform to perform hip joint prosthesis implantation simulation test; and counting hip joint prosthesis implantation simulation test data, processing the 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 integral CT image of the hip joint region.

Description

Technical Field

[0001] The present invention relates to the technical field of prostheses, and in particular to a method and device for selecting the position of a hip joint prosthesis based on mechanical simulation. Background Art

[0002] With the continuous progress of technology, significant development has been achieved in hip joint prosthesis technology, which is widely used in clinical practice. The placement position and angle of hip joint prostheses are crucial for the mechanical balance and recovery of motor function of postoperative patients. Therefore, targeted position selection needs to be carried out in the early stage of hip joint prosthesis placement.

[0003] For example, a control method for an arthroscopic surgery robot disclosed in the invention patent announcement with the publication number CN114305697B includes: taking images of the patient's bones and transmitting the images to a computer system, and the computer system obtains the image data of the bones; 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 between the joint prosthesis model and the bones is the part that needs to be resected and replaced; installing a signal source as a signal sending end, taking images of the patient's bones again, with the shooting range including the signal source and the joint part. In the computer system, a coordinate system is established with the position of the signal source as the origin. The computer system automatically calculates and memorizes the interface of the overlapping part between the joint prosthesis model and the bones, and its coordinates in the coordinate system. The surgical robot receives the coordinates and combines the signals sent by the received signal source to obtain the position coordinates of the part that needs to be resected from the diseased joint in the coordinate system, and the surgical robot performs surgical operations.

[0004] In the prior art, by taking images of the patient's bones and identifying and matching the diseased joints, and directly using the bone images for position selection. However, in actual application, there are often differences in the quality of bone images, and the quality of bone images will seriously affect the accuracy of positioning. Therefore, there is a problem in the prior art that the prosthesis position selection is unreasonable due to unqualified image quality processing. Summary of the Invention

[0005] Embodiments of the present invention provide a method and device for selecting the position of a hip joint prosthesis based on mechanical simulation, which solve the problem of unreasonable prosthesis position selection caused by unqualified image quality processing in the prior art, and achieve the accuracy and controllability of prosthesis position selection.

[0006] To achieve the above-mentioned invention purpose, the technical solutions provided by the present invention are as follows:

[0007] An embodiment of the present invention provides a method for selecting the position of a hip joint prosthesis based on mechanical simulation, including the following steps: S1. After the mechanical simulation platform receives the CT upload signal, obtain the CT image, segment the CT image to obtain the hip joint region image, and segment the hip joint region image to obtain the CT images of each hip joint sub-region; S2. Obtain the CT image performance parameters of each hip joint sub-region, analyze them to obtain the CT image performance indicators of each hip joint sub-region, thereby process to obtain the CT image analysis results of each hip joint sub-region, and perform corresponding image enhancement processing based on the CT image analysis results; S3. Obtain the CT images of each hip joint sub-region after image enhancement processing, perform stitching processing to obtain the overall CT image of the hip joint region, analyze to obtain the imaging quality index of the overall CT image of the hip joint region, and judge the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint region. The overall imaging verification result includes passing the overall imaging verification and failing the overall imaging verification; S4. Based on the overall imaging verification result, make a discrimination. If the overall imaging verification result is passing the overall imaging verification, input the overall CT image of the hip joint region into the three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implantation simulation test; S5. Statistically process the hip joint prosthesis implantation simulation test data to obtain a set of corrected region images, and process the set of corrected region images to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region.

[0008] Optionally, the specific steps of obtaining the CT image performance parameters of each hip joint sub-region, analyzing them to obtain the CT image performance indicators of each hip joint sub-region include: obtaining the CT image performance parameters of each hip joint sub-region, where the CT image performance parameters of each hip joint sub-region include the average gray value, average density, noise value, resolution, and contrast of each hip joint sub-region; obtaining the preset CT image reference set in the database, and comparing it with the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indicators of each hip joint sub-region; the CT image performance indicators of each hip joint sub-region are used to characterize the compliance degree of the CT image quality of each hip joint sub-region.

[0009] Optionally, the specific method for obtaining the CT image performance indicators of each hip joint sub-region is:

[0010]

[0011] In the formula, BX i represents the CT image performance indicator 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 the natural constant, HD irepresents the average gray value of the i-th hip sub-region, ΔHD represents the average gray reference value, ρ i represents the average density of the i-th hip sub-region, Δρ represents the average density reference value, ZS i represents the noise value of the i-th hip sub-region, ΔZS represents the noise reference value, FB i represents the resolution of the i-th hip sub-region, ΔFB represents the resolution reference value, DB i represents the contrast of the i-th hip sub-region, ΔDB represents the contrast reference value.

[0012] Optionally, analyze the CT image analysis results of each hip sub-region obtained thereby, and perform corresponding image enhancement processing based on the CT image analysis results. The specific steps include: obtaining the preset CT image performance threshold in the database, and comparing it with the CT image performance indicators of each hip sub-region. If the CT image performance indicator of a certain hip sub-region is less than the CT image performance threshold, the CT image analysis result of this hip sub-region is unqualified for CT image analysis. If the CT image performance indicator of a certain hip sub-region is above the CT image performance threshold, the CT image analysis result of this hip sub-region is qualified for CT image analysis; if the CT image analysis result is qualified for CT image analysis, do not perform image enhancement processing. If the CT image analysis result is unqualified for CT image analysis, perform image enhancement processing.

[0013] Optionally, obtain the imaging quality indicators of the overall CT image of the hip region, specifically including: obtaining the CT images of each hip sub-region after image enhancement processing, and splicing them in the order before segmentation to obtain the overall CT image of the hip region; obtaining the imaging quality parameters of the overall CT image of the hip region, and the imaging quality parameters of the overall CT image of the hip region include the artifact area ratio, gray value range, edge region gradient, image noise, and image granularity; obtaining the preset imaging quality reference set in the database, and comparing it with the imaging quality parameters of the overall CT image of the hip region to obtain the imaging quality indicators of the overall CT image of the hip region; the imaging quality indicators of the overall CT image of the hip region are used to characterize the compliance degree of the imaging quality of the overall CT image of the hip region.

[0014] Optionally, judging the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint area, the specific steps include: obtaining the imaging quality threshold of the overall CT image of the hip joint area preset in the database, and comparing it with the imaging quality index of the overall CT image of the hip joint area. If the imaging quality index of the overall CT image of the hip joint area is above the imaging quality threshold of the overall CT image of the hip joint area, the overall imaging verification result is that the overall imaging verification is qualified. If the imaging quality index of the overall CT image of the hip joint area is less than the imaging quality threshold of the overall CT image of the hip joint area, the overall imaging verification result is that the overall imaging verification is unqualified.

[0015] Optionally, the discrimination based on the overall imaging verification result specifically includes: if the overall imaging verification result is that the overall imaging verification is qualified, input the overall CT image of the hip joint area into the three-dimensional reconstruction module in the mechanical simulation platform for hip prosthesis implantation simulation inspection. If the overall imaging verification result is that the overall imaging verification is unqualified, re-obtain the CT image and repeat S1 and S2 within the preset number of times; if the overall imaging verification result after repeating S1 and S2 within the preset number of times is that the overall imaging verification is qualified, stop repeating the operation. If the overall imaging verification result after repeating S1 and S2 within the preset number of times is that the overall imaging verification is unqualified, issue a warning prompt message.

[0016] Optionally, the specific steps for obtaining the prosthesis position simulation analysis result of the overall CT image of the hip joint area after processing the corrected region image set include: obtaining the overall CT image of the hip joint area after hip prosthesis implantation, performing segmentation to obtain the CT images of each hip joint sub-region after hip prosthesis implantation, and marking them as the CT images of each implanted sub-region. Obtain the CT image performance parameters of each implanted sub-region and analyze them to obtain the CT image performance indicators of each implanted sub-region; obtain the implanted performance influencing factors preset in the database, and based on the implanted performance influencing factors, the CT image performance indicators of each hip joint sub-region, and the CT image performance indicators of each implanted sub-region, conduct comprehensive analysis to obtain the deviation adjustment indicators of each implanted sub-region; obtain the deviation adjustment threshold preset in the database and compare it with the deviation adjustment indicators of each implanted sub-region. If the deviation adjustment indicator of a certain implanted sub-region is above the deviation adjustment threshold, mark the implanted sub-region as the corrected region image set; after processing the corrected region image set, if the deviation adjustment indicator of the corrected region image set is above the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as unqualified for simulation verification. If the deviation adjustment indicator of the corrected region image set is less than the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as qualified for simulation verification.

[0017] Optionally, it further includes obtaining the force parameters of the hip prosthesis during each simulation, analyzing them to obtain the force index of the hip prosthesis, and thereby outputting and displaying the effective simulation information. The specific steps include: obtaining the force parameters of the hip prosthesis during each simulation; obtaining the set of force demand factors of the hip prosthesis preset in the database, and comprehensively analyzing them with the force parameters of the hip prosthesis during each simulation to obtain the force index of the hip prosthesis for each simulation; obtaining the preset standard force index range of the hip prosthesis in the database, and comparing it with the force index of the hip prosthesis for each simulation. If the force index of the hip prosthesis for a certain simulation is within the standard force index range of the hip prosthesis, then define this simulation of the hip prosthesis as an effective simulation, and output and display the effective simulation information.

[0018] Optionally, a device for selecting the position of a hip prosthesis based on mechanical simulation, characterized by comprising: a CT image acquisition module, a CT image result analysis module, an overall imaging verification module, a hip prosthesis implantation simulation inspection execution judgment module, and a result analysis module; the CT image acquisition module is used to obtain a CT image when the mechanical simulation platform receives a CT upload signal, segment the CT image to obtain a hip joint region image, and further segment the hip joint region image to obtain CT images of each hip joint sub-region; the CT image result analysis module is used to analyze 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 to splice the CT images of each hip joint sub-region after image enhancement processing to obtain an overall CT image of the hip joint region, analyze to obtain the imaging quality index of the overall CT image of the hip joint region, and judge the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint region. The overall imaging verification result includes overall imaging verification qualified and overall imaging verification unqualified; the hip prosthesis implantation simulation inspection execution judgment module is used to make a judgment based on the overall imaging verification result. If the overall imaging verification result is overall imaging verification qualified, then input the overall CT image of the hip joint region into the three-dimensional reconstruction module in the mechanical simulation platform for hip prosthesis implantation simulation inspection; the result analysis module is used to count the hip prosthesis implantation simulation inspection data, process it to obtain a corrected region image set, and process the corrected region image set to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region.

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

[0020] 1. A method for selecting the position of a hip joint prosthesis based on mechanical simulation provided by the present invention obtains the hip joint region image through CT image segmentation, extracts the CT image performance parameters of each hip joint sub-region for analysis and processing, thereby obtaining the imaging quality index of the overall CT image of the hip joint region, and further obtaining the overall imaging verification result, effectively solving the problem of unreasonable prosthesis position selection caused by unqualified image quality processing in the prior art.

[0021] 2. The present invention discriminates based on the overall imaging verification result. If the overall imaging verification result is qualified for overall imaging verification, the overall CT image of the hip joint region is input into the three-dimensional reconstruction module for hip joint prosthesis implantation simulation inspection, thereby obtaining the prosthesis position simulation analysis result of the overall CT image of the hip joint region, and further outputting and displaying the effective simulation information.

[0022] 3. The present invention statistically analyzes the hip joint prosthesis implantation simulation inspection data, combines the CT image performance indicators and implantation performance influencing factors of each implantation sub-region for comprehensive analysis, thereby obtaining the deviation adjustment indicators of each implantation sub-region, and further obtaining the corrected region image set for processing. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0024] Figure 1 It is a flowchart of a method for selecting the position of a hip joint prosthesis based on mechanical simulation of the present invention;

[0025] Figure 2 It is a structural diagram of a device for selecting the position of a hip joint prosthesis based on mechanical simulation of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0026] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the following will clearly and completely describe the technical solutions of the embodiments of the present invention in conjunction with the drawings of the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art without creative efforts based on the described embodiments of the present invention fall within the scope of protection of the present invention.

[0027] Such as Figure 1As shown in the figure, it is a flowchart of a method for selecting the position of a hip joint prosthesis based on mechanical simulation provided by an embodiment of the present invention. This method is applied to a device for selecting the position of a hip joint prosthesis based on mechanical simulation. The method includes the following steps: S1. After the mechanical simulation platform receives the CT upload signal, obtain the CT image, segment the CT image to obtain the hip joint region image, and segment the hip joint region image to obtain the CT images of each hip joint sub-region; S2. Obtain the CT image performance parameters of each hip joint sub-region, analyze them to obtain the CT image performance indicators of each hip joint sub-region, thereby process to obtain the CT image analysis results of each hip joint sub-region, and perform corresponding image enhancement processing based on the CT image analysis results; S3. Obtain the CT images of each hip joint sub-region after image enhancement processing, perform stitching processing to obtain the overall CT image of the hip joint region, analyze to obtain the imaging quality indicator of the overall CT image of the hip joint region, and judge the overall imaging verification result according to the imaging quality indicator of the overall CT image of the hip joint region. The overall imaging verification result includes that the overall imaging verification is qualified and the overall imaging verification is unqualified;

[0028] S4. Based on the overall imaging verification result, make a discrimination. If the overall imaging verification result is that the overall imaging verification is qualified, input the overall CT image of the hip joint region into the three-dimensional reconstruction module in the mechanical simulation platform for hip joint prosthesis implantation simulation test; S5. Statistically process the hip joint prosthesis implantation simulation test data to obtain the corrected region image set, and process the corrected region image set to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region.

[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] By segmenting and processing the hip joint region image, key CT image performance parameters can be extracted. The stitching of the overall CT image and the determination of the imaging quality indicator provide a basis for the accurate selection of the prosthesis position. The mechanical simulation platform can effectively analyze the stress state of the prosthesis at different positions, and through comprehensive analysis with the set of stress requirement factors of the hip joint prosthesis preset in the database, obtain the stress simulation indicators of the hip joint prosthesis for each simulation, and then analyze to obtain effective simulation information.

[0031] Among them, the CT image performance parameters of each hip joint sub-region are obtained for analysis to obtain the CT image performance indicators of each hip joint sub-region. The specific steps include: obtaining the CT image performance parameters of each hip joint sub-region, where the CT image performance parameters of each hip joint sub-region include the average gray value, average density, noise value, resolution, and contrast of each hip joint sub-region; obtaining the preset CT image reference set in the database and comparing it with the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indicators of each hip joint sub-region; the CT image performance indicators of each hip joint sub-region are used to characterize the compliance degree of the CT image quality of each hip joint sub-region.

[0032] In this 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 through statistical tools or programming libraries (such as MATLAB, Python, and ImageJ / Fiji) in image processing software.

[0033] Among them, the specific method for obtaining the CT image performance indicators of each hip joint sub-region is as follows:

[0034]

[0035] In the formula, BX i represents the CT image performance indicator 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 the 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, ΔDB represents the contrast reference value.

[0036] In this embodiment, by analyzing the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indicators of each hip joint sub-region, the mutual influence relationship between these parameters is considered. For example, a higher average gray value indicates a higher image brightness, and a lower contrast will make the gray value differences of different tissues less obvious, affecting the resolution of the 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 relatively high, it helps to distinguish bones from soft tissues and improve the 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 differences between tissues less obvious. A higher image resolution can usually capture more details, but if the noise is large, the image details will be masked by the noise. On the contrary, a low resolution cannot clearly display all tissue details, and a higher resolution can more clearly display the tissue boundaries, thereby enhancing the contrast. If the contrast is too low, even with a high resolution, 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 the average gray value, average density, noise value, resolution, and contrast, and comparing them with the preset CT image reference set in the database, the image quality of each sub-region is evaluated, which can quantify the degree of compliance of the image quality of each sub-region, ensuring that the CT image meets the accuracy requirements and avoiding the problem of inaccurate subsequent hip joint prosthesis analysis caused by unqualified image quality. By timely detecting and processing image problems before hip joint prosthesis implantation, the accuracy and reliability of prosthesis position selection are improved.

[0038] Among them, the CT image analysis results of each hip joint sub-region are obtained through such 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 it 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 this hip joint sub-region is unqualified for CT image analysis. 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 this hip joint sub-region is qualified for CT image analysis; if the CT image analysis result is qualified for CT image analysis, the image enhancement processing is not performed. If the CT image analysis result is unqualified for CT image analysis, the image enhancement processing is performed.

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

[0040] By performing image enhancement processing based on the CT image analysis results in the present invention, 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 sub-region with the preset thresholds, it is possible to quickly identify whether the CT images of each hip sub-region meet the quality standards, avoiding the waste of computing resources caused by using unqualified images for prosthesis position selection. For unqualified images, the clarity, contrast, and details of the images are enhanced through image enhancement processing, effectively reducing the impact of image quality problems on the accuracy of prosthesis positioning. By automatically determining whether image enhancement processing is required, 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] Among them, the imaging quality index of the overall CT image of the hip joint region is obtained, specifically including: obtaining the CT images of each hip sub-region after image enhancement processing, and splicing them in the order before 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, and the imaging quality parameters of the overall CT image of the hip joint region include the artifact area ratio, the gray value range, the edge region gradient, the image noise, and the image granularity; obtaining the preset imaging quality reference set in the database, and comparing it 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 characterize the degree of compliance of the imaging quality of the overall CT image of the hip joint region.

[0042] In this embodiment, it should be noted 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 method for obtaining the imaging quality index of the overall CT image of the hip joint region 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 ratio of the artifact area to the total area. Artifacts are usually caused by equipment errors or irregular noises during 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 area, the imaging quality index of the overall CT image of the hip joint area is obtained. Considering the mutual influence relationship between these parameters, for example: the greater the image noise, the more artifacts will appear, 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 in the edge area, which is more conducive to improving the clarity and tissue differentiation of the image. The edge region gradient represents the distinctness of the structures 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. Images with lower granularity help to clearly present edge details. Larger image granularity is likely to generate noise and artifacts, thus reducing the image quality. The increase in image noise will affect the stability of the gray value, resulting in uneven gray value changes.

[0047] Among them, judging the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint area, the specific steps include: obtaining the imaging quality threshold of the overall CT image of the hip joint area preset in the database and comparing it with the imaging quality index of the overall CT image of the hip joint area. If the imaging quality index of the overall CT image of the hip joint area is above the imaging quality threshold of the overall CT image of the hip joint area, the overall imaging verification result is that the overall imaging verification is qualified. If the imaging quality index of the overall CT image of the hip joint area is less than the imaging quality threshold of the overall CT image of the hip joint area, the overall imaging verification result is that the overall imaging verification is unqualified.

[0048] In this embodiment, by comparing the imaging quality index of the overall CT image of the hip joint area with the preset quality threshold, it is possible to effectively judge whether the image quality meets the standard. If the image quality meets the standard, further mechanical simulation can be carried out on the basis of high-quality images, improving the accuracy and reliability of the simulation results. If the image quality does not meet the standard, it can be promptly identified and corresponding image processing can be carried out, avoiding the influence of image quality problems on the accuracy of prosthesis position selection, effectively improving the automation level of image quality control, reducing manual intervention, and ensuring the efficiency and accuracy of subsequent processing.

[0049] Among them, discrimination is performed based on the overall imaging verification result, specifically including: if the overall imaging verification result is qualified for overall imaging verification, the overall CT image of the hip joint area is input into the three-dimensional reconstruction module in the mechanical simulation platform for hip prosthesis implantation simulation inspection; if the overall imaging verification result is unqualified for overall imaging verification, the CT image is re-obtained, and steps 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 qualified for overall imaging verification, the repeated execution operation is stopped; if the overall imaging verification result after repeatedly executing S1 and S2 within the preset number of times is unqualified for overall imaging verification, a warning prompt message is issued.

[0050] In this embodiment, when judging whether the image quality is qualified through the overall imaging verification result, if the image quality is qualified, the overall CT image of the hip joint area can be directly input into the three-dimensional reconstruction module in the mechanical simulation platform for simulation inspection of hip prosthesis implantation. Through this step, the rationality of prosthesis implantation can be deeply analyzed, and the rationality of prosthesis implantation can be improved. If the image quality fails to pass the overall imaging verification, the step of re-obtaining the CT image and performing image processing will be automatically triggered. This process will repeatedly execute the image acquisition and analysis steps (S1 and S2) according to the preset number of times to ensure qualified image data after multiple attempts. If the image quality still fails to meet the standard after multiple attempts, a warning prompt message will be issued, and the warning prompt message includes the warning reason, warning time, etc.

[0051] Among them, after processing the corrected region image set, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is obtained. The specific steps include: obtaining the overall CT image of the hip joint area after hip prosthesis implantation, performing segmentation to obtain the CT images of each hip joint sub-region after hip prosthesis implantation, and marking them as the CT images of each implanted sub-region; obtaining the CT image performance parameters of each implanted sub-region, and performing analysis to obtain the CT image performance indicators of each implanted sub-region; obtaining the preset implantation performance influencing factors in the database, and performing comprehensive analysis based on the implantation performance influencing factors, the CT image performance indicators of each hip joint sub-region, and the CT image performance indicators of each implanted sub-region to obtain the deviation adjustment indicators of each implanted sub-region; obtaining the preset deviation adjustment threshold in the database, and comparing it with the deviation adjustment indicators of each implanted sub-region. If the deviation adjustment indicator of a certain implanted sub-region is above the deviation adjustment threshold, the implanted sub-region is marked as the corrected region image set; after processing the corrected region image set, if the deviation adjustment indicator of the corrected region image set is above the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as unqualified for simulation verification; if the deviation adjustment indicator of the corrected region image set is less than the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as qualified for simulation verification.

[0052] In this embodiment, it should be noted that the prosthesis position simulation analysis results include qualified simulation verification and unqualified simulation verification. The corrected region image set includes the images of each implant sub-region corresponding to the deviation adjustment index of the implant sub-region being above the deviation adjustment threshold.

[0053] It should be noted that if the prosthesis position simulation analysis result of the overall CT image of the hip joint region is qualified for simulation verification, the force simulation of the hip joint prosthesis will be carried out in the next step. If the prosthesis position simulation analysis result of the overall CT image of the hip joint region is unqualified for simulation verification, the force simulation will not be carried out and a warning prompt message will be sent. The warning prompt message includes a prompt indicating that the simulation verification is unqualified and an alarm bell will ring.

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

[0055] It should also be noted that the CT image performance parameters of each implant sub-region can be processed in the same way as the CT image performance indicators of each hip joint sub-region to obtain the CT image performance indicators of each implant sub-region.

[0056] The deviation adjustment index of each implant sub-region is obtained. 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 the implant sub-region, 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 the hip joint sub-region, 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 the natural constant, and Y z represents the implant performance influence factor.

[0059] Through the processing of segmenting, analyzing, and correcting the regional image set of the overall CT image after hip prosthesis implantation, the present invention can effectively ensure the accuracy and stability of the prosthesis position. After obtaining the CT image after implantation and analyzing each implanted sub-region, combined with the preset implantation performance influencing factors in the database, it is possible to comprehensively evaluate the adaptability of the prosthesis to the image quality of the hip joint region, and optimize the deviation adjustment index according to the analysis results to ensure that the image quality of each sub-region is optimized. Among them, the preset implantation performance influencing factors in the database are used to represent the influence degree after hip prosthesis implantation. By comparing the deviation adjustment index with the preset threshold, the area that needs to be corrected can be automatically identified and marked. The automated correction and verification mechanism reduces manual intervention and improves the efficiency and accuracy of the operation.

[0060] Among them, it also includes obtaining the force parameters of the hip prosthesis during each simulation, and analyzing to obtain the force index of the hip prosthesis, and thus outputting and displaying the effective simulation information. The specific steps include: obtaining the force parameters of the hip prosthesis during each simulation; obtaining the set of force requirement factors of the hip prosthesis preset in the database, and comprehensively analyzing them with the force parameters of the hip prosthesis during each simulation to obtain the force index of the hip prosthesis during each simulation; obtaining the preset standard force index range of the hip prosthesis in the database, and comparing it with the force index of the hip prosthesis during each simulation. If the force index of the hip prosthesis during a certain simulation is within the standard force index range of the hip prosthesis, then this simulation of the hip prosthesis is defined as an effective simulation, and the effective simulation information is output and displayed. If the force index of the hip prosthesis during a certain simulation is not within the standard force index range of the hip prosthesis, then this simulation of the hip prosthesis is defined as an invalid simulation and is not displayed.

[0061] In this embodiment, the force parameters of the hip prosthesis include the reference value of the friction force between the prosthesis and the acetabular contact surface, the contact stress between the acetabulum and the prosthesis, the shear stress between the prosthesis and the bone tissue, the normal force between the prosthesis and the acetabulum, the torsional force received by the hip joint during flexion and extension, the abduction and adduction moment of the hip joint, and the internal and external rotation moment of the hip joint.

[0062] The set of force requirement factors of the hip prosthesis includes: the unit contact stress requirement factor between the acetabulum and the prosthesis, the unit shear stress requirement factor between the prosthesis and the bone tissue, the unit normal force requirement factor between the prosthesis and the acetabulum, the unit abduction and adduction moment requirement factor of the hip joint, the unit internal and external rotation moment requirement factor of the hip joint, the reference value of the friction force between the prosthesis and the acetabular contact surface, and the reference value of the torsional force received by the hip joint during flexion and extension.

[0063] Contact stress refers to the pressure applied between the prosthesis and the acetabular contact surface. Shear stress refers to the lateral stress generated at the contact surface between the prosthesis and the bone tissue due to movement and the action of forces. The normal force is the force perpendicular to the contact surface between the prosthesis and the acetabulum, representing the pressure distribution between the prosthesis and the acetabulum. Frictional force refers to the resistance generated between the contact surfaces of the prosthesis and the acetabulum due to contact and movement. The flexion-extension torsion moment refers to the moment experienced by the hip joint during flexion and extension movements, generated by the contact relationship and movement between the femoral head and the acetabulum. The abduction-adduction moment refers to the moment applied during abduction or adduction movements of the hip joint. The internal-external rotation moment refers to the moment applied during internal and external rotation of the hip joint, located between the femoral head and the acetabulum.

[0064] By analyzing the force parameters of the hip joint prosthesis during each simulation, the force indices for each simulation of the hip joint prosthesis are obtained. This is considering the mutual influence relationships among these parameters. For example, the greater the normal force, the larger the contact area between the prosthesis and the acetabulum, thus leading to an increase in contact stress. When the normal force increases, the influence of shear stress will be exacerbated. Excessive shear stress can cause separation between the prosthesis and the bone tissue, affecting contact stability. When the normal force is relatively large, the effect of shear stress is increased, leading to displacement of the prosthesis and the risk of fracture. Appropriate frictional force can provide stable joint movement, but excessive frictional force will cause wear, while too small frictional force will lead to prosthesis instability. A relatively large flexion-extension moment will increase the torsion moment, thus affecting the stability of the prosthesis and resulting in rotation or displacement of the prosthesis. Excessive abduction-adduction moment will generate an excessive rotation moment, leading to prosthesis displacement or dislocation. Contact stress and shear stress are usually in the same direction or interact with each other. Larger contact stress will exacerbate the influence of shear stress. When contact stress increases, shear stress can cause damage to bone tissue or dislocation of the prosthesis, affecting joint stability.

[0065] The method for obtaining the force indices for each simulation of the hip joint prosthesis is as follows:

[0066]

[0067] In the formula, L n represents the force index of the hip joint prosthesis for the nth simulation, n represents the numbering of the 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 for the nth simulation, ΔKY represents the unit contact stress requirement factor between the acetabulum and the prosthesis, KQ n represents the shear stress between the prosthesis and the bone tissue of the hip joint prosthesis for the nth simulation, ΔKQ represents the unit shear stress requirement factor between the prosthesis and the bone tissue of the hip joint prosthesis, KF nDenote the normal force between the prosthesis and the acetabulum of the hip prosthesis during the nth simulation, ΔKF denotes the unit normal force demand factor between the prosthesis and the acetabulum, and KM n Denote the frictional force between the prosthesis and the acetabulum contact surface of the hip prosthesis during the nth simulation, ΔKM denotes the reference value of the frictional force between the prosthesis and the acetabulum contact surface, and KN n Denote the torsional force exerted on the hip flexion and extension of the hip prosthesis during the nth simulation, ΔKN denotes the reference value of the torsional force exerted on the hip flexion and extension, and KS n Denote the abduction and adduction moment of the hip of the hip prosthesis during the nth simulation, ΔKS denotes the unit abduction and adduction moment demand factor of the hip, and KX n Denote the internal and external rotation moment of the hip of the hip prosthesis during the nth simulation, and ΔKX denotes the unit internal and external rotation moment demand factor of the hip.

[0068] It should be noted that the effective simulation information includes the force condition of the hip prosthesis, the installation position of the hip prosthesis, and the installation angle of the hip prosthesis, etc. The force indexes of each simulation of the hip prosthesis are used to characterize the force degree of the hip prosthesis during each simulation force process.

[0069] The processing of the corrected region image set can be carried out by using image registration and fusion algorithms. The specific method is to register and fuse the image of the corrected region with the CT image of any hip sub-region whose CT image analysis result is qualified for CT image analysis, and it is processed by using feature point-based registration methods such as SIFT (Full name: Scale Invariant Feature Transform), ORB (Full name: Oriented FAST and Rotated BRIEF).

[0070] As Figure 2 shown, it is the structural diagram of a hip prosthesis position selection device based on mechanical simulation provided by an embodiment of the present invention. A hip prosthesis position selection device based on mechanical simulation provided by an embodiment of the present invention includes:

[0071] A CT image acquisition module, a CT image result analysis module, an overall imaging verification module, a hip prosthesis implantation simulation test execution judgment module, and a result analysis module; the CT image acquisition module is used to acquire a CT image when the mechanical simulation platform receives a CT upload signal, segment the CT image to obtain a hip joint region image, and segment the hip joint region image to obtain CT images of each hip joint sub-region; the CT image result analysis module is used to analyze the CT image performance parameters of each hip joint sub-region to obtain the CT image performance indicators of each hip joint sub-region, thereby processing to obtain the CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; the overall imaging verification module is used to splice the CT images of each hip joint sub-region after image enhancement processing to obtain an overall CT image of the hip joint region, analyze to obtain the imaging quality indicators of the overall CT image of the hip joint region, and judge the overall imaging verification result according to the imaging quality indicators of the overall CT image of the hip joint region. The overall imaging verification result includes that the overall imaging verification is qualified and the overall imaging verification is unqualified; the hip prosthesis implantation simulation test execution judgment module is used to make a judgment based on the overall imaging verification result. If the overall imaging verification result is that the overall imaging verification is qualified, the overall CT image of the hip joint region is input into the three-dimensional reconstruction module in the mechanical simulation platform for hip prosthesis implantation simulation test; the result analysis module is used to count the hip prosthesis implantation simulation test data, process it to obtain a corrected region image set, and process 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 the specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. The protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A hip joint prosthesis position selection method based on mechanical simulation, characterized in that: include: S1. After receiving the CT upload signal, the mechanical simulation platform obtains a CT image, segments the CT image, obtains a hip joint region image, and segments the hip joint region image to obtain CT images of each hip joint sub-region; S2, obtaining CT image performance parameters of each hip joint sub-region for analysis to obtain CT image performance indicators of each hip joint sub-region, thereby obtaining CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; S3, obtaining the CT images of each hip joint sub-region after image enhancement processing, performing splicing processing to obtain an overall CT image of the hip joint region, and analyzing to obtain an imaging quality index of the overall CT image of the hip joint region, and judging an overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint region, wherein the overall imaging verification result includes an overall imaging verification qualified and an overall imaging verification unqualified; S4, making a judgment based on the overall imaging verification result, if the overall imaging verification result is qualified, inputting the overall CT image of the hip joint area into the three-dimensional reconstruction module in the mechanical simulation platform to perform hip joint prosthesis implantation simulation inspection; S5. Statistically process the hip joint prosthesis implantation simulation test data to obtain a correction region image set, and process the correction region image set to obtain a prosthesis position simulation analysis result of the overall CT image of the hip joint region.

2. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The CT image performance parameters of each hip joint sub-region are obtained and analyzed to obtain the CT image performance index of each hip joint sub-region, and the specific steps include: Acquire CT image performance parameters of each hip joint sub-region, wherein 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; Obtaining a preset CT image reference set in the database, and comparing it with the CT image performance parameters of each hip joint sub-region to obtain the CT image performance index of each hip joint sub-region; The CT image performance index of each hip joint sub-region is used to characterize the degree of compliance of the CT image quality of each hip joint sub-region.

3. The hip joint prosthesis position selection method based on mechanical simulation according to claim 2, characterized in that: The specific method for obtaining the CT image performance index of each hip joint sub-region is as follows: 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 the natural constant, HD i represents the average gray value of the i-th hip joint sub-region, ΔHD represents the average gray reference value, and ρ i represents the average density of the i-th hip joint sub-region, Δρ represents the average density reference value, and ZS i represents the noise value of the i-th hip joint sub-region, ΔZS represents the noise reference value, and 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.

4. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The CT image analysis results of each hip joint sub-region are obtained by this processing, and corresponding image enhancement processing is performed based on the CT image analysis results. The specific steps include: The CT image performance threshold preset in the database is obtained, and compared with the CT image performance index of each hip joint sub-region. If the CT image performance index of a hip joint sub-region is less than the CT image performance threshold, the CT image analysis result of the hip joint sub-region is unqualified; if the CT image performance index of a hip joint sub-region is above the CT image performance threshold, the CT image analysis result of the hip joint sub-region is qualified; If the CT image analysis result is that the CT image analysis is qualified, the image enhancement processing is not performed. If the CT image analysis result is that the CT image analysis is unqualified, the image enhancement processing is performed.

5. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The analysis obtains the imaging quality index of the overall CT image of the hip joint area, specifically including: Obtaining CT images of each hip joint sub-region after image enhancement processing, and splicing them in the order before segmentation to obtain an overall CT image of the hip joint region; Acquiring imaging quality parameters of the overall CT image of the hip joint region, wherein the imaging quality parameters of the overall CT image of the hip joint region include artifact area ratio, gray value extreme difference, edge region gradient, image noise and image granularity; Obtaining a preset imaging quality reference set in a database, and comparing it with the imaging quality parameters of the overall CT image of the hip joint region, to obtain an 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 characterize the conformity of the imaging quality of the overall CT image of the hip joint region.

6. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The step of judging the overall imaging verification result according to the imaging quality index of the overall CT image of the hip joint region specifically comprises: The imaging quality threshold of the overall CT image of the hip joint region preset in the database is obtained, and compared 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 that the overall imaging verification is qualified; 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 that the overall imaging verification is unqualified.

7. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The determination based on the overall imaging verification result specifically includes: If the overall imaging verification result is that the overall imaging verification is qualified, the overall CT image of the hip joint area is input into the three-dimensional reconstruction module in the mechanical simulation platform to perform hip prosthesis implantation simulation inspection; if the overall imaging verification result is that the overall imaging verification is unqualified, the CT image is re-acquired, and S1 and S2 are repeated within a preset number of times; If the overall imaging verification result after repeating S1 and S2 within the preset number of times is that the overall imaging verification is qualified, the repeated execution operation is stopped; if the overall imaging verification result after repeating S1 and S2 within the preset number of times is that the overall imaging verification is unqualified, an early warning prompt message is issued.

8. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: The process of processing the corrected region image set to obtain the prosthesis position simulation analysis result of the overall CT image of the hip joint region specifically comprises: Obtaining an overall CT image of the hip joint region after hip prosthesis implantation, segmenting to obtain CT images of each hip joint sub-region after hip prosthesis implantation, marking them as CT images of each implantation sub-region, obtaining CT image performance parameters of each implantation sub-region, and analyzing to obtain CT image performance indicators of each implantation sub-region; Obtaining implant performance influencing factors preset in the database, and performing comprehensive analysis based on the implant performance influencing factors, the CT image performance indicators of each hip joint sub-region, and the CT image performance indicators of each implant sub-region, to obtain the deviation adjustment indicator of each implant sub-region; Obtaining a deviation adjustment threshold preset in a database and comparing it with the deviation adjustment index of each implantation sub-region, if there is a deviation adjustment index of an implantation sub-region that is above the deviation adjustment threshold, marking the implantation sub-region as a correction region image set; After processing the correction area image set, if the deviation adjustment index of the correction area image set is above the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as unqualified simulation verification; if the deviation adjustment index of the correction area image set is less than the deviation adjustment threshold, the prosthesis position simulation analysis result of the overall CT image of the hip joint area is classified as qualified simulation verification.

9. The hip joint prosthesis position selection method based on mechanical simulation according to claim 1, characterized in that: It also includes obtaining the force parameters of the hip joint prosthesis during each simulation, and analyzing to obtain the force index of the hip joint prosthesis, thereby outputting and displaying the effective simulation information. The specific steps include: Obtain the force parameters of the hip joint prosthesis during each simulation; Obtain a set of force demand factors of the hip joint prosthesis preset in the database, and conduct a comprehensive analysis with the force parameters of the hip joint prosthesis during each simulation to obtain the force index of each simulation of the hip joint prosthesis; The preset standard hip prosthesis force index range in the database is obtained, and compared with each simulated force index of the hip prosthesis. If a simulated force index of the hip prosthesis is within the standard hip prosthesis force index range, the simulation of the hip prosthesis is defined as a valid simulation, and the valid simulation information is output and displayed.

10. A hip joint prosthesis position selection device based on mechanical simulation, characterized in that: include: CT image acquisition module, CT image result analysis module, overall imaging verification module, hip prosthesis implantation simulation test execution judgment module and result analysis module; The CT image acquisition module is used to acquire a CT image after the mechanical simulation platform receives the CT upload signal, and segment the CT image to obtain a hip joint area image, and segment the hip joint area image to obtain CT images of each hip joint sub-area; 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 obtaining the CT image analysis results of each hip joint sub-region, and performing corresponding image enhancement processing based on the CT image analysis results; The overall imaging verification module is used to obtain the CT images of each hip joint sub-region after image enhancement processing, perform splicing processing, obtain the overall CT image of the hip joint region, and analyze to obtain the imaging quality index of the overall CT image of the hip joint region, and judge 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 the overall imaging verification qualified and the overall imaging verification unqualified; The hip joint prosthesis implant simulation test execution judgment module is used to make a judgment based on the overall imaging verification result. If the overall imaging verification result is that the overall imaging verification is qualified, the overall CT image of the hip joint area is input into the three-dimensional reconstruction module in the mechanical simulation platform to perform hip joint prosthesis implant simulation test; The result analysis module is used to collect statistics on hip joint prosthesis implantation simulation inspection data for processing to obtain a correction area image set, and to process the correction area image set to obtain a prosthesis position simulation analysis result of the overall CT image of the hip joint area.

Images