CN109662716A - 软骨厚度测量方法、装置、计算机设备和存储介质 - Google Patents
软骨厚度测量方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN109662716A CN109662716A CN201811566467.7A CN201811566467A CN109662716A CN 109662716 A CN109662716 A CN 109662716A CN 201811566467 A CN201811566467 A CN 201811566467A CN 109662716 A CN109662716 A CN 109662716A
- Authority
- CN
- China
- Prior art keywords
- point
- cartilage
- interest
- location information
- edge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 210000000845 cartilage Anatomy 0.000 title claims abstract description 269
- 238000000691 measurement method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000004590 computer program Methods 0.000 claims description 17
- 210000000988 bone and bone Anatomy 0.000 description 23
- 238000004364 calculation method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 4
- 241000905957 Channa melasoma Species 0.000 description 2
- 210000001188 articular cartilage Anatomy 0.000 description 2
- 210000004394 hip joint Anatomy 0.000 description 2
- 210000003127 knee Anatomy 0.000 description 2
- 201000008482 osteoarthritis Diseases 0.000 description 2
- 208000025747 Rheumatic disease Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1079—Measuring physical dimensions, e.g. size of the entire body or parts thereof using optical or photographic means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Dentistry (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Analysis (AREA)
Abstract
本申请涉及一种软骨厚度测量方法、装置、计算机设备和存储介质。所述方法包括：根据3D软骨图像数据，获取软骨表面轮廓信息；根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。采用本方法能够对软骨表面任意一点的厚度进行测量。
Description
技术领域
本申请涉及医学影像技术领域，特别是涉及一种软骨厚度测量方法、装置、计算机设备和存储介质。
背景技术
骨关节炎是最常见的风湿性疾病，它以极高的致残率，给患者带来的沉重精神负担和经济负担，由此日益受到人们的重视。通过测量关节软骨的厚度和观测厚度的变化，可以辅助骨关节炎的诊断、病情的监测、治疗效果的评估等。
目前，对软骨的测量方法，一般是针对软骨的几个角度或者层面进行测量，而无法测量软骨表面任意一点的软骨厚度。
发明内容
基于此，有必要针对上述技术问题，提供一种能够测量软骨表面所有点的软骨厚度的软骨厚度方法、装置、计算机设备和存储介质。
一种软骨厚度测量方法，所述方法包括：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述根据3D软骨图像数据，获取软骨表面轮廓信息，包括：根据3D软骨图像数据，获得3D软骨模型信息，所述3D软骨模型信息包括所述3D软骨模型中目标点的位置信息；根据所述3D软骨模型中目标点的位置信息，判断所述目标点邻域内是否存在背景点；如果所述目标点邻域内存在背景点，则所述目标点为构成软骨表面轮廓的点，所述目标点的位置信息构成所述软骨表面轮廓中点的位置信息。
在其中一个实施例中，所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二表面的最短距离，得到所述第一感兴趣点位置的软骨厚度。
所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一表面中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨表面的第一切面；
根据所述第一切面在所述第一表面的第一法线，得到所述第一法线与所述第二表面的第一交点；
计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。
一种软骨厚度测量方法，所述方法包括：
获取软骨图像的2D软骨切面图像数据；
根据所述2D软骨切面图像数据，获取软骨边缘信息，所述软骨边缘信息中携带构成所述边缘的任意一点的位置信息；
根据所述软骨边缘信息，识别所述软骨边缘中的第一边缘和第二边缘；
根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二边缘的最短距离，得到所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一边缘中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨边缘的第一切线；
根据所述第一切线在所述第一边缘的第一法线，得到所述第一法线与所述第二边缘的第一交点；
计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。
一种软骨厚度测量装置，所述装置包括：
表面轮廓获取模块，用于根据3D软骨图像数据，获取软骨表面轮廓信息；
表面分类模块，用于根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
软骨厚度计算模块，用于根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
上述软骨厚度测量方法、装置、计算机设备和存储介质，通过将软骨表面从软骨图像中分离出来，再根据软骨其中一个表面的点的位置信息，计算所述点的位置的软骨厚度，能够对软骨表面任意一点的位置的软骨厚度进行计算，为后续软骨的诊断治疗提供了依据。
附图说明
图1为一个实施例中软骨厚度测量方法的流程示意图；
图2为一个实施例中获取软骨表面轮廓的方法的流程示意图；
图3为一个具体实施例中软骨表面轮廓的示意图；
图4为另一个实施例中软骨厚度测量方法的流程示意图；
图5为一个具体实施例中根据法线计算软骨厚度的示意图；
图6为一个实施例中软骨厚度测量装置的结构框图；
图7为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
在一个实施例中，如图1所示，提供了一种软骨厚度测量方法，包括以下步骤：
步骤S110，根据3D软骨图像数据，获取软骨表面轮廓信息。
其中，3D软骨图像数据可为MR医学影像或CT医学影像，所述软骨包括但不限于膝关节软骨、髋关节软骨、指关节软骨。拍摄3D软骨图像的过程中，其它骨头或组织的图像也会采集，因此，需要根据软骨的密度或形状特征，将软骨区域从3D软骨图像中分离出来。
步骤S120，根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面。
其中，软骨一般呈扁平状，第一表面、第二表面可为扁平状软骨的上下表面。具体的，第一表面、第二表面是与其它骨相连的面，如图3所示，月牙形状区域表面为软骨表面轮廓，第一表面、第二表面分别位于月牙形状区域两侧。
所述识别所述软骨表面中的第一表面和第二表面包括采用图像分割法、微分法、极值法、机器学习、曲面拟合。
步骤S130，根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度。
其中，所述第一感兴趣点位于所述第一表面，软骨表面轮廓信息是已知的，根据第一感兴趣点的位置信息，可以通过几何计算公式，计算出第一感兴趣点位置的软骨厚度。
步骤S140，根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
同理，根据第二感兴趣点的位置信息和软骨表面轮廓信息，通过几何计算公式，计算出第二感兴趣点位置的软骨厚度。
本实施例中，步骤S130和步骤S140的执行顺序不分先后，也可以仅执行步骤S130或步骤S140，或者两个步骤都执行。
在其中一个实施例中，在步骤S130和步骤S140中，所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二表面的最短距离，得到所述第一感兴趣点位置的软骨厚度。同理，根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度，包括：根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点到所述第一表面的最短距离，得到所述第二感兴趣点位置的软骨厚度。
其中，第一感兴趣点和第二感兴趣点的位置已知，同时，第一表面和第二表面中任意一点的位置信息已知，直接计算点到表面的最短距离。本实施例中，最短距离即为软骨厚度。
所述计算所述第一感兴趣点到所述第二表面的最短距离，计算方法包括遍历法、拉格朗日乘数法、跟踪法、切平面法，此最短距离对应的第二表面中的点称为第一感兴趣点对应的软骨厚度端点。同样可得第二感兴趣点到所述第一表面的最短距离，及其对应的软骨厚度端点。在其中一个实施例中，在步骤S130中，所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：根据所述第一表面中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨表面的第一切面；根据所述第一切面在所述第一表面的第一法线，得到所述第一法线与所述第二表面的第一交点；计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。其中，所述第一感兴趣点到所述第一交点的距离，即为所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，在步骤S140中，根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度，包括：根据所述第二表面中第二感兴趣点的位置信息，获取所述第二感兴趣点在所述软骨表面的第二切面；根据所述第二切面在所述第二表面的第二法线，得到所述第二法线与所述第一表面的第二交点；计算所述第二感兴趣点到所述第二交点的距离，得到所述第二感兴趣点位置的软骨厚度。其中，所述第二感兴趣点到所述第二交点的距离，即为所述第二感兴趣点位置的软骨厚度。
上述软骨厚度测量方法中，通过将软骨表面从软骨图像中分离出来，再根据软骨其中一个表面的点的位置信息，计算所述点的位置的软骨厚度，能够对软骨表面任意一点的位置的软骨厚度进行计算，为后续软骨的诊断治疗提供了依据。
在其中一个实施例中，在步骤S110中，如图2所示，所述根据3D软骨图像数据，获取软骨表面轮廓信息，包括步骤：
步骤S111，根据3D软骨图像数据，获得3D软骨模型信息，所述3D软骨模型信息包括所述3D软骨模型中目标点的位置信息。
其中，3D软骨模型信息为只包括软骨部分的图像数据。3D软骨模型可以还原软骨的形状和位置。
步骤S112，根据所述3D软骨模型中目标点的位置信息，判断所述目标点邻域内是否存在背景点。
其中，3D软骨模型由多个软骨图像点组成，位于3D软骨模型内部的点，其邻域内只包含软骨图像点，位于3D软骨模型表面的点，其邻域会存在其它骨头或组织的图像点。背景点即为除软骨以外其它骨头或组织的图像点。
步骤S113，如果所述目标点邻域内存在背景点，则所述目标点为构成软骨表面轮廓的点，所述目标点的位置信息构成所述软骨表面轮廓中点的位置信息。
在一个实施例中，如图4所示，提供了一种软骨厚度测量方法，包括以下步骤：
步骤S210，获取软骨图像的2D软骨切面图像数据。
其中，2D软骨切面图像数据只包括软骨图像的二维信息，所能展示也只是软骨平面图像。例如，2D软骨切面图像是垂直于软骨厚度方向拍摄，以完全展示软骨的厚度。
步骤S220，根据所述2D软骨切面图像数据，获取软骨边缘信息，所述软骨边缘信息中携带构成所述边缘的任意一点的位置信息。
其中，所述软骨包括但不限于膝关节软骨、髋关节软骨、指关节软骨。由于，软骨的密度或形状特征与周围的骨头和组织不同，可以将软骨边缘信息提取出来。
步骤S230，根据所述软骨边缘信息，识别所述软骨边缘中的第一边缘和第二边缘。
具体的，所述软骨边缘围成月牙形状，第一边缘和第二边缘位于月牙形状两侧。
步骤S240，根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度。
其中，所述第一感兴趣点位于所述第一边缘，软骨边缘信息是已知的，根据第一感兴趣点的位置信息，可以通过几何计算公式，计算出第一感兴趣点位置的软骨厚度。
步骤S250，根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
同理，根据第二感兴趣点的位置信息和软骨边缘信息，通过几何计算公式，计算出第二感兴趣点位置的软骨厚度。
本实施例中，步骤S240和步骤S250的执行顺序不分先后，也可以仅执行步骤S240或步骤S250，或者两个步骤都执行。
在其中一个实施例中，在步骤S240和步骤S250中，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二边缘的最短距离，得到所述第一感兴趣点位置的软骨厚度。同理，所述根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度，包括：根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点到所述第一边缘的最短距离，得到所述第二感兴趣点位置的软骨厚度。
其中，第一感兴趣点和第二感兴趣点的位置已知，同时，第一边缘和第二边缘中任意一点的位置信息已知，直接计算点到表面的最短距离。本实施例中，最短距离即为软骨厚度。
所述计算所述第一感兴趣点到所述第二边缘的最短距离，计算方法包括遍历法、拉格朗日乘数法、切线法，此最短距离对应的第二边缘中的点称为第一感兴趣点对应的软骨厚度端点。同样可得第二感兴趣点到所述第一边缘的最短距离，及其对应的软骨厚度端点。
在其中一个实施例中，在步骤S240中，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：根据所述第一边缘中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨边缘的第一切线；根据所述第一切线在所述第一边缘的第一法线，得到所述第一法线与所述第二边缘的第一交点；计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。其中，所述第一感兴趣点到所述第一交点的距离，即为所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，在步骤S250中，所述根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度，包括：根据所述第二边缘中第二感兴趣点的位置信息，获取所述第二感兴趣点在所述软骨边缘的第二切线；根据所述第二切线在所述第二边缘的第二法线，得到所述第二法线与所述第一边缘的第二交点；计算所述第二感兴趣点到所述第二交点的距离，得到所述第二感兴趣点位置的软骨厚度。其中，所述第二感兴趣点到所述第二交点的距离，即为所述第二感兴趣点位置的软骨厚度。
具体的，如图5所示，月牙形状为软骨区域，第一边缘和第二边缘位于月牙形状两侧，第二边缘中一点B做垂直于第二边缘的法线MN，法线MN与第一边缘相交于点A，计算点B到点A的距离，AB之间的距离即为点B位置的软骨厚度。
上述软骨厚度测量方法中，通过将软骨边缘分为第一边缘和第二边缘，再根据第一边缘和第二边缘中的点的位置信息，计算所述点的位置的软骨厚度，能够对软骨表面任意一点的位置的软骨厚度进行计算，为后续软骨的诊断治疗提供了依据。
应该理解的是，虽然图1、2、4的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1、2、4中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中，如图6所示，提供了一种软骨厚度测量装置，包括：表面轮廓获取模块310、表面分类模块320和软骨厚度计算模块330，其中：
表面轮廓获取模块310，用于根据3D软骨图像数据，获取软骨表面轮廓信息。
表面分类模块320，用于根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面。
软骨厚度计算模块330，用于根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块330，具体用于根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二表面的最短距离，得到所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块330，具体用于根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点到所述第一表面的最短距离，得到所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块330，包括：第一切面获取单元，用于根据所述第一表面中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨表面的第一切面；第一交点获取单元，用于根据所述第一切面在所述第一表面的第一法线，得到所述第一法线与所述第二表面的第一交点；厚度计算单元，用于计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。其中，所述第一感兴趣点到所述第一交点的距离，即为所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块330，包括：第二切面获取单元，用于根据所述第二表面中第二感兴趣点的位置信息，获取所述第二感兴趣点在所述软骨表面的第二切面；第二交点获取单元，用于根据所述第二切面在所述第二表面的第二法线，得到所述第二法线与所述第一表面的第二交点；厚度计算单元，用于计算所述第二感兴趣点到所述第二交点的距离，得到所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述表面轮廓获取模块310，包括：软骨模型获取单元，用于根据3D软骨图像数据，获得3D软骨模型信息，所述3D软骨模型信息包括所述3D软骨模型中目标点的位置信息；背景点判断单元，用于根据所述3D软骨模型中目标点的位置信息，判断所述目标点邻域内是否存在背景点，如果所述目标点邻域内存在背景点，则所述目标点为构成软骨表面轮廓的点，所述目标点的位置信息构成所述软骨表面轮廓中点的位置信息。
在一个实施例中，提供了一种软骨厚度测量装置，包括：
切面图像获取模块，用于获取软骨图像的2D软骨切面图像数据。
软骨边缘获取模块，用于根据所述2D软骨切面图像数据，获取软骨边缘信息，所述软骨边缘信息中携带构成所述边缘的任意一点的位置信息。
边缘分类模块，用于根据所述软骨边缘信息，识别所述软骨边缘中的第一边缘和第二边缘。
软骨厚度计算模块，用于根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块，具体用于根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二边缘的最短距离，得到所述第一感兴趣点位置的软骨厚度；和/或根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点到所述第一边缘的最短距离，得到所述第二感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块包括：第一切线获取单元，用于根据所述第一边缘中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨边缘的第一切线；第一交点获取单元，用于根据所述第一切线在所述第一边缘的第一法线，得到所述第一法线与所述第二边缘的第一交点；距离计算单元，用于计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。其中，所述第一感兴趣点到所述第一交点的距离，即为所述第一感兴趣点位置的软骨厚度。
在其中一个实施例中，所述软骨厚度计算模块包括：第二切线获取单元，用于根据所述第二边缘中第二感兴趣点的位置信息，获取所述第二感兴趣点在所述软骨边缘的第二切线；第二交底获取单元，用于根据所述第二切线在所述第二边缘的第二法线，得到所述第二法线与所述第一边缘的第二交点；距离计算单元，用于计算所述第二感兴趣点到所述第二交点的距离，得到所述第二感兴趣点位置的软骨厚度。其中，所述第二感兴趣点到所述第二交点的距离，即为所述第二感兴趣点位置的软骨厚度。
关于软骨厚度测量装置的具体限定可以参见上文中对于软骨厚度测量方法的限定，在此不再赘述。上述软骨厚度测量装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图7所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储软骨图像数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种软骨厚度测量方法。
本领域技术人员可以理解，图7中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种软骨厚度测量方法，其特征在于，所述方法包括：
根据3D软骨图像数据，获取软骨表面轮廓信息；
根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
2.根据权利要求1所述的方法，其特征在于，所述根据3D软骨图像数据，获取软骨表面轮廓信息，包括：
根据3D软骨图像数据，获得3D软骨模型信息，所述3D软骨模型信息包括所述3D软骨模型中目标点的位置信息；
根据所述3D软骨模型中目标点的位置信息，判断所述目标点邻域内是否存在背景点；
如果所述目标点邻域内存在背景点，则所述目标点为构成软骨表面轮廓的点，所述目标点的位置信息构成所述软骨表面轮廓中点的位置信息。
3.根据权利要求1所述的方法，其特征在于，所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二表面的最短距离，得到所述第一感兴趣点位置的软骨厚度。
4.根据权利要求1所述的方法，其特征在于，所述根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一表面中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨表面的第一切面；
根据所述第一切面在所述第一表面的第一法线，得到所述第一法线与所述第二表面的第一交点；
计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。
5.一种软骨厚度测量方法，其特征在于，所述方法包括：
获取软骨图像的2D软骨切面图像数据；
根据所述2D软骨切面图像数据，获取软骨边缘信息，所述软骨边缘信息中携带构成所述边缘的任意一点的位置信息；
根据所述软骨边缘信息，识别所述软骨边缘中的第一边缘和第二边缘；
根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或
根据所述第二边缘中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
6.根据权利要求5所述的方法，其特征在于，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点到所述第二边缘的最短距离，得到所述第一感兴趣点位置的软骨厚度。
7.根据权利要求5所述的方法，其特征在于，所述根据所述第一边缘中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度，包括：
根据所述第一边缘中第一感兴趣点的位置信息，获取所述第一感兴趣点在所述软骨边缘的第一切线；
根据所述第一切线在所述第一边缘的第一法线，得到所述第一法线与所述第二边缘的第一交点；
计算所述第一感兴趣点到所述第一交点的距离，得到所述第一感兴趣点位置的软骨厚度。
8.一种软骨厚度测量装置，其特征在于，所述装置包括：
表面轮廓获取模块，用于根据3D软骨图像数据，获取软骨表面轮廓信息；
表面分类模块，用于根据所述软骨表面轮廓信息，识别所述软骨表面中的第一表面和第二表面；
软骨厚度计算模块，用于根据所述第一表面中第一感兴趣点的位置信息，计算所述第一感兴趣点位置的软骨厚度；和/或根据所述第二表面中第二感兴趣点的位置信息，计算所述第二感兴趣点位置的软骨厚度。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811566467.7A CN109662716B (zh) | 2018-12-19 | 2018-12-19 | 软骨厚度测量方法、装置、计算机设备和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811566467.7A CN109662716B (zh) | 2018-12-19 | 2018-12-19 | 软骨厚度测量方法、装置、计算机设备和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109662716A true CN109662716A (zh) | 2019-04-23 |
| CN109662716B CN109662716B (zh) | 2021-10-22 |
Family
ID=66144125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811566467.7A Active CN109662716B (zh) | 2018-12-19 | 2018-12-19 | 软骨厚度测量方法、装置、计算机设备和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109662716B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110827285A (zh) * | 2019-11-15 | 2020-02-21 | 上海联影智能医疗科技有限公司 | 软骨厚度检测方法、装置、计算机设备和可读存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020177770A1 (en) * | 1998-09-14 | 2002-11-28 | Philipp Lang | Assessing the condition of a joint and assessing cartilage loss |
| CN1735373A (zh) * | 2003-01-07 | 2006-02-15 | 成像治疗仪股份有限公司 | 预测肌和骨骼疾病的方法 |
| CN1780594A (zh) * | 2002-11-07 | 2006-05-31 | 康复米斯公司 | 用于确定半月板的大小和形状以及设计治疗的方法 |
| CN102438559A (zh) * | 2009-02-25 | 2012-05-02 | 穆罕默德·拉什万·马赫福兹 | 定制矫形植入物及相关方法 |
| CN105361883A (zh) * | 2014-08-22 | 2016-03-02 | 方学伟 | 膝关节置换三维空间下肢生物力线的确定方法 |
-
2018
- 2018-12-19 CN CN201811566467.7A patent/CN109662716B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020177770A1 (en) * | 1998-09-14 | 2002-11-28 | Philipp Lang | Assessing the condition of a joint and assessing cartilage loss |
| CN1780594A (zh) * | 2002-11-07 | 2006-05-31 | 康复米斯公司 | 用于确定半月板的大小和形状以及设计治疗的方法 |
| CN1735373A (zh) * | 2003-01-07 | 2006-02-15 | 成像治疗仪股份有限公司 | 预测肌和骨骼疾病的方法 |
| CN102438559A (zh) * | 2009-02-25 | 2012-05-02 | 穆罕默德·拉什万·马赫福兹 | 定制矫形植入物及相关方法 |
| CN105361883A (zh) * | 2014-08-22 | 2016-03-02 | 方学伟 | 膝关节置换三维空间下肢生物力线的确定方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110827285A (zh) * | 2019-11-15 | 2020-02-21 | 上海联影智能医疗科技有限公司 | 软骨厚度检测方法、装置、计算机设备和可读存储介质 |
| CN110827285B (zh) * | 2019-11-15 | 2022-07-26 | 上海联影智能医疗科技有限公司 | 软骨厚度检测方法、装置、计算机设备和可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109662716B (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111325739B (zh) | 肺部病灶检测的方法及装置，和图像检测模型的训练方法 | |
| CN109767461B (zh) | 医学影像配准方法、装置、计算机设备和存储介质 | |
| CN111080573B (zh) | 肋骨图像检测方法、计算机设备和存储介质 | |
| JP3932303B2 (ja) | 臓器動態の定量化方法、装置、臓器位置の予測方法、装置、放射線照射方法、装置及び臓器異常検出装置 | |
| Kabus et al. | B-spline registration of 3D images with Levenberg-Marquardt optimization | |
| US20060188139A1 (en) | Tree structure based 2D to 3D registration | |
| CN110599526A (zh) | 图像配准方法、计算机设备和存储介质 | |
| US20190130572A1 (en) | Registration method and system for non-rigid multi-modal medical image | |
| Kim et al. | Automation of spine curve assessment in frontal radiographs using deep learning of vertebral-tilt vector | |
| US20220036575A1 (en) | Method for measuring volume of organ by using artificial neural network, and apparatus therefor | |
| CN110751187A (zh) | 异常区域图像生成网络的训练方法和相关产品 | |
| CN113935943A (zh) | 颅内动脉瘤识别检测的方法、装置、计算机设备和存储介质 | |
| Jiang et al. | B‐spline Based Multi‐organ Detection in Magnetic Resonance Imaging | |
| US11138736B2 (en) | Information processing apparatus and information processing method | |
| JP2016209399A (ja) | 画像処理装置及び方法、並びにコンピュータプログラム | |
| CN109662716A (zh) | 软骨厚度测量方法、装置、计算机设备和存储介质 | |
| Motta et al. | Fundus image transformation revisited: Towards determining more accurate registrations | |
| CN115552458A (zh) | 测量医学图像中肿瘤体积的变化 | |
| CN105551042A (zh) | 一种扫描床标记点位置的确定方法和装置 | |
| Xu et al. | A symmetric 4D registration algorithm for respiratory motion modeling | |
| CN118229975A (zh) | 一种乳腺mri的乳腺、肿瘤及腺体分割方法、系统及介质 | |
| CN112488982A (zh) | 一种超声图像检测方法及装置 | |
| Neeraja et al. | A review on automatic cephalometric landmark identification using artificial intelligence techniques | |
| CN109598751B (zh) | 一种医学影像图片处理的方法、设备及装置 | |
| Pinheiro et al. | Deep homography based localization on videos of endoscopic capsules |
Legal Events
| 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 | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |