

技术领域technical field
本发明涉及医疗信息化技术领域,尤其涉及一种基于股骨旋转轴与内翻轴的分析模型的构建方法及系统。The invention relates to the technical field of medical informatization, in particular to a method and system for constructing an analysis model based on a femoral rotation axis and a varus axis.
背景技术Background technique
股骨头坏死是一个病理演变过程,初始发生在股骨头的负重区,应力作用下坏死骨骨小梁结构发生损伤即显微骨折以及随后针对损伤骨组织的修复过程。造成骨坏死的原因不消除,修复不完善,损伤-修复的过程继续,导致股骨头结构改变、股骨头塌陷、变形,关节炎症,功能障碍。对这一疾病如果单纯的从临床症状及体征上进行确诊比较困难,容易出现误诊或者漏诊的情况,现阶段临床上主要采取的检测方法是X线检查及CT检查、磁共振检查3种诊断方式。在股骨头坏死的临床诊断中,X线检查是应用最广泛的一种方法,通过检查主要能够对患者的病情进展及骨坏死程度进行了解,CT检查也主要是对患者骨坏死的情况进行了解,通过对影像进行分析选择合适的治疗方法。Necrosis of the femoral head is a pathological evolution process that initially occurs in the weight-bearing area of the femoral head. Under the action of stress, the trabecular structure of necrotic bone is damaged, that is, microfracture and the subsequent repair process for damaged bone tissue. If the cause of osteonecrosis is not eliminated, the repair is not perfect, and the process of injury-repair continues, resulting in structural changes of the femoral head, collapse, deformation of the femoral head, joint inflammation, and dysfunction. It is difficult to diagnose this disease purely from clinical symptoms and signs, and it is prone to misdiagnosis or missed diagnosis. At this stage, the main clinical detection methods are X-ray examination, CT examination, and magnetic resonance examination. . In the clinical diagnosis of femoral head necrosis, X-ray examination is the most widely used method. The examination can mainly understand the progress of the patient's disease and the degree of osteonecrosis. CT examination is also mainly used to understand the situation of the patient's osteonecrosis. , select the appropriate treatment method by analyzing the images.
股骨头坏死(avascular necrosis of the femoral head,ANFH)主要采用以下方式治疗:非手术治疗、保髋手术治疗和人工关节置换术。随着技术计算机技术、空间导航技术以及图像处理技术的迅速发展,计算机辅助导航技术在传统的TKR手术方法中得到了较好应用。计算机辅助导航技术系统是将术前和术后获得的影像学资料经过计算机处理后形成三维可视影像。其可以进行模拟手术,术前可自动选出最适合的尺寸关节零件、模拟开刀,手术中在骨头中放置感应器,导航系统上有红外线摄像头,动态追踪手术器械相对患者解剖的实时空间位置关系,手术医生据此可精确地放置人工关节。不过导航技术存在如下地不足:导航技术无法确定负重片区,即导航技术是在非负重的情况下测量下肢力线,其结果与负重下的存在差异。Avascular necrosis of the femoral head (ANFH) is mainly treated by the following methods: non-surgical treatment, hip-sparing surgery and artificial joint replacement. With the rapid development of technical computer technology, spatial navigation technology and image processing technology, computer-aided navigation technology has been well applied in traditional TKR surgical methods. The computer-aided navigation technology system is to form a three-dimensional visual image after the imaging data obtained before and after the operation are processed by the computer. It can perform simulated surgery. Before surgery, it can automatically select the most suitable size joint parts and simulate surgery. During the surgery, sensors are placed in the bones. There is an infrared camera on the navigation system to dynamically track the real-time spatial position relationship of the surgical instruments relative to the patient's anatomy. , the surgeon can accurately place the artificial joint accordingly. However, the navigation technology has the following shortcomings: the navigation technology cannot determine the weight-bearing area, that is, the navigation technology measures the lower extremity alignment under the condition of non-weight-bearing, and the results are different from those under the weight-bearing condition.
例如,公开号为CN108542408A的中国专利公开的一种三维立体股骨头尺寸测量装置。其主要包括:扫描床、防护装置、测量装置。扫描床第一段安装有防护装置,另一端安装测量装置,扫面装置安装在扫面床的侧面。该发明的三维立体股骨头尺寸测量装置将X射线扫描、双目立体摄像、三维建模技术相结合,准确地测定患者的股骨头的三维立体尺寸,达到精确测量的目的,在进行人工股骨头置换手术之前,对患者的股骨头进行三维建模,在根据患者的实际治疗情况,医护人员制定出最佳手术方案,结构巧妙,操作简单快捷,提高医疗人员的工作效率和患者的治愈率。For example, a three-dimensional three-dimensional femoral head size measurement device is disclosed in Chinese Patent Publication No. CN108542408A. It mainly includes: scanning bed, protective device, measuring device. A protective device is installed on the first section of the scanning bed, a measuring device is installed at the other end, and the scanning device is installed on the side of the scanning bed. The three-dimensional three-dimensional femoral head size measuring device of the invention combines X-ray scanning, binocular stereo camera, and three-dimensional modeling technology to accurately measure the three-dimensional three-dimensional size of the patient's femoral head to achieve the purpose of accurate measurement. Before the replacement surgery, 3D modeling of the patient's femoral head is carried out. According to the actual treatment situation of the patient, the medical staff formulates the best surgical plan. The structure is ingenious, the operation is simple and fast, and the work efficiency of the medical staff and the cure rate of the patient are improved.
例如,公开号为CN107545578A的中国专利公开的一种CT图像中股骨头区域分割方法、装置及设备。该方法包括获取被检体股骨头区域CT序列图像,通过对所述CT序列图像的每一幅图像的股骨头区域进行定位,构建股骨头包围盒和股骨颈包围盒;基于空间辐射线法在股骨头包围盒和股骨颈包围盒内获取多个采样点;根据多个采样点生成二维采样图像;在二维采样图像上确定股骨头区域的骨皮质外边界。该发明旨在准确地分割CT序列图像中股骨头区域的骨皮质外边界。For example, Chinese Patent Publication No. CN107545578A discloses a method, device and device for segmenting femoral head region in CT images. The method includes acquiring a CT sequence image of the femoral head region of a subject, and constructing a femoral head bounding box and a femoral neck bounding box by locating the femoral head region of each image of the CT sequence image; Obtain multiple sampling points in the femoral head bounding box and the femoral neck bounding box; generate a two-dimensional sampling image according to the multiple sampling points; determine the outer boundary of the bone cortex of the femoral head region on the two-dimensional sampling image. The invention aims to accurately segment the outer cortical boundary of the femoral head region in CT sequence images.
例如,公开号为CN104462636B的中国专利公开的一种基于伞状股骨头支撑器的坏死股骨头修复模型的建模方法。其方法包括:一、待修复股骨头的三维模型获取;获取待修复股骨头NURBS曲面模型,待修复股骨头为存在骨组织坏死区域且预采用股骨头支撑器修复的股骨头;股骨头支撑器由伞状支撑器和支撑套筒组成;二、根据伞状支撑器形状,确定需分离坏死区域,建立坏死股骨头模型;三、股骨头支撑器模型建立;四、坏死股骨头植入模型建立:建立带植入通道的坏死股骨头植入模型和植入骨的三维模型;五、坏死股骨头修复模型建立。该发明能够简便地、快速建立股骨头支撑器植入坏死股骨头的修复模型,且所建立的修复模型质量高。For example, Chinese Patent Publication No. CN104462636B discloses a modeling method for a necrotic femoral head repair model based on an umbrella-shaped femoral head supporter. The method includes: 1. obtaining a three-dimensional model of the femoral head to be repaired; obtaining a NURBS surface model of the femoral head to be repaired, the femoral head to be repaired is a femoral head with an area of bone tissue necrosis and pre-repaired with a femoral head supporter; It consists of an umbrella-shaped supporter and a support sleeve; 2. According to the shape of the umbrella-shaped supporter, determine the necrotic area to be separated, and establish a necrotic femoral head model; 3. Establish a femoral head supporter model; 4. Establish a necrotic femoral head implantation model : Establish the necrotic femoral head implantation model with implantation channel and the three-dimensional model of the implanted bone; 5. Establishment of the necrotic femoral head repair model. The invention can simply and quickly establish a repair model for implanting a femoral head supporter into a necrotic femoral head, and the established repair model is of high quality.
在科技论文研究方面,例如,《大连理工大学学报》2009年第49卷(第二期)的文章——基于CT图像反求技术的人体股骨头修复建模。利用基于模型的分割方法从人体髋关节CT图像中分别提取股骨和髋臼骨轮廓三维点数据,通过对离散点云数据进行规则化处理、精简、分割等,提取目标区域点云数据,最后利用最小二乘法拟合对因病变而发生塌陷的股骨头表面进行修复重建,恢复其健康状态的球面形态特征。利用两例不同程度的股骨头坏死病例进行修复实验,均得到较精确的修复还原效果。该方法不仅可以很好地重构人体髋关节模型、还原缺损股骨头的原始形态,而且可以推广到其他关节或骨组织,为骨科手术的精确定位和生物力学有限元分析提供了理论模型,并为关节假体的个体化快速制造提供了一种有效的方法。In terms of scientific and technological paper research, for example, the article in "Journal of Dalian University of Technology", 2009, volume 49 (second issue) - Modeling of human femoral head repair based on CT image reverse technology. The model-based segmentation method is used to extract the three-dimensional point data of the femur and acetabular bone contours from the CT images of the human hip, respectively. The least-squares fitting was used to repair and reconstruct the collapsed femoral head surface due to lesions, and restore the spherical morphological characteristics of its healthy state. Using two cases of different degrees of femoral head necrosis to carry out repair experiments, more accurate repair and restoration effects were obtained. This method can not only reconstruct the human hip joint model well and restore the original shape of the defective femoral head, but also can be extended to other joints or bone tissues, providing a theoretical model for accurate positioning and biomechanical finite element analysis of orthopaedic surgery, and also It provides an effective method for individualized and rapid manufacture of joint prostheses.
东北大学硕士论文提供了一种CT和MRI图像配准及融合方法。The master thesis of Northeastern University provides a method of CT and MRI image registration and fusion.
在其他相关技术领域方面,公开号为CN104622572A的基于医学图像的个性化骨科定位片。该专利旨在提供一种能够精确确定膝关节置换中关节定位孔和切削面的位置。首先,基于医学图像重建骨骼模型,对重构模型进行术前规划,确定下肢力线,股骨旋转轴和截骨参考点等参数,模拟截骨和假体置入。利用以上结果,确定和设计定位片;医生登录数据管理模块,查询数据。In other related technical fields, the publication number is CN104622572A, a medical image-based personalized orthopedic positioning sheet. This patent aims to provide a method that can accurately determine the location of joint positioning holes and cutting surfaces in knee replacements. First, reconstruct the skeletal model based on medical images, carry out preoperative planning on the reconstructed model, determine the lower limb force line, femoral rotation axis and osteotomy reference point and other parameters, and simulate osteotomy and prosthesis placement. Using the above results, determine and design the positioning sheet; the doctor logs in the data management module to query the data.
中国专利(公开号为CN108711187A)公开了一种配准融合CT和MRI信号建立人体腰椎三维仿真模型的方法,包括:采集计算机断层扫描CT图像;采集磁共振成像MRI图像;建立计算机断层扫描图像三维模型;建立磁共振成像图像三维模型;配准融合计算机断层扫描图像三维模型和磁共振成像图像三维模型,包括根据腰椎解剖结构进行简单配准和进行全局计算配准。该专利充分利用现有常规检查如计算机断层扫描和磁共振成像结合优选磁共振扫描序列,建立可以在各个磁共振序列内相互验证准确度的腰椎间盘、神经根、黄韧带的重要软组织三维模型,建立了腰椎间盘新的医学影像高准确度的建模方式,同时大大提高了医学影像检查的数据的利用率。The Chinese patent (publication number CN108711187A) discloses a method for establishing a three-dimensional simulation model of human lumbar spine by registering and fusing CT and MRI signals, including: acquiring computed tomography CT images; acquiring magnetic resonance imaging MRI images; establishing three-dimensional computed tomography images model; establishing a three-dimensional model of a magnetic resonance imaging image; registration and fusion of the three-dimensional model of the computed tomography image and the three-dimensional model of the magnetic resonance imaging image, including simple registration and global calculation registration according to the anatomical structure of the lumbar spine. This patent makes full use of existing routine examinations such as computed tomography and magnetic resonance imaging combined with the optimal magnetic resonance scanning sequence to establish a three-dimensional model of important soft tissue of lumbar intervertebral disc, nerve root, and ligamentum flavum that can mutually verify the accuracy within each magnetic resonance sequence. A new high-accuracy modeling method for medical imaging of lumbar intervertebral disc has been established, and the utilization rate of medical imaging data has been greatly improved.
对于截骨术而言,在手术前,外科医生需要获得对手术对象身体结构尽可能进行精确地解剖学图像,针对股骨头坏死的诊断主要来自于CT/MR的影像数据,才能确定手术方案,例如确定切开的大小、位置,如何避开其他脏器和神经等等。骨科医生只有根据CT/MR的影像来预估旋转角度以及旋转相应的角度后大致估计是否有将大部分坏死骨质区域移出负重区和将正常的骨质区域移入负重区的可能,而且即使得到了预估的结果,在实际操作过程中,截骨后,一些医生还可能要旋转到不同的角度去尝试是否可能有更佳的手术方案。因此,外壳医生基于CT或者MR影响数据难以准确地判定股骨头坏死的空间立体关系,这对于手术方案的确定、实施以及手术的成功率是不利的。For osteotomy, before the operation, the surgeon needs to obtain the most accurate anatomical images of the body structure of the operation object. The diagnosis of femoral head necrosis mainly comes from the image data of CT/MR to determine the operation plan. For example, determine the size and location of the incision, how to avoid other organs and nerves, and so on. Orthopedic surgeons can only estimate the rotation angle based on CT/MR images and roughly estimate whether it is possible to move most of the necrotic bone area out of the weight-bearing area and the normal bone area into the weight-bearing area after rotating the corresponding angle. According to the estimated results, during the actual operation, after the osteotomy, some doctors may also rotate to different angles to try to see if there may be a better surgical plan. Therefore, it is difficult for the shell doctor to accurately determine the spatial three-dimensional relationship of femoral head necrosis based on CT or MR influence data, which is unfavorable for the determination and implementation of the operation plan and the success rate of the operation.
发明内容SUMMARY OF THE INVENTION
针对现有技术之不足,本发明旨在使用计算机断层扫面技术,获取断层检查数据,生成三维可视化图像,用于能够在股骨头坏死的空间立体关系,从而能够基于该立体关系确定手术方案以及在手术中存在的风险进行评估。本发明能够在虚拟空间内,通过导入检查数据的形式,构造一个能够模拟人体构造的动态模型,通过三维切割、三维交互以及三维测量等手段,便于医生能够模拟手术,为患者手术提供一种合理的手术方案。In view of the deficiencies of the prior art, the present invention aims to use the computed tomography scanning technology to obtain the tomographic examination data and generate a three-dimensional visual image, which is used for the spatial three-dimensional relationship of the femoral head necrosis, so that the operation plan can be determined based on the three-dimensional relationship. Risks present during surgery are assessed. The present invention can construct a dynamic model capable of simulating human body structure in the virtual space by importing the form of examination data, and by means of three-dimensional cutting, three-dimensional interaction, three-dimensional measurement, etc., it is convenient for doctors to simulate surgery, and provides a reasonable operation for patients. surgical plan.
基于此,本发明提供了一种基于股骨旋转轴与内翻轴分析模型的构建系统,包括:模型构建模块,基于检查数据构建股骨-髋臼模型;图像显示模块,用于呈现所述股骨-髋臼模型;模型分割模块,能够基于所述股骨-髋臼模型生成髋臼模型和股骨模型;所述系统还包括:模型切割模块,基于所述股骨模型生成能够用于分割股骨头部模型和股骨体模型的切割面,旋转轴构建模块,其能够基于所述股骨头部模型和所述切割面构建旋转轴,其中,在所述股骨头部模型绕所述旋转轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动,内翻轴构建模块,其能够依据所述股骨头部模型的骨表面和所述旋转轴构建内翻轴,其中,在所述股骨头部模型绕所述内翻轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移除负重区的方式与所述髋关节面生成倾角。Based on this, the present invention provides a construction system based on the analysis model of the femoral rotation axis and varus axis, including: a model construction module, which constructs a femur-acetabular model based on inspection data; an image display module, which is used to present the femur-acetabular model. An acetabular model; a model segmentation module capable of generating an acetabular model and a femoral bone model based on the femur-acetabular model; the system further comprises: a model cutting module, based on the femoral model to generate a model that can be used to segment the femoral head and The cutting plane of the femoral body model, the rotation axis building block, which can construct the rotation axis based on the femoral head model and the cutting plane, wherein, when the femoral head model rotates around the rotation axis, the The femoral head model can be rotated in the hip joint surface in a manner of moving a partial area of the femoral head out of the weight-bearing area, and the varus axis builds a module, which can build varus based on the bone surface of the femoral head model and the rotation axis. an axis, wherein the femoral head model is capable of generating an inclination angle with the hip articular surface in a manner that removes a partial area of the femoral head from the weight bearing zone when the femoral head model is rotated about the varus axis.
根据一种优选的实施方式,所述系统还包括判别模块,用于判别所述旋转轴是否满足第一特征需求;如果所述旋转轴与所述切割面的法向量的角度小于第一角度阈值,则所述旋转轴满足所述第一特征需求;否则,所述旋转轴不满足所述第一特征需求;在所述判别模块判别出所述旋转轴不满足所述第一特征需求的情况下,所述模型切割模块基于所述股骨模型对所述股骨头部模型和股骨体模型进行切割以重新生成切割面。According to a preferred embodiment, the system further includes a judging module for judging whether the rotation axis meets the first characteristic requirement; if the angle between the rotation axis and the normal vector of the cutting plane is less than a first angle threshold , then the rotation axis satisfies the first feature requirement; otherwise, the rotation axis does not meet the first feature requirement; the discriminating module determines that the rotation axis does not meet the first feature requirement Next, the model cutting module cuts the femoral head model and the femoral body model based on the femoral model to regenerate cutting surfaces.
根据一种优选的实施方式,所述判别模块还用于判别所述切割面是否满足第二特征需求;如果所述切割面的法向量与所述髋关节面形成的空间体的重心线小于第二角度阈值,则所述切割面满足所述第二特征需求;否则,所述切割面不满所述第二特征需求;在所述判别模块判别出所述切割面不满足所述第二特征需求的情况下,所述模型切割模块基于所述股骨模型对所述股骨头部模型和股骨体模型进行切割以重新生成切割面。According to a preferred embodiment, the judging module is further configured to judge whether the cutting plane meets the second characteristic requirement; Two angle thresholds, the cutting surface meets the second characteristic requirement; otherwise, the cutting surface does not meet the second characteristic requirement; the discriminating module determines that the cutting surface does not meet the second characteristic requirement In the case of the femoral head model and the femoral body model, the model cutting module cuts the femoral head model and the femoral body model based on the femur model to regenerate the cutting surface.
根据一种优选的实施方式,所述判别模块还用于判别所述切割面是否满足第三特征需求;如果所述切割面与所述髋关节面形成的空间体的空间距离满足距离阈值,则所述切割面满足所述第三特征需求;否则,所述切割面不满所述第三特征需求;在所述判别模块判别出所述切割面不满足所述第三特征需求的情况下,所述模型切割模块基于所述股骨模型对所述股骨头部模型和股骨体模型进行切割以重新生成切割面。According to a preferred embodiment, the judging module is further configured to judge whether the cutting plane meets the third feature requirement; if the spatial distance between the cutting plane and the spatial body formed by the hip joint surface meets the distance threshold, then The cutting surface meets the third feature requirement; otherwise, the cutting surface does not meet the third feature requirement; in the case that the discriminating module determines that the cutting surface does not meet the third feature requirement, the The model cutting module cuts the femoral head model and the femoral body model based on the femur model to regenerate cutting surfaces.
根据一种优选的实施方式,所述系统还包括回转轴构建模块,用于将至少两个符合所述第一特征需求的旋转轴生成一个回转轴,其中,在所述股骨头部模型绕所述回轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动。According to a preferred embodiment, the system further comprises a rotational axis building module for generating a rotational axis from at least two rotational axes that meet the requirements of the first feature, wherein the femoral head model is formed around the femoral head model. In the case of rotation of the gyrus axis, the femoral head model can be rotated within the hip joint surface in such a manner that the partial region of the femoral head is moved out of the weight-bearing area.
根据一种优选的实施方式,所述内翻轴构建模块能够在所述回转轴构建模块构建了所述回转轴的情况下基于所述回转轴构建内翻轴而不基于旋转轴构建内翻轴,以使得所述内翻轴能够与所述股骨头部模型的旋转中心轴配套。According to a preferred embodiment, the varus axis building module is capable of constructing the varus axis based on the rotational axis and not based on the rotational axis if the rotational axis building module constructs the rotational axis , so that the varus shaft can be matched with the rotation center shaft of the femoral head model.
根据一种优选的实施方式,所述旋转轴构建模块能够按照如下方式构建所述旋转轴:计算所述切割面与所述股骨头部模型的点集合,并基于所述点集合形成的平面图形生成所述平面图形的重心作为第一标志点;以所述股骨头部模型的重心作为第二标志点;构建同时穿过所述第一标志点和所述第二标志点的直线作为所述旋转轴。According to a preferred embodiment, the rotation axis building module can construct the rotation axis as follows: calculating a point set of the cutting plane and the femoral head model, and based on the plane figure formed by the point set Generating the center of gravity of the plane figure as the first marker point; taking the center of gravity of the femoral head model as the second marker point; constructing a straight line passing through the first marker point and the second marker point as the axis of rotation.
根据一种优选的实施方式,所述内翻轴构建模块按照如下方式构建所述内翻轴:生成能包围完整的股骨头部模型表面的方向包围盒;提取构成方向包围盒的所有线段中的最短线段;对最短线段进行归一化处理得到第一线段;将第一线段移动到所述第一标志点处得到第二线段;对所述旋转轴与所述第二线段进行叉积运算得到归一化的第三线段;求第三线段与完整的股骨模型的表面的交点作为第三标志点;和将第二线段沿第三线段平移到第三标志点处作为内翻轴。According to a preferred embodiment, the varus axis building module constructs the varus axis as follows: generating a directional bounding box that can enclose the surface of the complete femoral head model; extracting all the line segments constituting the directional bounding box The shortest line segment; the first line segment is obtained by normalizing the shortest line segment; the second line segment is obtained by moving the first line segment to the first mark point; the cross product is performed on the rotation axis and the second line segment The operation obtains the normalized third line segment; finds the intersection of the third line segment and the surface of the complete femur model as the third marker point; and translates the second line segment along the third line segment to the third marker point as the varus axis.
根据一种优选的实施方式,本发明还提供了一种基于股骨旋转轴与内翻轴分析模型的构建方法,包括:基于检查数据构建股骨-髋臼模型;呈现所述股骨-髋臼模型;基于所述股骨-髋臼模型生成髋臼模型和股骨模型;所述方法还包括:基于所述股骨模型生成能够用于分割股骨头部模型和股骨体模型的切割面,基于所述股骨头部模型和所述切割面构建旋转轴,其中,在所述股骨头部模型绕所述旋转轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动,基于所述股骨头部模型的骨表面和所述旋转轴构建内翻轴,其中,在所述股骨头部模型绕所述内翻轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移除负重区的方式与所述髋关节面生成倾角。According to a preferred embodiment, the present invention also provides a method for constructing an analytical model based on a femoral rotation axis and a varus axis, comprising: constructing a femur-acetabular model based on inspection data; presenting the femoral-acetabular model; generating an acetabular model and a femoral bone model based on the femur-acetabular model; the method further includes: generating a cutting plane that can be used to segment a femoral head model and a femoral body model based on the femoral head model, based on the femoral head The model and the cutting plane establish an axis of rotation, wherein the femoral head model can be positioned at the hip joint in such a way as to move the femoral head sub-region out of the weight-bearing zone when the femoral head model is rotated about the rotation axis. In-plane rotation, a varus axis is constructed based on the bone surface of the femoral head model and the axis of rotation, wherein the femoral head model is rotated about the varus axis An inclination can be created with the hip articular surface in such a way that the partial region of the femoral head is removed from the weight bearing zone.
根据一种优选的实施方式,所述方法还包括:判别所述旋转轴是否满足第一特征需求;如果所述旋转轴与所述切割面的法向量的角度小于第一角度阈值,则所述旋转轴满足所述第一特征需求;否则,所述旋转轴不满足所述第一特征需求;在判别出所述旋转轴不满足所述第一特征需求的情况下,基于所述股骨模型对所述股骨头部模型和股骨体模型进行切割以重新生成切割面。According to a preferred embodiment, the method further includes: judging whether the rotation axis satisfies the first characteristic requirement; if the angle between the rotation axis and the normal vector of the cutting plane is less than a first angle threshold, then the The rotation axis satisfies the first feature requirement; otherwise, the rotation axis does not meet the first feature requirement; when it is determined that the rotation axis does not meet the first feature requirement, based on the femur model The femoral head model and femoral body model are cut to regenerate cut surfaces.
基于本发明提供的技术方案,本发明相比较于现有技术方案至少具有如下优势:Based on the technical solution provided by the present invention, the present invention has at least the following advantages compared with the prior art solution:
(1)在本发明中,在股骨头部模型绕旋转轴转动的情况下,股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动。这种模拟操作,与实际手术的操作是相近或者相同的,即在股骨头部分区域移除负重区后,股骨头仍然保持在髋关节面内;这有利于医生在手术前对于手术方案的确定是具有重要价值的;(1) In the present invention, when the femoral head model is rotated about the rotation axis, the femoral head model can be rotated within the hip joint surface so as to move the femoral head partial region out of the weight-bearing area. This simulated operation is similar to or the same as the actual operation, that is, after the weight-bearing area is removed from the femoral head, the femoral head remains in the hip joint surface; this is helpful for doctors to determine the surgical plan before the operation. is of great value;
(2)能够明确地限定出旋转轴和切割面的关系位置,从而能够提高旋转轴的构建精度,这对于手术方案的确定是十分有利的;(2) The relationship between the rotation axis and the cutting surface can be clearly defined, so that the construction accuracy of the rotation axis can be improved, which is very beneficial for the determination of the surgical plan;
(3)本发明能够对旋转轴的位置进行约束,使得股骨头部模型在绕旋转轴模型转动时,股骨头模型能够在髋关节面的窝形构造体内转动;其次,这种判别是基于髋关节面形成的空间体的重心线的角度为基础进行的,有利于提高股骨头部分区域移出负重区所需旋转角度的准确度。(3) The present invention can constrain the position of the rotation axis, so that when the femoral head model rotates around the rotation axis model, the femoral head model can rotate in the socket-shaped structure of the hip joint surface; secondly, this judgment is based on the hip Based on the angle of the center of gravity of the space body formed by the articular surface, it is beneficial to improve the accuracy of the rotation angle required for the partial area of the femoral head to move out of the weight-bearing area.
附图说明Description of drawings
图1是本发明提供的一种系统的模块之间的优选的关系示意图;和1 is a schematic diagram of a preferred relationship between modules of a system provided by the present invention; and
图2是本发明提供的一种方法的一种优选的流程示意图。Fig. 2 is a preferred schematic flow chart of a method provided by the present invention.
附图标记列表List of reference signs
1:模型构建模块 5:旋转轴构建模块1: Model building block 5: Rotation axis building block
2:图像显示模块 6:内翻轴构建模块2: Image Display Module 6: Inverted Shaft Building Module
3:模型分割模块 7:判别模块3: Model segmentation module 7: Discrimination module
4:模型切割模块 8:旋转轴拟合模块4: Model cutting module 8: Rotation axis fitting module
具体实施方式Detailed ways
下面结合附图1-2进行详细说明。A detailed description will be given below in conjunction with accompanying drawings 1-2.
实施例1Example 1
本实施例公开一种基于股骨旋转轴与内翻轴的分析模型的构建系统。This embodiment discloses a construction system based on the analysis model of the femoral rotation axis and the varus axis.
该构建系统中,涉及到重要术语如下,In this build system, the important terms involved are as follows,
股骨-髋臼模型:由股骨-髋臼部位真实检查数据构建的三维模型。Femur-acetabular model: A 3D model constructed from actual inspection data of the femoral-acetabular site.
股骨模型:由股骨-髋臼模型分离出来的股骨部位的模型。Femur Model: A model of the femoral part isolated from the femoral-acetabular model.
髋臼模型:由股骨-髋臼模型分离出来的髋臼部位的模型Acetabular model: a model of the acetabular part isolated from the femur-acetabular model
切割面:用于将股骨头部模型和股骨体模型分离的虚拟平面。Cutting plane: The virtual plane used to separate the femoral head model and the femoral body model.
旋转轴:是本实施例中构建一条虚拟轴线。Rotation axis: a virtual axis is constructed in this embodiment.
内翻轴:是本实施例中构建的另一条虚拟轴线。Inversion axis: is another virtual axis constructed in this embodiment.
股骨头部模型:包括股骨头模型和股骨颈模型。股骨头模型对应人体的股骨头。股骨颈模型对应人体的股骨颈。Femoral head model: including femoral head model and femoral neck model. The femoral head model corresponds to the femoral head of the human body. The femoral neck model corresponds to the femoral neck of the human body.
股骨体模型:对应人体的股骨体。Femoral body model: corresponds to the femoral body of the human body.
髋关节面:股骨头与髋臼通过髋关节配合,即髋臼上具有一个容纳股骨头的窝形构造。髋关节面是与股骨头表面相互吻合的窝形面。股骨头/股骨头部模型能够在髋关节面内绕旋转轴转动。Hip articular surface: The femoral head and the acetabulum cooperate through the hip joint, that is, the acetabulum has a socket-shaped structure that accommodates the femoral head. The hip articular surface is a socket-shaped surface that fits with the surface of the femoral head. The femoral head/femoral head model can be rotated about the axis of rotation within the hip joint surface.
本实施例提供一种基于股骨旋转轴与内翻轴的分析模型的构建系统,如图1所示,包括模型构建模块1、图像显示模块2、模型切割模块4、旋转轴构建模块5和内翻轴构建模块6。模型构建模块1和图像显示模块2数据交互。优选地,本系统还包括模型分割模块3。图像显示模块2和模型分割模块3数据交互。模型切割模块4和模型分割模块3数据交互。模型切割模块4分别与旋转轴构建模块5和内翻轴构建模块6数据交互。具体地,This embodiment provides a construction system based on the analysis model of the femoral rotation axis and varus axis. As shown in FIG. 1 , it includes a model construction module 1, an image display module 2, a model cutting module 4, a rotation
模型构建模块1,基于检查数据构建股骨-髋臼模型。检查数据至少来源于两种不同的扫面成像方式。例如,检查数据主要来源于CT、MRI等医疗扫描数据。当然,检查数据不仅限于CT和MRI扫描数据,还可以是PET等扫描数据。CT,是指计算机断层造影或者计算机断层扫描。MRI是Nuclear Magnetic Resonance Imaging的缩写,是指磁共振断层造影或者磁共振断层扫描。本实施例通过配准融合方法将CT采集的股骨-髋臼的计算机断层扫面数据和MR采集的股骨-髋臼的计算机断层扫描数据构建股骨-髋臼模型。医学图像配准与融合技术主要包括两方面的技术:医学图像配准和图像融合技术。医学图像配准旨在确定两个图像之间的几何变换关系,通过这一变换可以把一个图像中的任意一点的坐标变换到另一个图像中对应位置的坐标。医学图像配准主要包括三个步骤:提取图像特征、根据特征向量对并确定几何变换和执行变换。图像融合技术主要分为:以图像像素为基础的融合方法和以图像特征为基础的融合方法。Model building module 1, builds a femur-acetabular model based on the inspection data. Examination data are derived from at least two different scanning imaging modalities. For example, inspection data mainly comes from medical scan data such as CT and MRI. Of course, the examination data is not limited to CT and MRI scan data, but may also be scan data such as PET. CT stands for computed tomography or computed tomography. MRI is the abbreviation of Nuclear Magnetic Resonance Imaging, which refers to magnetic resonance tomography or magnetic resonance tomography. In this embodiment, a femur-acetabular model is constructed by using the registration fusion method to construct the femur-acetabular computed tomography scan data acquired by CT and the femoral-acetabular computed tomography scan data acquired by MR. Medical image registration and fusion technology mainly includes two technologies: medical image registration and image fusion technology. Medical image registration aims to determine the geometric transformation relationship between two images, through which the coordinates of any point in one image can be transformed to the coordinates of the corresponding position in the other image. Medical image registration mainly consists of three steps: extracting image features, determining geometric transformations based on feature vector pairs, and performing transformations. Image fusion technology is mainly divided into: fusion methods based on image pixels and fusion methods based on image features.
图像显示模块2,用于呈现股骨-髋臼模型。例如,将MRI扫描数据和CT数据导入构建系统中并基于模型构建模块1生成股骨-髋臼模型,图像显示模块2可以将其显示出来。图像显示模块2可以是液晶显示屏、VR眼镜和计算机屏幕等显示装置。Image display module 2, used to present the femur-acetabular model. For example, MRI scan data and CT data are imported into the construction system and a femur-acetabular model is generated based on the model construction module 1, and the image display module 2 can display it. The image display module 2 may be a display device such as a liquid crystal display, VR glasses, and a computer screen.
模型分割模块3,能够基于股骨-髋臼模型生成髋臼模型和股骨模型。髋臼和股骨是两个相对独立的骨组织。本实施例中,髋臼模型和股骨模型的生成可以按照如下方法生成:有界面生成、轮廓提取和轮廓图分割。有界面生成用于根据股骨-髋臼模型的三角形网络数据对该模型中包含的所有三角形进行处理,以生成适用于对该三维模型进行分割的有界面。轮廓提取用于通过有有界面提取三维模型的轮廓图。轮廓图分割,用于根据生成的有界面信息以及三维模型的顶点邻接图的信息,将提取出的轮廓图分割成满足预定条件的一个子图或者至少两个相互不重叠的子图。顶点邻接图可以通过以下方式构建:根据三维模型的三角形网络数据,以三维模型的顶点为节点,在被一个或多个三角形所共有的每两个顶点之间添加一条边,从而构建顶点邻接图。预定条件包括以下几种中的任意一个:子图的顶点数量小于预定的第一阈值;子图含有的有界面数量小于预定的第二阈值;子图含有的所有有界平面中不存在候选分割面;或子图中不存在可将与该子图对应的轮廓图分割为至少两个子图的候选分割面。又如,模型分割模块3能够基于布尔运算将股骨-髋臼模型分割为髋臼模型和股骨模型。再如,模型分割模块3还可以采用医学图像三维分割技术——基于结构的分割技术。基于结构的分割技术是通过检测边缘或体数据中的边界面进行图像分割,其包括:通过区分不同属性检测出边缘点;将这些边缘点组合成连续的轮廓,将感兴趣的区域的体素从其他体素中分离出来。The model segmentation module 3 can generate an acetabular model and a femur model based on the femur-acetabular model. The acetabulum and femur are two relatively independent bone tissues. In this embodiment, the generation of the acetabular model and the femur model can be generated according to the following methods: there are interface generation, contour extraction and contour map segmentation. The interface generation is used to process all triangles contained in the femoral-acetabular model according to the triangle network data of the model to generate interfaces suitable for segmentation of the 3D model. Contour extraction is used to extract the contour map of the 3D model through an interface. The contour map segmentation is used to divide the extracted contour map into one sub-image or at least two non-overlapping sub-images that satisfy a predetermined condition according to the generated interface information and the information of the vertex adjacency graph of the three-dimensional model. A vertex adjacency graph can be constructed in the following ways: According to the triangle network data of the 3D model, taking the vertices of the 3D model as nodes, adding an edge between every two vertices shared by one or more triangles, thereby constructing a vertex adjacency graph . The predetermined condition includes any one of the following: the number of vertices of the subgraph is less than a predetermined first threshold; the number of bounded surfaces contained in the subgraph is less than a predetermined second threshold; no candidate segmentation exists in all bounded planes contained in the subgraph face; or there is no candidate dividing face in the subgraph that can divide the contour graph corresponding to the subgraph into at least two subgraphs. For another example, the model segmentation module 3 can segment the femur-acetabular model into an acetabular model and a femur model based on Boolean operations. For another example, the model segmentation module 3 may also adopt a three-dimensional segmentation technology of medical images—a structure-based segmentation technology. Structure-based segmentation technology performs image segmentation by detecting edges or boundary surfaces in volume data, which includes: detecting edge points by distinguishing different attributes; separated from other voxels.
模型切割模块4,基于股骨模型生成能够用于分割股骨头部模型和股骨体模型的切割面。切割面还能够作为用于旋转轴生成的基础面。优选地,切割面是平面。例如,切割面的生成可以按照基于结构的分割技术生成。首先,操作人员可以选择股骨头部上的一个点作为分割面的一个点,然后模型切割模块4即与结构特征和该点生成分割面。又如,切割面的生成还可以基于统计学方法生成,例如基于阈值分割算法、马尔科夫随机场模型、聚类算法或者分类器算法。再如,切割面的生成还可以基于混合分割算法,例如区域增长算法、分裂合并算法或图谱引导算法。The model cutting module 4 generates a cutting plane that can be used to divide the femoral head model and the femoral body model based on the femur model. Cut planes can also be used as base planes for rotation axis generation. Preferably, the cutting plane is a plane. For example, the generation of cut planes can be generated according to structure-based segmentation techniques. First, the operator can select a point on the femoral head as a point of the segmentation plane, and then the model cutting module 4 generates the segmentation plane with the structural features and this point. For another example, the generation of the cutting plane can also be generated based on a statistical method, for example, based on a threshold segmentation algorithm, a Markov random field model, a clustering algorithm or a classifier algorithm. For another example, the generation of the cutting plane may also be based on a hybrid segmentation algorithm, such as a region growing algorithm, a split-merge algorithm, or a map-guided algorithm.
旋转轴构建模块5,其能够基于股骨头部模型和切割面构建旋转轴。在本发明中,旋转轴是一条虚拟轴线,主要使得股骨头部模型能够模拟股骨头部旋转,用于外科医生观察部分区域移除负重区域所需的最大旋转角度。该部分区域是由医生选定的股骨头坏死区域。医生在构建系统中操作该模型时能够对角度进行测量,以判定其最大旋转角度。在真实截骨术过程中,需要将股骨头部从股骨体上截取下来,然后将股骨头部分区域的大部分或者全部旋转出负重区。由于股骨头部周边存在大量的血管和神经,在实际手术过程中,股骨头能转动的角度有限,其不能完全的旋转一周。据临床数据表明,股骨头逆时针转动仅能旋转90°,逆时针旋转仅能旋转130°。在本发明中,在股骨头部模型绕旋转轴转动的情况下,股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动。这种模拟操作,与实际手术的操作是相近或者相同的,即在股骨头部分区域移除负重区后,股骨头仍然保持在髋关节面内。这有利于医生在手术前对于手术方案的确定是具有重要价值的。优选地,旋转轴构建模块5可以按照如下方式生成旋转轴。在切割面的位置关系符合上述第二特征需求和第三特征需求的情况下,旋转轴构建模块5计算切割面与股骨头部模型的点集合。此时,切割面与股骨头部模型的点集合形成一个封闭平面图形。旋转轴构建模块5将该封闭平面图形的重心作为构建旋转轴的第一标志点。切割面与股骨头部模型形成了一个以空间体。旋转轴构建模块5将该空间体的重心作为第二标志点。第一标志点与第二标志点构成的直线记为旋转轴。Rotation
内翻轴构建模块6,其能够依据股骨头部模型的骨表面和旋转轴构建内翻轴。其中,在股骨头部模型绕内翻轴转动的情况下,股骨头部模型能够以将股骨头部分区域移除负重区的方式与髋关节面生成倾角。优选地,内翻轴构建模块6可以按照如下方式生成内翻轴:在旋转轴满足第一特征需求的情况下,内翻轴构建模块6生成能包围完整的股骨头部模型表面的方向包围盒。提取构成方向包围盒的所有线段中的最短线段。对最短线段进行归一化处理得到第一线段。将第一线段移动到第一标志点处得到第二线段。对旋转轴与第二线段进行叉积运算得到归一化的第三线段。求第三线段与完整的股骨模型的表面的交点作为第三标志点。将第二线段沿第三线段平移到第三标志点处作为内翻轴。The varus axis building block 6, which can construct the varus axis from the bone surface of the femoral head model and the axis of rotation. Wherein, when the femoral head model is rotated around the varus axis, the femoral head model can generate an inclination angle with the hip joint surface in a manner of removing a partial area of the femoral head from the weight-bearing area. Preferably, the varus axis building module 6 can generate the varus axis as follows: in the case that the rotation axis satisfies the first feature requirement, the varus axis building module 6 generates a direction bounding box that can enclose the complete femoral head model surface . Extracts the shortest line segment among all the line segments that make up the bounding box of the direction. The first line segment is obtained by normalizing the shortest line segment. Move the first line segment to the first marker point to get the second line segment. A normalized third line segment is obtained by performing a cross product operation on the rotation axis and the second line segment. Find the intersection of the third line segment and the surface of the complete femur model as the third marker point. Translate the second line segment along the third line segment to the third mark point as the inversion axis.
旋转轴和内翻轴的构建精度对于临床方案的确定是至关重要的。The construction accuracy of the rotational and varus axes is critical for clinical protocol determination.
第一方面,在该构建系统中,为了能够通过该系统模拟出股骨头部分区域大部分旋转移除负重区所需要的角度,为手术决策提供精确的数据支撑,旋转轴的构建精度是至关重要的。为了能够有效地控制构建的旋转轴能够与实际手术过程中的旋转轴的误差,优选地,本发明提供的系统还包括了判别模块7,主要用于判别旋转轴的构建精度。结合图2,J3:判别模块7用于判别旋转轴是否满足第一特征需求。如果旋转轴与切割面的法向量的角度小于第一角度阈值,则旋转轴满足第一特征需求。否则,旋转轴不满足第一特征需求。第一角度阈值可以根据部分区域的面积来决定。部分区域的面积与第一角度阈值呈负相关。例如,部分区域占股骨头的10%,则第一角度阈值为1.5°。部分区域占股骨头的12%,则第一角度阈值为1.4°。部分区域的面积能够根据检查数据由医生确定,然后输入该构建系统。由于在一种优选的实施方式过程中,旋转轴是要穿过切割面的重心以及切割面与股骨头部模型围成的三维空间体的重心。因而,在判别模块7判别出旋转轴不满足第一特征需求的情况下,模型切割模块4基于股骨模型对股骨头部模型和股骨体模型进行切割以重新生成切割面。通过这种设置:能够明确地限定出旋转轴和切割面的关系位置,从而能够提高旋转轴的构建精度,这对于手术方案的确定是十分有利的。First, in this construction system, in order to be able to simulate the angle required for most of the rotation of the femoral head to remove the weight-bearing area through the system, and to provide accurate data support for surgical decision-making, the construction accuracy of the rotation axis is crucial. important. In order to effectively control the error between the constructed rotation axis and the rotation axis during the actual operation, preferably, the system provided by the present invention further includes a
第二方面,旋转轴构建的准确度与切割面的构建精度和准确度是正相关的。如果切割面构造具有明显的偏差,这对于将股骨头部分区域移除负重区所需要的旋转角度的判定是十分不利的。因此,优选地,判别模块7在切割面生成之后就需要对切割面的位置关系进行判定。具体地,结合图2,J1:判别模块7还用于判别切割面是否满足第二特征需求。第二特征需求是用于判定切割面相对于髋关节面的位置是否满足旋转轴的构建需求。优选地,第二特征需求是用于判定切割面相对于髋关节面的窝形构造形成的空间体的位置是否满足旋转轴的构建需求。如果切割面的法向量与髋关节面形成的空间体的重心线的角度小于第二角度阈值,则切割面满足第二特征需求。否则,切割面不满第二特征需求。第二角度阈值同样与部分区域的面积呈负相关。部分区域面积越大,第二角度阈值越小。部分区域面积越小,第二角度阈值越大。第二角度阈值还与个体股骨-髋臼的结构特征有关。因此,具体的第二角度阈值的由医生设定。在判别模块7判别出切割面不满足第二特征需求的情况下,模型切割模块4基于股骨模型对股骨头部模型和股骨体模型进行切割以重新生成切割面。从对切割面的约束条件可以看出,该约束是对旋转轴相对髋关节面的位置进行约束。通过这种判别设置,本发明能够对旋转轴的位置进行约束,使得股骨头部模型在绕旋转轴模型转动时,股骨头模型能够在髋关节面的窝形构造体内转动;其次,这种判别是基于髋关节面形成的空间体的重心线的角度为基础进行的,有利于提高股骨头部分区域移出负重区所需旋转角度的准确度。In the second aspect, the accuracy of the rotation axis construction is positively correlated with the construction accuracy and accuracy of the cutting plane. If the cutting plane configuration has obvious deviation, it is very disadvantageous for the determination of the rotation angle required to remove the partial area of the femoral head from the load-bearing area. Therefore, preferably, the
第三方面,切割面距离髋关节面的距离对于手术方案的确定是十分重要的。在实际手术过程中,截骨术是将股骨颈从股骨体上的大转子截下后进行旋转,将股骨头上的部分区域移出负重区。基于此,优选地,切割面是股骨颈模型与股骨体模型的分界平面。因此,切割面与髋关节面的距离关系是有必要进行精确要求的。优选地,结合图2,J2:判别模块7还用于判别切割面是否满足第三特征需求。第三特征需求是对切割面与髋关节面之间的距离进行精确限定。切割面不能距离髋关节面太远,以使得切割面不与股骨颈模型有交点。切割面不能距离髋关节面太远,以使得切割面不予股骨颈模型有交点或者与股骨头模型有交点。每一个人体的股骨颈与大转子交接的面距离髋关节的距离都是不一样的。因此,该距离阈值与人体的股骨-髋臼结构相关。医生能够根据实际情况对距离阈值进行设定。如果切割面与髋关节面形成的空间体的空间距离满足距离阈值,则切割面满足第三特征需求。否则,切割面不满第三特征需求。在判别模块7判别出切割面不满足第三特征需求的情况下,模型切割模块4基于股骨模型对股骨头部模型和股骨体模型进行切割以重新生成切割面。例如,距离阈值可以是的2cm~12cm。Thirdly, the distance between the cutting surface and the hip joint surface is very important for the determination of the surgical plan. In actual surgery, osteotomy is to cut the femoral neck from the greater trochanter on the femoral body and then rotate it to move part of the femoral head out of the weight-bearing area. Based on this, preferably, the cutting plane is the boundary plane between the femoral neck model and the femoral body model. Therefore, the distance relationship between the cutting surface and the hip joint surface is necessary to be precise. Preferably, with reference to Fig. 2, J2: the
实施例2Example 2
本实施例是实施例1的进一步补充或者技术方案的进一步优化。This embodiment is a further supplement of Embodiment 1 or a further optimization of the technical solution.
由于旋转轴是在切割面的基础上生成的。由于满足第一特征需求、第二特征需求和第三特征需求的切割面是无穷多个。从而,旋转轴也是无穷多个的。为此,本实施例为了进一步地提高旋转轴的构建精度。Since the axis of rotation is generated on the basis of the cutting plane. Since there are infinitely many cutting planes that satisfy the first feature requirement, the second feature requirement and the third feature requirement. Therefore, the number of rotation axes is also infinite. Therefore, in this embodiment, the construction accuracy of the rotating shaft is further improved.
模型切割模块4至少会生成两个符合第一特征需求、第二特征需求和第三特征需求的切割面。此时,旋转轴构建模块5构建的旋转轴符合第一特征需求。优选地,旋转轴构建模块5能够生成至少两个与切割面彼此对应的旋转轴。此时,结合图2,步骤S7:旋转轴拟合模块8会基于至少两个旋转轴生成一个回转轴。The model cutting module 4 generates at least two cutting surfaces that meet the first feature requirement, the second feature requirement and the third feature requirement. At this time, the rotation axis constructed by the rotation
优选地,旋转轴拟合模块8按照如下方式基于至少两个旋转轴生成一个回转轴:Preferably, the rotation axis
标号,将至少两个旋转轴标号:L1、L2……Ln。Label, label at least two axes of rotation: L1 , L2 . . . Ln .
构造直线,选取至少两条能够同时与直线L1、L2……Ln有交点的平行直线。To construct a straight line, select at least two parallel straight lines that can intersect with the straight lines L1 , L2 . . . Ln at the same time.
交点求取:求取第一条直线与旋转轴L1、L2……Ln的交点,然后求取同一条直线上交点坐标的平均值,求得一个平均交点。依次类推,求得每一条直线与旋转轴L1、L2……Ln形成的平均交点。将这些交点进行连线获得回转轴。股骨头部模型绕该回转轴旋转。Obtaining the intersection point: Obtain the intersection point between the first straight line and the rotation axis L1 , L2 ......Ln , and then obtain the average of the coordinates of the intersection points on the same straight line to obtain an average intersection point. By analogy, the average point of intersection formed by each straight line and the rotation axes L1 , L2 . . . Ln is obtained. Connect these intersection points to obtain the rotary axis. The femoral head model rotates about this axis of rotation.
例如,旋转轴构建模块5生成了两根旋转轴L1和L2。空间上存在两条平行直线S1和S2。S1与L1具有交点A(x1,y1,z1)。S1与L2具有交点B(x2,y2,z2)。S1与L1、L2形成的平均交点为同理,S2与L1具有交点C(x3,y3,z3)。S2与L2具有交点D(x4,y4,z4)。S2与L1、L2形成的平均交点为那么回转轴以为EF。For example, the rotational
优选地,旋转轴拟合模块8按照如下方式基于至少两个旋转轴生成一个回转轴:Preferably, the rotation axis
首先,对旋转轴在直角坐标系下进行空间几何建模描述:First, the spatial geometric modeling description of the rotation axis in the Cartesian coordinate system:
L1:a1x+b1y+c1z=d1;L2:a2x+b2y+c2z=d2;···Ln:anx+bny+cnz=dn;L1 : a1 x+b1 y+c1 z=d1 ; L2 : a2 x+b2 y+c2 z=d2 ;...Ln : anx +b ny+ cn z=dn ;
其次,将旋转轴的数学方程进行叠加运算可得数学模型:Secondly, the mathematical model can be obtained by superposing the mathematical equation of the rotation axis:
该叠加后的数学模型为一条直线,将该直线作为回转轴的空间方程。The superimposed mathematical model is a straight line, and the straight line is used as the space equation of the rotation axis.
优选地,股骨头部模型的旋转中心轴应该与内翻轴是相互配套使用的。因此,在股骨头模型的旋转中心轴是由若干个旋转轴构建的回转轴时,内翻轴构建模块6能够基于所述回转轴构建内翻轴而不基于旋转轴构建内翻轴。Preferably, the rotation center axis of the femoral head model should be matched with the varus axis. Therefore, when the central axis of rotation of the femoral head model is a rotational axis constructed by several rotational axes, the varus axis building module 6 can construct the varus axis based on the rotational axis without constructing the varus axis based on the rotational axis.
实施例3Example 3
本实施例提供了一种基于股骨旋转轴与内翻轴的分析模型的构建方法,包括:This embodiment provides a method for constructing an analysis model based on a femoral rotation axis and a varus axis, including:
S1,基于检查数据构建股骨-髋臼模型。S1, build a femur-acetabular model based on the inspection data.
S2,呈现所述股骨-髋臼模型。S2, the femur-acetabular model is presented.
S3,所述股骨-髋臼模型生成髋臼模型和股骨模型。S3, the femur-acetabular model generates an acetabular model and a femur model.
S4,基于所述股骨模型生成能够用于分割股骨头部模型和股骨体模型的切割面。S4, based on the femur model, generate a cutting plane that can be used to segment the femoral head model and the femoral body model.
S5,基于所述股骨头部模型和所述切割面构建旋转轴,其中,在所述股骨头部模型绕所述旋转轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移出负重区的方式在髋关节面内转动。S5 , constructing a rotation axis based on the femoral head model and the cutting plane, wherein, when the femoral head model rotates around the rotation axis, the femoral head model can transform the femoral head partial area Move out of the weight-bearing zone in a way that rotates within the hip joint surface.
S6,基于所述股骨头部模型的骨表面和所述旋转轴构建内翻轴,其中,在所述股骨头部模型绕所述内翻轴转动的情况下,所述股骨头部模型能够以将股骨头部分区域移除负重区的方式与所述髋关节面生成倾角。S6, constructing a varus axis based on the bone surface of the femoral head model and the rotation axis, wherein, when the femoral head model is rotated around the varus axis, the femoral head model can be The manner in which the partial area of the femoral head is removed from the weight-bearing area creates an inclination with the hip articular surface.
优选地,在步骤S5中,需要对旋转轴的构建精度进行判定。判别所述旋转轴是否满足第一特征需求。如果所述旋转轴与所述切割面的法向量的角度小于第一角度阈值,则所述旋转轴满足所述第一特征需求。否则,所述旋转轴不满足所述第一特征需求。在判别出所述旋转轴不满足所述第一特征需求的情况下,基于所述股骨模型对所述股骨头部模型和股骨体模型进行切割以重新生成切割面。Preferably, in step S5, the construction accuracy of the rotation axis needs to be determined. It is judged whether the rotation axis satisfies the first characteristic requirement. If the angle between the rotation axis and the normal vector of the cutting plane is smaller than a first angle threshold, the rotation axis satisfies the first feature requirement. Otherwise, the rotation axis does not meet the first feature requirement. In the case that it is determined that the rotation axis does not meet the first feature requirement, the femoral head model and the femoral body model are cut based on the femoral model to regenerate a cut surface.
如本发明中所用的词语“模块”描述任一种硬件、软件或软硬件组合,其能够执行与“模块”相关联的功能。例如,模型构建模块1、模型分割模块3、模型切割模块4、旋转轴构建模块5、内翻轴构建模块6判别模块7和旋转轴拟合模块8均是一种数据运算处理器以及配置了执行其各自模块中具体步骤的算法。例如,模型构建模块1、模型分割模块3、模型切割模块4、旋转轴构建模块5、内翻轴构建模块6判别模块7和旋转轴拟合模块8可以是处理器、服务器或者专用集成芯片中的一种或几种。The term "module" as used in this disclosure describes any hardware, software or combination of hardware and software capable of performing the functions associated with the "module". For example, the model building module 1, the model segmentation module 3, the model cutting module 4, the rotation
需要注意的是,上述具体实施例是示例性的,本领域技术人员可以在本发明公开内容的启发下想出各种解决方案,而这些解决方案也都属于本发明的公开范围并落入本发明的保护范围之内。本领域技术人员应该明白,本发明说明书及其附图均为说明性而并非构成对权利要求的限制。本发明的保护范围由权利要求及其等同物限定。It should be noted that the above-mentioned specific embodiments are exemplary, and those skilled in the art can come up with various solutions inspired by the disclosure of the present invention, and these solutions also belong to the disclosure scope of the present invention and fall within the scope of the present invention. within the scope of protection of the invention. It should be understood by those skilled in the art that the description of the present invention and the accompanying drawings are illustrative rather than limiting to the claims. The protection scope of the present invention is defined by the claims and their equivalents.
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011325087.1ACN112451090B (en) | 2019-11-06 | 2019-11-06 | An Analytical Model Construction System |
| CN201911080296.1ACN110811829B (en) | 2019-11-06 | 2019-11-06 | A method and system for constructing an analytical model based on femoral rotation axis and varus axis |
| CN202011324990.6ACN112472292B (en) | 2019-11-06 | 2019-11-06 | Construction method and system for orthopedic surgery analysis model |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911080296.1ACN110811829B (en) | 2019-11-06 | 2019-11-06 | A method and system for constructing an analytical model based on femoral rotation axis and varus axis |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011325087.1ADivisionCN112451090B (en) | 2019-11-06 | 2019-11-06 | An Analytical Model Construction System |
| CN202011324990.6ADivisionCN112472292B (en) | 2019-11-06 | 2019-11-06 | Construction method and system for orthopedic surgery analysis model |
| Publication Number | Publication Date |
|---|---|
| CN110811829Atrue CN110811829A (en) | 2020-02-21 |
| CN110811829B CN110811829B (en) | 2020-12-08 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011325087.1AActiveCN112451090B (en) | 2019-11-06 | 2019-11-06 | An Analytical Model Construction System |
| CN201911080296.1AActiveCN110811829B (en) | 2019-11-06 | 2019-11-06 | A method and system for constructing an analytical model based on femoral rotation axis and varus axis |
| CN202011324990.6AActiveCN112472292B (en) | 2019-11-06 | 2019-11-06 | Construction method and system for orthopedic surgery analysis model |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011325087.1AActiveCN112451090B (en) | 2019-11-06 | 2019-11-06 | An Analytical Model Construction System |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011324990.6AActiveCN112472292B (en) | 2019-11-06 | 2019-11-06 | Construction method and system for orthopedic surgery analysis model |
| Country | Link |
|---|---|
| CN (3) | CN112451090B (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115272615A (en)* | 2022-08-10 | 2022-11-01 | 新乡学院 | Medical semantic segmentation method for femoral point cloud model and application of medical semantic segmentation method |
| CN117442395A (en)* | 2023-09-06 | 2024-01-26 | 北京长木谷医疗科技股份有限公司 | Method, device and equipment for acquiring femoral head rotation center based on clustering algorithm |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114758074A (en)* | 2022-04-14 | 2022-07-15 | 北京斫月医疗科技有限公司 | Method and device for judging distance between acetabulum and femoral head |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002061688A2 (en)* | 2001-01-29 | 2002-08-08 | The Acrobot Company Limited | Modelling for surgery |
| US20030176860A1 (en)* | 2002-03-18 | 2003-09-18 | Fuji Photo Film Co., Ltd. | Operation aid system |
| CA2427590C (en)* | 1999-11-01 | 2008-07-08 | Arthrovision, Inc. | Evaluating disease progression using magnetic resonance imaging |
| CN101711695A (en)* | 2008-09-30 | 2010-05-26 | 德普伊产品公司 | Customized patient-specific acetabular orthopaedic surgical instrument and method of use and fabrication |
| CN101969869A (en)* | 2007-12-06 | 2011-02-09 | 史密夫和内修有限公司 | Systems and methods for determining the mechanical axis of a femur |
| WO2013158416A2 (en)* | 2012-04-20 | 2013-10-24 | Bausch & Lomb Incorporated | System and method for creating a customized anatomical model of an eye |
| CN107274389A (en)* | 2017-05-25 | 2017-10-20 | 中国科学院苏州生物医学工程技术研究所 | Femur and Acetabular dissection parameter acquiring method based on CT three-dimensional series images |
| CN108348193A (en)* | 2015-08-05 | 2018-07-31 | F·贝特纳 | A Fluoroscopy-Based Technique for Measuring Intraoperative Cup Anteversion |
| CN108697471A (en)* | 2016-03-02 | 2018-10-23 | 思想外科有限公司 | The arthroplasty of automation is planned |
| CN108921832A (en)* | 2018-06-26 | 2018-11-30 | 陈卫衡 | Femoral head image analysis method, device, server and medium |
| CN109496143A (en)* | 2016-05-27 | 2019-03-19 | 马科外科公司 | It is preoperative to plan and be registrated in the related art to surgery systems |
| US20190262081A1 (en)* | 2013-03-14 | 2019-08-29 | Biomet Manufacturing, Llc | Method for implanting a hip prosthesis and related system |
| CN110215281A (en)* | 2019-06-11 | 2019-09-10 | 北京和华瑞博科技有限公司 | A kind of femur or shin bone method for registering and device based on total knee replacement |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7468075B2 (en)* | 2001-05-25 | 2008-12-23 | Conformis, Inc. | Methods and compositions for articular repair |
| JP2006509609A (en)* | 2002-10-04 | 2006-03-23 | オルトソフト インコーポレイテッド | Computer-aided hip replacement surgery |
| CN1187026C (en)* | 2003-05-29 | 2005-02-02 | 上海交通大学 | Fully automatic femur reference axis determining method |
| US7780672B2 (en)* | 2006-02-27 | 2010-08-24 | Biomet Manufacturing Corp. | Femoral adjustment device and associated method |
| KR101973101B1 (en)* | 2009-05-29 | 2019-04-26 | 스미스 앤드 네퓨, 인크. | Methods and apparatus for performing knee arthroplasty |
| CN104622559B (en)* | 2014-12-23 | 2017-03-15 | 河海大学常州校区 | A kind of construction method of parametrization femur template |
| CN104614200A (en)* | 2015-02-10 | 2015-05-13 | 唐佩福 | Computer CT and 3D printing technology-based method for precisely sampling any cancellous bone area in human femoral head |
| CN105741354A (en)* | 2016-03-22 | 2016-07-06 | 陈继营 | Preoperative hip joint deformity bone model manufacturing method based on 3D printing technology |
| CN106388978B (en)* | 2016-09-13 | 2018-07-20 | 中南大学湘雅医院 | Preparation method of acetabulum side model and guide plate based on three-dimensional reconstruction technology |
| US20180360540A1 (en)* | 2017-06-16 | 2018-12-20 | Episurf Ip-Management Ab | System and method for creating a decision support material indicating damage to an anatomical joint |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2427590C (en)* | 1999-11-01 | 2008-07-08 | Arthrovision, Inc. | Evaluating disease progression using magnetic resonance imaging |
| WO2002061688A2 (en)* | 2001-01-29 | 2002-08-08 | The Acrobot Company Limited | Modelling for surgery |
| US20030176860A1 (en)* | 2002-03-18 | 2003-09-18 | Fuji Photo Film Co., Ltd. | Operation aid system |
| CN101969869A (en)* | 2007-12-06 | 2011-02-09 | 史密夫和内修有限公司 | Systems and methods for determining the mechanical axis of a femur |
| CN101711695A (en)* | 2008-09-30 | 2010-05-26 | 德普伊产品公司 | Customized patient-specific acetabular orthopaedic surgical instrument and method of use and fabrication |
| WO2013158416A2 (en)* | 2012-04-20 | 2013-10-24 | Bausch & Lomb Incorporated | System and method for creating a customized anatomical model of an eye |
| US20190262081A1 (en)* | 2013-03-14 | 2019-08-29 | Biomet Manufacturing, Llc | Method for implanting a hip prosthesis and related system |
| CN108348193A (en)* | 2015-08-05 | 2018-07-31 | F·贝特纳 | A Fluoroscopy-Based Technique for Measuring Intraoperative Cup Anteversion |
| CN108697471A (en)* | 2016-03-02 | 2018-10-23 | 思想外科有限公司 | The arthroplasty of automation is planned |
| CN109496143A (en)* | 2016-05-27 | 2019-03-19 | 马科外科公司 | It is preoperative to plan and be registrated in the related art to surgery systems |
| CN107274389A (en)* | 2017-05-25 | 2017-10-20 | 中国科学院苏州生物医学工程技术研究所 | Femur and Acetabular dissection parameter acquiring method based on CT three-dimensional series images |
| CN108921832A (en)* | 2018-06-26 | 2018-11-30 | 陈卫衡 | Femoral head image analysis method, device, server and medium |
| CN110215281A (en)* | 2019-06-11 | 2019-09-10 | 北京和华瑞博科技有限公司 | A kind of femur or shin bone method for registering and device based on total knee replacement |
| Title |
|---|
| 季文婷: "人工全髋关节理论活动度的参数化动态仿真", 《系统仿真学报》* |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115272615A (en)* | 2022-08-10 | 2022-11-01 | 新乡学院 | Medical semantic segmentation method for femoral point cloud model and application of medical semantic segmentation method |
| CN115272615B (en)* | 2022-08-10 | 2023-09-08 | 新乡学院 | A medical semantic segmentation method for femur point cloud model and its application |
| CN117442395A (en)* | 2023-09-06 | 2024-01-26 | 北京长木谷医疗科技股份有限公司 | Method, device and equipment for acquiring femoral head rotation center based on clustering algorithm |
| Publication number | Publication date |
|---|---|
| CN112451090A (en) | 2021-03-09 |
| CN112451090B (en) | 2022-02-18 |
| CN112472292A (en) | 2021-03-12 |
| CN110811829B (en) | 2020-12-08 |
| CN112472292B (en) | 2022-01-18 |
| Publication | Publication Date | Title |
|---|---|---|
| AU2019396551C1 (en) | Systems and methods for preoperative planning and postoperative analysis of surgical procedures | |
| CN107016666B (en) | Image Post-processing Method for 3D Reconstruction of Pelvic CT Based on Coordinate System | |
| CN110811829B (en) | A method and system for constructing an analytical model based on femoral rotation axis and varus axis | |
| JP2016532475A (en) | Method for optimal visualization of bone morphological regions of interest in X-ray images | |
| CN117078840A (en) | Automatic quantitative calculation method for three-dimensional modeling of hip joint based on CT image | |
| Bartels et al. | Computed tomography-based joint locations affect calculation of joint moments during gait when compared to scaling approaches | |
| CN110680353B (en) | Analysis device for femoral head integrity rate | |
| US12062183B2 (en) | Closed surface fitting for segmentation of orthopedic medical image data | |
| Gomes et al. | Patient-specific modelling in orthopedics: from image to surgery | |
| CN110811830B (en) | Analysis device and method based on femoral model | |
| Chen | Reconstruction individual three-dimensional model of fractured long bone based on feature points | |
| Chang et al. | A pre‐operative approach of range of motion simulation and verification for femoroacetabular impingement | |
| van Loon et al. | Automatic identification of radius and ulna bone landmarks on 3D virtual models | |
| Mahfouz et al. | Reconstruction of 3D Patient-Specific Bone Models From Biplanar X-Ray Images Utilizing Morphometric Measurements. | |
| CN116797634A (en) | Image registration method for three-dimensional broken bone registration and splicing oriented to anatomical reduction | |
| Handels et al. | Simulation of hip operations and design of custom-made endoprostheses using virtual reality techniques | |
| Song et al. | Computer-aided modeling and morphological analysis of hip joint | |
| Stindel et al. | Bone morphing: 3D reconstruction without pre-or intraoperative imaging-concept and applications | |
| Valenti | Statistical shape models based 2D/3D registration methods for knee orthopaedic surgery | |
| CN119856980B (en) | Navigation method, device and system for orthopedic surgery, electronic equipment and storage medium | |
| Albers et al. | Computer assisted diagnosis and treatment planning of femoroacetabular impingement (FAI) | |
| Guezou-Philippe et al. | Statistical shape modeling to determine the anterior pelvic plane for total hip arthroplasty | |
| Negrillo-Cárdenas et al. | Automatic Detection of Distal Humerus Features: First Steps. | |
| Atmani et al. | From medical data to simple virtual mock-up of scapulo-humeral joint | |
| Zhang et al. | Video-Driven Comprehensive 3D Hip Joint Motion Model for FAI Auxiliary Diagnosis |
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information | Inventor after:Cheng Hui Inventor after:Luo Dianzhong Inventor after:Zhang Hong Inventor after:Xie Zhiheng Inventor after:Ming Wenhua Inventor after:Huang Yupu Inventor before:Cheng Hui Inventor before:Luo Dianzhong Inventor before:Zhang Hong Inventor before:Xie Zhiheng Inventor before:Ming Wenhua Inventor before:Huang Yupu | |
| GR01 | Patent grant | ||
| GR01 | Patent grant |