CN107016666B - 基于坐标体系的骨盆ct三维重建的图像后处理方法 - Google Patents
基于坐标体系的骨盆ct三维重建的图像后处理方法 Download PDFInfo
- Publication number
- CN107016666B CN107016666B CN201710114304.4A CN201710114304A CN107016666B CN 107016666 B CN107016666 B CN 107016666B CN 201710114304 A CN201710114304 A CN 201710114304A CN 107016666 B CN107016666 B CN 107016666B
- Authority
- CN
- China
- Prior art keywords
- pelvis
- coordinate system
- dimensional
- dimensional model
- axis
- 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.)
- Expired - Fee Related
Links
- 210000004197 pelvis Anatomy 0.000 title claims abstract description 98
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000012805 post-processing Methods 0.000 title claims abstract description 18
- 238000006073 displacement reaction Methods 0.000 claims abstract description 60
- 206010017076 Fracture Diseases 0.000 claims abstract description 26
- 206010034246 Pelvic fractures Diseases 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000013178 mathematical model Methods 0.000 claims abstract description 5
- 210000000988 bone and bone Anatomy 0.000 claims description 58
- 230000005484 gravity Effects 0.000 claims description 49
- 208000010392 Bone Fractures Diseases 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000011218 segmentation Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 210000001694 thigh bone Anatomy 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000004445 quantitative analysis Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 230000000873 masking effect Effects 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 2
- 238000012067 mathematical method Methods 0.000 claims 1
- 238000005728 strengthening Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 abstract description 15
- 230000006378 damage Effects 0.000 description 9
- 208000027418 Wounds and injury Diseases 0.000 description 7
- 208000014674 injury Diseases 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004451 qualitative analysis Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 206010053206 Fracture displacement Diseases 0.000 description 3
- 208000011092 Hand injury Diseases 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004393 prognosis Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 208000035965 Postoperative Complications Diseases 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 210000001981 hip bone Anatomy 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000968 intestinal effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 230000002746 orthostatic effect Effects 0.000 description 1
- 210000001696 pelvic girdle Anatomy 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 210000004061 pubic symphysis Anatomy 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 230000002485 urinary effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开的基于坐标体系的骨盆CT三维重建的图像后处理方法,属于医用图像处理领域。包括以下步骤:101)采集病人骨盆CT数据:采集病人术前骨盆二维CT的病例数据;102)建立病例三维模型:利用病例数据,通过三维重建,得到病例的骨盆三维模型;103)在病例的骨盆三维模型中建立坐标系;104)获得骨盆骨折的距离移位和旋转移位:根据103)步骤中的坐标系中建立数学模型,精确计算骨折部投射在X、Y、Z轴的距离移位和旋转移位。本发明能帮助骨科医生方便直观、精确地了解患者骨盆骨折的移位方式及程度,并可为患者的术中复位进行精确指导。
Description
技术领域
本发明涉及一种医用图像处理方法,是一种基于坐标体系的骨盆CT三维重建的图像后处理方法,特别涉及一种骨盆骨折的移位分析及其移位程度的精确定量方法。
背景技术
骨盆骨折常由高能量损伤导致,不稳定的骨盆骨折常合并各种并发症,诸如泌尿、肠道及神经血管等的损伤,其致残率和死亡率都较高。近年来,随着电子导航设备、内固定系统以及手术技术的进步,骨盆骨折的治疗逐渐向微创化和精准化手术转变,如何依靠现有的影像资料,掌握骨盆骨折的移位方式和移位程度是手术复位成功、减少并发症和获得良好预后的关键。但由于骨盆解剖结构的深在和特殊性,损伤机制的复杂性,不稳定骨盆骨折后其移位方式也多种多样。传统的影像资料诸如骨盆X片,因伪影遮挡、体位受限等的影响,较难确定骨折的具体移位及旋转方向。而二维CT图像可以显示骨折移位的位置,但要求有较强的空间三维想象能力,同时在显示移位的距离和角度,尤其是旋转移位时表现欠佳,无法用于指导复位。目前常用的骨盆CT三维重建图像可以动态、清晰地显示骨盆骨折的位置,并且观察损伤时移位及旋转方向。但是,因为缺乏具体的标准参照,没有具体的移位数值,导致遵循“逆力学损伤机制”进行复位的过程中施力的具体力度与方向无法把控与监督,使得靠现有的影像图像技术进行闭合复位甚至是有限切开复位治疗骨盆仍然十分困难。
为了弥补影像资料目前骨盆骨折移位显示和骨折复位方面的不足,国内外医师进行了很多努力,并根据影像资料建立了很多定性分析,而在定量分析方面,国内外暂无相关文献报道。目前国内外骨科医师最常用的定性分型为Tile分型和Young-Burgess骨盆分型,是根据患者的受伤机制、解剖位置以及骨盆X片来制定的。在骨科临床工作中我们发现,这两种分型对于评价患者的骨折损伤情况、损伤机制,评估治疗方式、推测患者的病程转归及预后有重大作用。但对于骨折三维显示、尤其是不稳定骨折的精准复位指导帮助不够。
发明内容
本发明旨在提供一种骨盆CT三维图像的图像后处理及三维骨盆移位程度定量分析的建立方法,利用三维立体坐标体系对骨盆CT三维重建的图像进行后处理,建立几种操作简便、结构简单、重复性强、可准确显示和定量计算骨盆骨折移位程度和角度的骨盆移位分类方法。
为达到上述目的,本发明采用的技术方案如下:
首先建立一种均值骨盆CT三维坐标体系,包括以下步骤:
201)采集CT数据:采集N个健康成人骨盆二维CT的N组原始数据;
202)建立三维模型:利用N组原始数据,分别建立N个三维模型Ai,去除三维模型Ai的股骨和腰椎,得到N个骨盆的三维模型Bi;
203)装配三维模型:获取三维模型Bi的重心Bi,以重心Bi为原点,建立坐标系Bi;将N个三维模型Bi连同坐标系Bi一起重叠,得到基准坐标系B0和基准原点B0;
204)拟合对齐:利用最小二乘法原理,固定任意一个三维模型Bj为参考模型,将所有其他三维模型Bi和三维模型Bj的骶尾骨部位进行最佳拟合对齐,获得一个拟合对齐后的三维模型C;
205)重建均值骨盆模型:以基准坐标系B0的XY平面为起始平面沿着Z轴将三维模型C进行等层距平行双向切割,获得截面轮廓线,将获得的所有轮廓线统一装配,在每一个截面上均构建平面网格,记录Z轴位点坐标,并获取和记录网格与轮廓线交点的点云数据,通过数学运算中的均值计算来获取均值点云数据,形成一系列平行于XY面的均值离散坐标点(Xi,Yi),利用正交多项拟合曲线法排除偏离主点的点云,并将离散坐标点按照多项拟合曲线走向进行有次序地连接和封闭,将所有平行于XY平面的曲线通过点云逆向三维重建形成均值三维模型D;
206)建立均值骨盆三维坐标系D:获取三维模型D的重心D,以重心D为原点,建立坐标系D;
所述N≥5,i为1-N的任意自然数,j为1-N的任意自然数。
优选的,所述202)步骤中,建立三维模型Ai后,进行阈值分割,并予区域增长将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除三维模型Ai的股骨和腰椎。
优选的,所述203)步骤中,所述坐标系Bi面向骨盆正位水平面上过重心Bi指向右侧为X轴,经过重心Bi垂直于骨盆正位水平面向上为Z轴,经过重心Bi垂直于X、Z轴指向骶骨为Y轴。
进一步的,在进行203)步骤之前,先对三维模型Bi进行表面松弛、重划网格优化处理。
进一步的,所述204)步骤中,通过以下公式处理,
设骨盆三维模型Bi上一点P(x1,x2,x3),对应参照模型Bj上的P'(y1,y2,y3),
求解上述方程,可解的未知数aij的值,即可得到模型Bi和Bj的转换矩阵,完成三维模型的最佳拟合对齐。
优选的,所述205)步骤中,所述层距为1mm。
优选的,所述206)步骤中,所述坐标系D面向骨盆正位水平面上过重心D指向右侧为X轴,经过重心D垂直于水平面向上为Z轴,经过重心D垂直于X、Z轴指向骶骨为Y轴。
基于坐标体系的骨盆CT三维重建的图像后处理方法,包括以下步骤:
101)采集病人骨盆CT数据:采集病人术前骨盆二维CT的病例数据;
102)建立病例三维模型:利用病例数据,建立病例三维模型,去除病例三维模型的股骨和腰椎,得到病例的骨盆三维模型;
103)在病例的骨盆三维模型中选择均值骨盆重心D为原点建立笛卡尔三维坐标系D或选择髂前下棘连线中点E为原点建立笛卡尔三维坐标系E;
104)获得骨盆骨折的距离移位和旋转移位:根据103)步骤中的坐标系分别计算骨折部投射在X、Y、Z轴的距离移位和旋转移位。
优选的,所述102)步骤中,建立病例三维模型后,进行阈值分割,并予区域增长将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除病例三维模型的股骨和腰椎。
优选的,所述103)步骤中,以206)步骤中的重心D为原点,建立坐标系D。
优选的,所述103)步骤中,以髂前下棘连线中点为原点,建立坐标系E,所述坐标系E面向骨盆正位水平面上过髂前下棘连线中点E指向右侧为X轴,经过髂前下棘连线中点E垂直于骨盆正位水平面向上为Z轴,经过髂前下棘连线中点E垂直于X、Z轴指向骶骨为Y轴。
进一步的,所述104)步骤中,计算骨折的距离移位及角度移位:
方法一、通过模拟游标卡尺的方法即“两平面法”:即在建模软件中,建立两个相互平行且垂直于移位轴的平面,通过小心平移两个平面,使得两平面均与移位骨块断端两侧恰好相切,且相切的两点为相互复位点,此时通过参数计算的两平面之间的距离即为移位的距离。计算骨折的角度移位:即移位骨块和健侧骨块相切的平面与移位轴平面形成的夹角之差。
方法二、401)利用建模软件以重心D为中心,以移位轴为对称轴,在对侧建立健侧骨盆的虚拟镜像U;
402)通过阈值分割、编辑蒙版,分离出移位骨块A,扫描移位骨块上的几何特征,并将移位骨块W的空间信息以点集的形式记录Ai={(xi,yi,zi)}。将移位骨块A通过平移或者旋转使得与健侧镜像U恰好完全重合,得到B,采集此时的几何特征,以点集Bi={(xi,yi,zi)}记录。
403)建立数学模型,利用采集的点集信息,分别求取Ai和Bi的均分平面La和Lb,则骨块移位程度的定量分析,可以通过La和Lb之间的距离及旋转移位来计算。具体计算方法为:设La的方程为a1x+b1y+c1z+d1=0,Lb为a2x+b2y+c2z+d2=0;根据La和Lb分别均分Ai和Bi,可知点集Ai和Bi上所有点分别到两平面的距离的平方之和最短,可利用点到面的距离公式及偏导数进行求解,即:可知A和B上的任一点到La和Lb上的距离为:
方法三:501)同方法二得到移位骨块A,及可完全与镜像U完全重合的骨块B;502)利用Geomagic Studio软件执行特征选择-均分平面操作,则可较粗略地得到均分移位骨块的平面La和Lb的方程矩阵,继而求解。
本发明具有以下有益效果:
1、本发明的图像的数据来源为患者术前骨盆CT三维重建图像,该技术成熟、数据准确、获取容易,而采取的图像后处理方法中的定性及其中两种定量方法也较为简单,前期复杂的医学重建软件过程是为确定最佳原点及坐标系,实际运用过程中只需要将图像上建立我们确定的三轴坐标系即可完成骨折的三轴移位定性分析,用于指导术中的复位,尤其是闭合复位治疗骨盆骨折的治疗中。此方法简单易行,就算基层医院也可以实施,而且结果显示精确,符合目前精准化治疗的趋势,具有较大的推广价值。
2、本发明区别于传统的图像处理技术,将临床问题、放射影像图像、医学软件的配准技术、数学三维坐标系和精准医学有效地综合利用,在骨盆三维重建图片中将复杂的骨盆骨折分解成X、Y、Z三个坐标轴的移位,可将所有骨盆骨折都标准化一致化,做到操作简单、图像清晰、三维可视性和可重复性强。个体化地通过数学几何、偏导数的方法计算移位程度,计算定量结果准确。
3、创新性地将均值骨盆三维模型的重心骨盆原点建立坐标系。重心作为物体能够保持平衡的点,重心稳固不变是其能够维持稳定、受力平衡的重要因素。我们将重心作为原点,也作为骨盆骨折的移位参考点,参考骨盆原点及坐标轴位置,可对骨折骨块的移位情况进行可靠评估。重心也可作为骨盆术中的复位中心,这符合使骨盆环稳定平衡的复位标准,用于指导骨盆骨折的术中复位。因此,将重心定义为骨盆原点建立的三轴坐标系,可重复地应用于骨盆骨折倾斜度及旋转度的评估,以及骨盆骨折的移位程度的精确定量计算中。在骨盆移位的定性分析或较粗略地定量计算中,首次将逼近重心的两髂前下棘中点确定为骨盆原点建立了坐标系,极大的简化了运算过程和坐标系统,使之可应用于无相关基础的基层医院医生中,具有较大的推广价值和广阔的市场前景。
4、我科已初步将此图像后处理方法应用于骨盆骨折的手术复位中,相比于传统的手术治疗,具有手术切口小、复位速度快且正确率高、手术时间短、术中受辐射时间少、术后并发症较少及术后骨盆功能恢复良好等优点,Matta评分标准优良率达94.54%,已取得了好的微创固定效果。
附图说明
图1本发明基于坐标体系的骨盆CT三维重建的图像后处理方法流程图;
图2为均值骨盆三维模型确定的骨盆原点和骨盆三维坐标系示意图;
图3为骨盆骨折距离及角度移位定量计算示意图;
图4为骨盆骨折位移定性分析示意图;
图2中,a为均值骨盆三维模型D;b为医学建模软件中显示的均值骨盆的重心D的位置O;c-f为在均值骨盆上建立的骨盆三维立体坐标系D;c为骨盆正位,d为入口位,e为出口位,f为左侧位;
图3中,a-d为患者1,为X轴上的距离移位,a为骨盆正位,b为左侧位,c和d示距离移位的定量计算方法。c为使用方法一示意图,建立两个相互平行且垂直于X轴的平面M1和N1,通过平移两个平面,使得两平面均与移位骨块两侧相切,且保证相切的两点为相互之间为复位点。d为使用方法二示意图,具体操作方法见实施方案104)。
e-g为患者2,表示的为旋转移位;e示骨盆正位,f为骨盆出口位,g和h示旋转移位的定量计算方法。g为方法一,计算移位骨块相切的平面M2与移位轴所在的平面Y形成的夹角为α,而健侧骨块(或复位骨块)相切的平面N2与面Y的夹角为β,则此骨盆骨折沿着Z轴旋转移位的角度为|α-β|。h为方法二,通过求取均分平面La和Lb的方程来计算旋转角度,利用求取旋转角度θ。
图4中,a为距离移位;b为X轴旋转移位,用右手法则,该示意图为X旋负移位;c为Y轴旋转移位,左侧损伤用左手法则,c为Y旋正;d为Y轴旋转移位,右侧损伤用右手法则,为Y旋负;e为Z轴旋转移位,左侧损伤用左手法则,E为Z旋正;f为Z轴旋转移位,右侧损伤用右手法则,F为Z旋负。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。
实施例一:
201)采集CT数据:
在同一台CT机上筛选N个健康成人骨盆CT的原始数据,扫描参数为:电压120V,层厚1mm,层距1mm;扫描平面:使十字瞄准器对准耻骨联合上缘和双髂前上棘连线中点,下方到股骨中上段,上方位于中上腹平面;以DICOM格式存储并刻录于光盘CD-ROM上。
202)建立三维模型:
打开三维重建软件,选中目标的所有图像文件,将原始二维CT数据按默认方式导入,按照解剖位置设定图像方位得到三维模型Ai。在菜单栏上根据设定的骨组织阈值(>226HU)执行阈值分割,并予区域增长将图像做去除噪点处理。执行Edit Masks,利用工具Erase和布尔运算分离各关节,去除股骨和腰椎,得到骨盆(左右髋骨及骶尾骨)蒙版;然后手动分别在矢状面、水平面和冠状面逐层封闭蒙版缺损,填充空隙,并通过膨胀和腐蚀操作关闭小空洞和边缘。执行Calculate 3D操作,并设定重建质量为最佳,对上述骨盆蒙板进行三维重建,得到对应的完全封闭的正常骨盆三维模型Bi。
203)装配三维模型:
利用center of cavity操作对操作202)中获得的骨盆三维模型Bi求解其模型重心Bi,并建立以骨盆重心Bi为原点的右手坐标系:即面向骨盆正位水平面上过重心Bi指向右侧为X轴,经过重心Bi垂直于水平面向上为Z轴,经过重心Bi垂直于X、Z轴指向骶骨为Y轴。对纳入的每一个模型均执行上述操作,并进行模型的刚性变换即使得所有模型的重心坐标轴系和建模软件的世界坐标系匹配,重新定义骨盆的重心B0为骨盆坐标系原点,建立右手坐标系B0。
204)拟合对齐:
利用最小二乘法原理,固定一个正常骨盆三维模型Bj为参考模型的基础上,将所有其他模型Bi和参照模型Bj的骶尾骨部位进行最佳拟合对齐,目的是为了使得所有模型骶尾骨部位重叠部分最多,后期进行最小二乘法时相对应的采样点之间的距离最短,获得一个拟合对齐后的三维模型C,具体的数学运算如下:
设骨盆三维模型Bi上一点P(x1,x2,x3),对应参照模型Bj上的P'(y1,y2,y3),线性转换公式为则P到P'的转换矩阵为本发明中模型Bi和Bj上拟合采样点为Q个(Q→∞),并需使得Q个对应采样点之间的距离之和最小,即最小,其中d为对应点之间的距离,
求解上述方程,可解的未知数aij的值,即可得到模型Bi和Bj的转换矩阵,完成三维模型的最佳拟合对齐。
205)重建均值骨盆模型:
根据步骤的204)中建立的三维模型C。然后以基准坐标系B0的XY平面(即通过原点B0的水平面)为起始平面沿着Z轴进行等层距(1mm)平行双向切割,使得其与CT二维断层图像相互一致,获得截面曲线即为模型的轮廓线。然后将获得的所有轮廓线统一装配,在每一个截面上均构建平面网格,记录Z轴位点坐标,并获取和记录网格与曲线交点的点云数据。通过数学运算中的均值计算来获取均值点云数据,形成一系列平行于XY面的均值离散坐标点,利用正交多项拟合曲线法排除偏离主点的点云,并将离散坐标点按照多项拟合曲线走向进行有次序地连接和封闭。最后将所有平行于XY平面的曲线通过点云逆向三维重建形成均值三维模型D。
206)建立均值骨盆三维坐标系D:
通过步骤205),建立了均值骨盆三维模型D,将其导入医学建模软件中,再次使用Center of cavity求取其重心D,确定重心D的坐标位置位于骨盆上口与骨盆下口内,骶正中嵴上方平第四骶前孔(图2b),即面向骨盆正位水平面上过重心D指向右侧为X轴,经过重心D垂直于水平面向上为Z轴,经过重心D垂直于X、Z轴指向骶骨为Y轴。
所述N≥10,i为1-N的任意自然数,j为1-N的任意自然数。
101)采集病人骨盆CT数据:
采集病人术前骨盆二维CT的病例数据。
102)建立病例三维模型:
利用病例数据,建立病例三维模型,去除病例三维模型的股骨和腰椎,得到病例的骨盆三维模型,具体的,建立病例三维模型后,进行阈值分割,并予区域增长将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除病例三维模型的股骨和腰椎。
103)在病例的骨盆三维模型中选择均值骨盆重心D为原点建立笛卡尔三维坐标系D:
以206)步骤中的重心D为原点,建立坐标系D。
104)获得骨盆骨折的距离移位和旋转移位:
根据103)步骤中的坐标系D进一步的,计算骨折的距离移位及角度移位:方法一、通过模拟游标卡尺的方法即“两平面法”:即在建模软件中,建立两个相互平行且垂直于移位轴的平面,通过小心平移两个平面,使得两平面均与移位骨块断端两侧恰好相切,且相切的两点为相互复位点,此时通过参数计算的两平面之间的距离即为移位的距离。计算骨折的角度移位:即移位骨块和健侧骨块相切的平面与移位轴平面形成的夹角之差。
方法二、(1)利用建模软件以重心D为中心,以移位轴为对称轴,在对侧建立健侧骨盆的虚拟镜像U;
(2)通过阈值分割、编辑蒙版,分离出移位骨块A,扫描移位骨块上的几何特征,并将移位骨块W的空间信息以点集的形式记录Ai={(xi,yi,zi)}。将移位骨块A通过平移或者旋转使得与健侧镜像U恰好完全重合,得到B,采集此时的几何特征,以点集Bi={(xi,yi,zi)}记录。
(3)建立数学模型,利用采集的点集信息,分别求取Ai和Bi的均分平面La和Lb,则骨块移位程度的定量分析,可以通过La和Lb之间的距离及旋转移位来计算。具体计算方法为:设La的方程为a1x+b1y+c1z+d1=0,Lb为a2x+b2y+c2z+d2=0;根据La和Lb分别均分Ai和Bi,可知点集Ai和Bi上所有点分别到两平面的距离的平方之和最短,可利用点到面的距离公式及偏导数进行求解,即:可知A和B上的任一点到La和Lb上的距离为:
方法三:(1)同方法二得到移位骨块A,及可完全与镜像U完全重合的骨块B;(2)利用Geomagic Studio软件执行特征选择-均分平面操作,则可较粗略地得到均分移位骨块的平面La和Lb的方程矩阵,继而求解。
实施例二:由于实施例一中的坐标原点D的获得方法较为复杂,在临床上可采用靠近D点的髂前下棘连线中点E代替,尤其适用于在骨盆移位的定性分析或较粗略地定量计算中,具体数据见表1。
101)采集病人骨盆CT数据:
采集病人术前骨盆二维CT的病例数据。
102)建立病例三维模型:
利用病例数据,建立病例三维模型,去除病例三维模型的股骨和腰椎,得到病例的骨盆三维模型,具体的,建立病例三维模型后,进行阈值分割,并予区域增长将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除病例三维模型的股骨和腰椎。
103)在病例的骨盆三维模型中选择髂前下棘连线中点E为原点建立笛卡尔三维坐标系E:以髂前下棘连线中点为原点E,建立坐标系E,即面向骨盆正位水平面上过原点E指向右侧为X轴,经过原点E垂直于水平面向上为Z轴,经过原点E垂直于X、Z轴指向骶骨为Y轴。
104)获得骨盆骨折的距离移位和旋转移位:
方法一、通过模拟游标卡尺的方法即“两平面法”:即在建模软件中,建立两个相互平行且垂直于移位轴的平面,通过小心平移两个平面,使得两平面均与移位骨块断端两侧恰好相切,且相切的两点为相互复位点,此时通过参数计算的两平面之间的距离即为移位的距离。计算骨折的角度移位:即移位骨块和健侧骨块相切的平面与移位轴平面形成的夹角之差。
方法二、(1)利用建模软件以重心D为中心,以移位轴为对称轴,在对侧建立健侧骨盆的虚拟镜像U;
(2)通过阈值分割、编辑蒙版,分离出移位骨块A,扫描移位骨块上的几何特征,并将移位骨块W的空间信息以点集的形式记录Ai={(xi,yi,zi)}。将移位骨块A通过平移或者旋转使得与健侧镜像U恰好完全重合,得到B,采集此时的几何特征,以点集Bi={(xi,yi,zi)}记录。
(3)建立数学模型,利用采集的点集信息,分别求取Ai和Bi的均分平面La和Lb,则骨块移位程度的定量分析,可以通过La和Lb之间的距离及旋转移位来计算。具体计算方法为:设La的方程为a1x+b1y+c1z+d1=0,Lb为a2x+b2y+c2z+d2=0;根据La和Lb分别均分Ai和Bi,可知点集Ai和Bi上所有点分别到两平面的距离的平方之和最短,可利用点到面的距离公式及偏导数进行求解,即:可知A和B上的任一点到La和Lb上的距离为:
方法三:(1)同方法二得到移位骨块A,及可完全与镜像U完全重合的骨块B;(2)利用Geomagic Studio软件执行特征选择-均分平面操作,则可较粗略地得到均分移位骨块的平面La和Lb的方程矩阵,继而求解。
表1各标志点左右侧及连线中点的坐标及其到骨盆重心D的距离(单位:mm)
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (8)
1.基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于,包括以下步骤:
101)采集病人骨盆CT数据:采集病人术前骨盆二维CT的病例数据;
102)建立病例三维模型:利用病例数据,建立病例三维模型,去除病例三维模型的股骨和腰椎,得到病例的骨盆三维模型;
103)在病例的骨盆三维模型中选择均值骨盆重心D’为原点建立笛卡尔三维坐标系D”或选择髂前下棘连线中点E’为原点建立笛卡尔三维坐标系E”;
104)获得骨盆骨折的距离移位和旋转移位:根据103)步骤中的坐标系,利用数学方法分别精确计算骨折部投射在X、Y、Z轴的距离移位和旋转移位;
所述103)步骤中,以均值骨盆的重心D’为原点,建立坐标系D”,所述均值骨盆的重心D’以及坐标系D”的建立方法如下:
201)采集CT数据:采集N个健康成人骨盆二维CT的N组原始数据;
202)建立三维模型:利用N组原始数据,分别建立N个三维模型Ai,去除三维模型Ai的股骨和腰椎,得到N个骨盆的三维模型Bi;
203)装配三维模型:获取三维模型Bi的重心Bi’,以重心Bi’为原点,建立坐标系Bi”;将N个三维模型Bi连同坐标系Bi”一起重叠,得到基准坐标系B0’和基准原点B0;
204)拟合对齐:利用最小二乘法原理,固定任意一个三维模型Bj为参考模型,将所有其他三维模型Bi和三维模型Bj的骶尾骨部位进行最佳拟合对齐,获得一个拟合对齐后的三维模型C;
205)重建均值骨盆模型:以基准坐标系B0’的XY平面为起始平面沿着Z轴将三维模型C进行等层距平行双向切割,获得截面轮廓线,将获得的所有轮廓线统一装配,在每一个截面上均构建平面网格,记录Z轴位点坐标,并获取和记录网格与轮廓线交点的点云数据,通过数学运算中的均值计算来获取均值点云数据,形成一系列平行于XY面的均值离散坐标点(Xi,Yi),利用正交多项拟合曲线法排除偏离主点的点云,并将离散坐标点按照多项拟合曲线走向进行有次序地连接和封闭,将所有平行于XY平面的曲线通过点云逆向三维重建形成均值三维模型D;
206)建立均值骨盆三维坐标系D”:获取三维模型D的重心D’,以重心D’为原点,建立坐标系D”;
所述N≥10,i为1-N的任意自然数,j为1-N的任意自然数。
2.根据权利要求1所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述102)步骤中,建立病例三维模型后,进行阈值分割,将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除病例三维模型的股骨和腰椎。
3.根据权利要求1所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述202)步骤中,建立三维模型Ai后,进行阈值分割,将原始数据图像做去噪点处理,所述骨组织阈值>226HU,再进行编辑蒙版,最后再去除三维模型Ai的股骨和腰椎。
4.根据权利要求1所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述203)步骤中,所述坐标系Bi”面向骨盆正位水平面上过重心Bi’指向右侧为X轴,经过重心Bi’垂直于骨盆正位水平面向上为Z轴,经过重心Bi’垂直于X、Z轴指向骶骨为Y轴,在进行203)步骤之前,先对三维模型Bi进行表面松弛、重划网格优化处理。
6.根据权利要求1所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述205)步骤中,所述层距为1mm,所述206)步骤中,所述坐标系D”面向骨盆正位水平面上过重心D’指向右侧为X轴,经过重心D’垂直于水平面向上为Z轴,经过重心D’垂直于X、Z轴指向骶骨为Y轴。
7.根据权利要求1所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述103)步骤中,以髂前下棘连线中点E’为原点,建立坐标系E”,所述坐标系E”面向骨盆正位水平面上过髂前下棘连线中点E’指向右侧为X轴,经过髂前下棘连线中点E’垂直于骨盆正位水平面向上为Z轴,经过髂前下棘连线中点E’垂直于X、Z轴指向骶骨为Y轴。
8.根据权利要求7所述的基于坐标体系的骨盆CT三维重建的图像后处理方法,其特征在于:所述104)步骤中,计算骨折的距离移位及角度移位:
方法一、通过模拟游标卡尺的方法即“两平面法”:即在建模软件中,建立两个相互平行且垂直于移位轴的平面,通过小心平移两个平面,使得两平面均与移位骨块断端两侧恰好相切,且相切的两点为相互复位点,此时通过参数计算的两平面之间的距离即为移位的距离,计算骨折的角度移位:即移位骨块和健侧骨块相切的平面与移位轴平面形成的夹角之差;
方法二、包括以下步骤:
401)利用建模软件以重心D’为中心,以移位轴为对称轴,在对侧建立健侧骨盆的虚拟镜像U;
402)通过阈值分割、编辑蒙版,分离出移位骨块A,扫描移位骨块上的几何特征,并将移位骨块W的空间信息以点集的形式记录Ai={(xi,yi,zi)},将移位骨块A通过平移或者旋转使得与健侧镜像U恰好完全重合,得到B,采集此时的几何特征,以点集Bi={(xi,yi,zi)}记录;
403)建立数学模型,利用采集的点集信息,分别求取Ai和Bi的均分平面La和Lb,则骨块移位程度的定量分析,可以通过La和Lb之间的距离及旋转移位来计算,具体计算方法为:设La的方程为a1x+b1y+c1z+d1=0,Lb为a2x+b2y+c2z+d2=0;根据La和Lb分别均分Ai和Bi,可知点集Ai和Bi上所有点分别到两平面的距离的平方之和最短,可利用点到面的距离公式及偏导数进行求解,即:可知A和B上的任一点到La和Lb上的距离为:
方法三:501)同方法二得到移位骨块A,及可完全与镜像U完全重合的骨块B;502)利用Geomagic Studio软件执行特征选择-均分平面操作,则可较粗略地得到均分移位骨块的平面La和Lb的方程矩阵,继而求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710114304.4A CN107016666B (zh) | 2017-02-28 | 2017-02-28 | 基于坐标体系的骨盆ct三维重建的图像后处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710114304.4A CN107016666B (zh) | 2017-02-28 | 2017-02-28 | 基于坐标体系的骨盆ct三维重建的图像后处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016666A CN107016666A (zh) | 2017-08-04 |
CN107016666B true CN107016666B (zh) | 2020-09-18 |
Family
ID=59440301
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710114304.4A Expired - Fee Related CN107016666B (zh) | 2017-02-28 | 2017-02-28 | 基于坐标体系的骨盆ct三维重建的图像后处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016666B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945157B (zh) * | 2017-11-15 | 2021-07-30 | 上海交通大学医学院附属第九人民医院 | 骶骨上终板测量方法/系统、计算机可读存储介质及设备 |
CN113272869A (zh) * | 2018-12-05 | 2021-08-17 | 西门子医疗有限公司 | 医学成像中从定位片进行三维形状重构 |
CN109820590B (zh) * | 2019-02-15 | 2024-04-12 | 中国人民解放军总医院 | 一种骨盆骨折复位智能监控系统 |
CN109984843B (zh) * | 2019-03-20 | 2021-08-06 | 中国科学院苏州生物医学工程技术研究所 | 骨折闭合复位导航系统及方法 |
CN111590584B (zh) * | 2020-05-27 | 2021-12-10 | 京东方科技集团股份有限公司 | 安全限位区的确定方法和装置、复位方法和医用机器人 |
CN112270744A (zh) * | 2020-11-03 | 2021-01-26 | 广州三瑞医疗器械有限公司 | 一种基于骨盆模型的产轴线构建方法 |
CN113469935B (zh) * | 2021-04-15 | 2023-06-13 | 卡本(深圳)医疗器械有限公司 | 一种基于ct图像的髂后上棘自动检测与定位方法 |
CN113298784B (zh) * | 2021-05-25 | 2023-03-21 | 中南大学 | 基于视角的骨折线提取映射方法 |
CN115249290B (zh) * | 2022-06-10 | 2023-05-16 | 首都医科大学附属北京友谊医院 | 针对单侧颞骨的空间数据处理方法、空间定位方法及设备 |
CN115908951B (zh) * | 2023-01-06 | 2023-05-26 | 长江勘测规划设计研究有限责任公司 | 一种机器视觉实时识别预紧螺栓状态的方法 |
CN117132747B (zh) * | 2023-10-25 | 2024-03-19 | 北京爱康宜诚医疗器材有限公司 | 基于骨骼模型的骨骼复位方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102314712A (zh) * | 2011-09-23 | 2012-01-11 | 暨南大学 | 基于三维成像技术的骨盆建模与骨盆参数测算方法 |
CN105139442A (zh) * | 2015-07-23 | 2015-12-09 | 昆明医科大学第一附属医院 | 一种结合ct和mri二维图像建立人体膝关节三维仿真模型的方法 |
CN106264731A (zh) * | 2016-10-11 | 2017-01-04 | 昆明医科大学第附属医院 | 一种基于点对点配准技术虚拟膝关节单髁置换术模型构建的方法 |
CN106447764A (zh) * | 2016-09-08 | 2017-02-22 | 福州大学 | 一种三维可视化的人体骨盆参数自动测量方法 |
CN106420053A (zh) * | 2016-09-05 | 2017-02-22 | 上海逸动医学科技有限公司 | 骶髂关节置入螺钉固定s2椎体的通道定位方法及系统 |
-
2017
- 2017-02-28 CN CN201710114304.4A patent/CN107016666B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102314712A (zh) * | 2011-09-23 | 2012-01-11 | 暨南大学 | 基于三维成像技术的骨盆建模与骨盆参数测算方法 |
CN105139442A (zh) * | 2015-07-23 | 2015-12-09 | 昆明医科大学第一附属医院 | 一种结合ct和mri二维图像建立人体膝关节三维仿真模型的方法 |
CN106420053A (zh) * | 2016-09-05 | 2017-02-22 | 上海逸动医学科技有限公司 | 骶髂关节置入螺钉固定s2椎体的通道定位方法及系统 |
CN106447764A (zh) * | 2016-09-08 | 2017-02-22 | 福州大学 | 一种三维可视化的人体骨盆参数自动测量方法 |
CN106264731A (zh) * | 2016-10-11 | 2017-01-04 | 昆明医科大学第附属医院 | 一种基于点对点配准技术虚拟膝关节单髁置换术模型构建的方法 |
Non-Patent Citations (2)
Title |
---|
"不稳定骨盆复位架的研发与应用";赵晶鑫;《中国博士学位论文全文数据库 医药卫生科技辑》;20151015(第10期);正文第3.1-3.2节 * |
"基于CT及逆向工程软件构建正常人体骨盆三维有限元模型";刘小龙等;《中国组织工程研究》;20140924;第18卷(第40期);正文第1节 * |
Also Published As
Publication number | Publication date |
---|---|
CN107016666A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016666B (zh) | 基于坐标体系的骨盆ct三维重建的图像后处理方法 | |
US11826111B2 (en) | Surgical navigation of the hip using fluoroscopy and tracking sensors | |
CN106780715B (zh) | 一种骨盆ct三维坐标体系的建立方法 | |
US11553969B1 (en) | System for computation of object coordinates accounting for movement of a surgical site for spinal and other procedures | |
Chapuis et al. | A new system for computer-aided preoperative planning and intraoperative navigation during corrective jaw surgery | |
Taylor et al. | Computer-integrated revision total hip replacement surgery: concept and preliminary results | |
AU2009241518B2 (en) | System and method for image segmentation in generating computer models of a joint to undergo arthroplasty | |
CN105832415B (zh) | 用于股骨精准去旋转的导航装置及其制作方法和使用方法 | |
US20010021806A1 (en) | System and method for fusing three-dimensional shape data on distorted images without correcting for distortion | |
TW201801682A (zh) | 影像增強真實度之方法與應用該方法在可穿戴式眼鏡之手術導引 | |
EP4086854A1 (en) | Surface and image integration for model evaluation and landmark determination | |
CN105361883A (zh) | 膝关节置换三维空间下肢生物力线的确定方法 | |
CN111166474B (zh) | 一种关节置换手术术前的辅助诊查方法和装置 | |
Pietruski et al. | Accuracy of experimental mandibular osteotomy using the image-guided sagittal saw | |
JP6806655B2 (ja) | 放射線撮像装置、画像データ処理装置及び画像処理プログラム | |
CN105894508A (zh) | 一种医学图像自动定位质量的评估方法 | |
JP4512824B2 (ja) | 画像処理装置及びそれに用いられるプログラム | |
CN115153835A (zh) | 基于特征点配准与增强现实的髋臼假体放置引导系统及方法 | |
Irwansyah et al. | Algorithm for segmentation and reduction of fractured bones in computer-aided preoperative surgery | |
Liu et al. | Virtual plate pre-bending for the long bone fracture based on axis pre-alignment | |
WO2023198118A1 (en) | Method and apparatus for determining acetabulum-to-femoral-head distance | |
CN112472292B (zh) | 一种用于骨科手术分析模型的构建方法及系统 | |
Roner et al. | Three-dimensional automated assessment of the distal radioulnar joint morphology according to sigmoid notch surface orientation | |
US20240024035A1 (en) | Preoperative imaging combined with intraoperative navigation before and after alteration of a surgical site to create a composite surgical three dimensional structural dataset | |
CN114403911B (zh) | 一种基于ct三维重建图像的尺骨近端参数测量方法及装置 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200918 |