CN108510580A - 一种椎骨ct图像三维可视化方法 - Google Patents
一种椎骨ct图像三维可视化方法 Download PDFInfo
- Publication number
- CN108510580A CN108510580A CN201810266578.XA CN201810266578A CN108510580A CN 108510580 A CN108510580 A CN 108510580A CN 201810266578 A CN201810266578 A CN 201810266578A CN 108510580 A CN108510580 A CN 108510580A
- Authority
- CN
- China
- Prior art keywords
- vertebra
- image
- voxel
- tri patch
- vertex
- 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.)
- Pending
Links
Classifications
-
- 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
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种椎骨CT图像三维可视化方法,属于医学领域,为了弥补传统Marching Cubes算法在进行三维重建可视化的过程中,处理体素和等值点数量过多的问题,本发明提供一种椎骨CT图像三维可视化方法,在传统Marching Cubes算法理论基础之上进行改进,确定体素内三角面片剖分构型时,利用黄金分割方法计算等值面与体素棱边交点,将三角面片拟合成等值面,再计算三角面片的法向量,采用中心差分方法求解体素内8个顶点的梯度值,并通过可视化工具包VTK(Visualization Toolkit)完成显示、旋转、缩放等交互操作,实现透明、半透明可视化功能设置,以便于医生和研究人员能够从任意角度观察和分析,为观察椎骨内、外结构信息和后期的治疗、研究提供理论依据和技术支持。
Description
技术领域
本发明属于医学领域,特别涉及了椎骨临床外科和科学研究中一种椎骨CT图像三维可视化方法。
背景技术
Marching Cubes算法是三维重建可视化方法面绘制中经典的算法,原理简单、重建效率高,但其重建结果的质量和速度还有待进一步提高,国内外的研究人员相继提出了区间树硬件加速,提高算法处理效率,Lopes提出了提高绘制精度的方法,针对传统Marching Cubes算法在进行三维重建可视化的过程中,处理体素和等值点数量过多的问题还有待进一步改进,三维重建可视化的效果直接影响临床手术的成功率和科学研究的精确度,尤其在椎弓根螺钉精准置入手术中,对椎骨三维可视化过程的时间和精度都有更高的要求,所以,设计出一种更适合的椎骨CT图像三维可视化方法至关重要。
发明内容
本发明目标在于针对传统Marching Cubes算法在进行三维重建可视化的过程中,处理体素和等值点数量过多的问题,提供一种高效的椎骨CT图像三维可视化方法以解决上述背景技术中的问题。
为实现上述目标,本发明提供如下技术方案:
一种椎骨CT图像三维可视化方法,进行图像预处理后,在传统Marching Cubes算法基础之上进行改进,在确定体素内三角面片剖分构型时,利用黄金分割方法计算等值面与体素棱边的交点,即三角面片的顶点,将获得的三角面片近似拟合成等值面,为了使重建后的图像更加逼真加入光照模型,再计算三角面片的法向量,直接计算比较复杂,采用中心差分方法求解出体素内8个顶点的梯度值,并通过可视化工具包VTK (Visualization Toolkit)完成显示、旋转、缩放等交互操作,实现透明、半透明可视化功能设置,以便于医生和研究人员能够从任意角度观察和分析,为观察椎骨内、外结构信息和后期的治疗、研究提供理论依据和技术支持。
作为本发明优选:确定三角剖分构型时,利用黄金分割方法代替传统线性插值法,公共棱边上的顶点计算次数由4减少到1。
作为本发明优选:利用VTK进行个性化的可视化操作,并且可以根据操作者需要进行个性化功能设置。
与现有技术相比,本发明的有益效果为:所述椎骨CT图像三维可视化方法,与传统Marching Cubes算法相比,在计算等值面与体素棱边交点过于冗杂的问题上进行相应改进,对于相邻体素,公共棱边只进行一次计算,这样在不影响绘制效果的同时,大大降低了处理时间,提高了算法的执行效率。
附图说明
图1为椎骨CT图像三维可视化方法流程图;
图2为四个相邻的体素三维结构示意图;
图3为重建结果示意图。
具体实施方式:
下面将结合本发明实施例中的附图,对本发明中的具体实施方案进行详细的说明,很明显,所说明的实施例只是本发明中的一部分,并不是全部实施例,对于本发明中的实施例,在本领域的普通技术人员在没有付出创造性劳动而获得的,都属于本发明所保护的范围。
请参阅图1~3,本发明实施例中:一种椎骨CT图像三维可视化方法,进行图像预处理后,在传统Marching Cubes算法基础之上进行改进,在确定体素内三角面片剖分构型时,利用黄金分割方法计算等值面与体素棱边的交点,即三角面片的顶点,将获得的三角面片近似拟合成等值面,为了使重建后的图像更加逼真加入光照模型,再计算三角面片的法向量,直接计算比较复杂,采用中心差分方法求解出体素内8个顶点的梯度值,并通过VTK完成显示、旋转、缩放等交互操作,实现透明、半透明可视化功能设置,以便于医生和研究人员能够从任意角度观察和分析,为观察椎骨内、外结构信息和后期的治疗、研究提供理论依据和技术支持。
具体的工作流程为:进行图像预处理后,在传统Marching Cubes算法基础之上进行改进,在确定体素内三角面片剖分构型时,利用黄金分割方法计算等值面与体素棱边的交点,即三角面片的顶点,将获得的三角面片近似拟合成等值面,为了使重建后的图像更加逼真加入光照模型,再计算三角面片的法向量,直接计算比较复杂,采用中心差分方法求解出体素内8个顶点的梯度值,并通过VTK完成显示、旋转、缩放等交互操作,实现透明、半透明可视化功能设置,以便于医生和研究人员能够从任意角度观察和分析,为观察椎骨内、外结构信息和后期的治疗、研究提供理论依据和技术支持。
1、确定三角面片顶点
传统的Marching Cubes算法中8个相邻的立方体中一个公共的顶点需要计算8次,四个相邻立方体中的一条公共棱边,需要计算4次,严重影响运算的速度,在确定等值面与棱边交点位置时利用棱边的黄金分割点来代替线性插值点,再分别计算交点的坐标,不需要进行重复计算,公共棱边计算次数由4减少到1,计算过程简化,时间降低,算法执行效率整体提高。
2、VTK可视化设置
VTK提供图像处理和可视化的能力,在VTK中,二维图像数据和三维体数据表示为vtkImagineData对象,VTK将常用的重建算法封装到类中用于三维重建,VTK类库作为一个开源可视化软件,可以根据不同研究需要进行相应的改进,对试验的成功至关重要。
利用VTK将三维重建可视化后的每一节椎骨根据不同需要,可以设定不同的不透明度,也可以将病灶区域所在的椎骨设置成透明,有利于观察椎骨内部结构,对后期的椎弓根螺钉置入轨迹分析至关重要。
通过改进的Marching Cubes移动立方体算法,进行椎骨CT图像三维可视化,经过测试获得实验结果,且效率更高。
上述描述清楚的说明了本发明的技术方案、流程及优势,本领域的技术人员显然理解,本发明不因上述实施例而受到限制,上述描述的实施例与说明书只是本发明的技术方案及原理并不代表全部,在不背弃本发明精神和内容的前提下,本发明进行相应算法的改进,都在本发明要求保护的范围之内,以特有的形式实现本发明的实验结果。本发明所保护的范围由所附的权利要求书及等同要件限定。
Claims (3)
1.一种椎骨CT图像三维可视化方法,进行图像预处理后,在传统Marching Cubes算法基础之上进行改进,在确定体素内三角面片剖分构型时,利用黄金分割方法计算等值面与体素棱边的交点,即三角面片的顶点,将获得的三角面片近似拟合成等值面,为了使重建后的图像更加逼真加入光照模型,再计算三角面片的法向量,直接计算比较复杂,采用中心差分方法求解出体素内8个顶点的梯度值,并通过可视化工具包VTK (VisualizationToolkit)完成显示、旋转、缩放等交互操作,实现透明、半透明可视化功能设置,以便于医生和研究人员能够从任意角度观察和分析,为观察椎骨内、外结构信息和后期的治疗、研究提供理论依据和技术支持。
2.根据权利要求1所述一种椎骨CT图像三维可视化方法,其特征在于,确定三角剖分构型时,利用黄金分割方法代替传统线性插值法,公共棱边上的顶点计算次数由4减少到1。
3.根据权利要求1所述一种椎骨CT图像三维可视化方法,其特征在于,利用VTK进行个性化的可视化操作,并且可以根据操作者需要进行个性化功能设置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810266578.XA CN108510580A (zh) | 2018-03-28 | 2018-03-28 | 一种椎骨ct图像三维可视化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810266578.XA CN108510580A (zh) | 2018-03-28 | 2018-03-28 | 一种椎骨ct图像三维可视化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108510580A true CN108510580A (zh) | 2018-09-07 |
Family
ID=63378960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810266578.XA Pending CN108510580A (zh) | 2018-03-28 | 2018-03-28 | 一种椎骨ct图像三维可视化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108510580A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109727306A (zh) * | 2018-11-20 | 2019-05-07 | 哈尔滨理工大学 | 一种基于vtk的脊柱医学图像三维可视化方法 |
CN110264504A (zh) * | 2019-06-28 | 2019-09-20 | 北京国润健康医学投资有限公司 | 一种用于增强现实的三维配准方法和系统 |
CN110796693A (zh) * | 2019-09-11 | 2020-02-14 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN111476764A (zh) * | 2020-03-27 | 2020-07-31 | 西安科技大学 | 一种运动模糊ct图像三维重建的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625766A (zh) * | 2009-08-03 | 2010-01-13 | 深圳先进技术研究院 | 医学图像处理方法 |
US20110254845A1 (en) * | 2010-04-16 | 2011-10-20 | Hitachi Medical Corporation | Image processing method and image processing apparatus |
CN103824328A (zh) * | 2014-03-27 | 2014-05-28 | 东北大学 | 一种基于OpenCL的移动立方体方法 |
CN106821496A (zh) * | 2016-12-28 | 2017-06-13 | 妙智科技(深圳)有限公司 | 一种经皮椎间孔镜手术精准规划系统及方法 |
-
2018
- 2018-03-28 CN CN201810266578.XA patent/CN108510580A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625766A (zh) * | 2009-08-03 | 2010-01-13 | 深圳先进技术研究院 | 医学图像处理方法 |
US20110254845A1 (en) * | 2010-04-16 | 2011-10-20 | Hitachi Medical Corporation | Image processing method and image processing apparatus |
CN103824328A (zh) * | 2014-03-27 | 2014-05-28 | 东北大学 | 一种基于OpenCL的移动立方体方法 |
CN106821496A (zh) * | 2016-12-28 | 2017-06-13 | 妙智科技(深圳)有限公司 | 一种经皮椎间孔镜手术精准规划系统及方法 |
Non-Patent Citations (3)
Title |
---|
帅仁俊 等: ""一种改进的MC三维重建算法"", 《中国数字医学》 * |
杨琴: ""基于VTK的虚拟手可视化模型的初步研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
马凯: ""腰骶段三维重建与测量技术的研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109727306A (zh) * | 2018-11-20 | 2019-05-07 | 哈尔滨理工大学 | 一种基于vtk的脊柱医学图像三维可视化方法 |
CN110264504A (zh) * | 2019-06-28 | 2019-09-20 | 北京国润健康医学投资有限公司 | 一种用于增强现实的三维配准方法和系统 |
CN110264504B (zh) * | 2019-06-28 | 2021-03-30 | 北京国润健康医学投资有限公司 | 一种用于增强现实的三维配准方法和系统 |
CN110796693A (zh) * | 2019-09-11 | 2020-02-14 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN110796693B (zh) * | 2019-09-11 | 2023-03-21 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN111476764A (zh) * | 2020-03-27 | 2020-07-31 | 西安科技大学 | 一种运动模糊ct图像三维重建的方法 |
CN111476764B (zh) * | 2020-03-27 | 2021-08-24 | 西安科技大学 | 一种运动模糊ct图像三维重建的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104794758B (zh) | 一种三维图像的裁剪方法 | |
CN108510580A (zh) | 一种椎骨ct图像三维可视化方法 | |
CN101271574B (zh) | 三维可视化的方法、装置 | |
US8929635B2 (en) | Method and system for tooth segmentation in dental images | |
CN103106685B (zh) | 一种基于gpu的腹部脏器三维可视化方法 | |
CN107808156A (zh) | 感兴趣区域提取方法 | |
US8970581B2 (en) | System and method for interactive contouring for 3D medical images | |
WO2017020580A1 (zh) | 一种医学组织的切片图像三维重建的方法及设备 | |
US20110254845A1 (en) | Image processing method and image processing apparatus | |
CN108537750A (zh) | 一种快速的脊柱ct图像三维重建方法 | |
CN102074039B (zh) | 一种体绘制裁剪面绘制方法 | |
CN106960439B (zh) | 一种脊椎骨识别装置及方法 | |
CN106548476B (zh) | 利用医学图像统计肺部三维特征形状方法 | |
CN109360219A (zh) | 一种增强现实辅助手术方法及系统 | |
Shi et al. | A new volumetric geometric model for cutting procedures in surgical simulation | |
JP2012085833A (ja) | 3次元医用画像データの画像処理システム、その画像処理方法及びプログラム | |
CN107170009B (zh) | 一种基于医学图像的护目镜基弯数据测量方法 | |
CN102074051B (zh) | 一种体绘制平行投影的数据包围体快速定位方法 | |
CN107967682A (zh) | 一种两点式网格模型裁切方法 | |
Cheng et al. | Research on medical image three dimensional visualization system | |
Schein et al. | Adaptive extraction and visualization of silhouette curves from volumetric datasets | |
Liu et al. | Optimization of reconstruction of 2D medical images based on computer 3D reconstruction technology | |
CN113679417A (zh) | 模型引导优化并行的超声影像3d重建方法 | |
Segars et al. | The MCAT, NCAT, XCAT, and MOBY computational human and mouse phantoms | |
Lv et al. | Interactive curved planar reformation based on snake model |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180907 |
|
WD01 | Invention patent application deemed withdrawn after publication |