一种基于GPU的腹部脏器三维可视化方法技术

一种基于GPU的腹部脏器三维可视化的方法，属于医学图像处理技术领域，计算机X射线断层扫描设备对人体腹部进行CT扫描；获取腹部CT图像，并对获取的腹部CT图像进行三维可视化处理，得到三维可视化处理的结果图像；对腹部CT图像进行对应点匹配插值得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据；进行光线投射，得到腹部CT图像光线投射结果图像；将腹部CT图像的三维可视化处理的结果图像在装有独立GPU的计算机屏幕上显示。在对医学图像三维可视化处理时先进行插值，使断层图像的层间隔减小；基于GPU进行光线投射，利用GPU高并发的特性、有效地降低计算时间。

【技术实现步骤摘要】

本专利技术属于医学图像处理
，具体涉及一种基于GPU的腹部脏器三维可视化方法。
技术介绍
医学影像三维可视化技术是指基于科学计算可视化技术，从一系列二维断层医学图像中重建出三维图像，利用计算机显示并完成三维图像操作的技术。首先，对CT、MRI等医学体数据进行三维重建得到三维模型，对体数据的绘制结果从各个方向进行投影显示，便于医生从不同方向对人体内部结构进行观察，能辅助医生对感兴趣区域或感兴趣器官的大小形状和空间位置获得定量描述，帮助医生了解复杂的解剖细节，提供直观的视觉印象。在此基础上结合实际临床应用，实现计算机与用户的交互，对可视化结果进行任意平面切害I]，可以帮助实现手术模拟、手术规划和虚拟内窥镜等辅助诊断操作。提高医生诊断的准确率和医院的工作效率。人体内部结构十分复杂，通常医生无法观察到手术进行的实际情况，而且手术具有高危险性，不能预先在人体上进行模拟，利用可视化技术，在重构出的三维图像基础上进行计算机模拟，可以设计、选择最佳手术方案。还可以在手术过程中对手术进行情况在屏幕上予以监视，使医生们做到心中有数，从而提高手术的成功率。其次，目前国内多数医院仍采用传统的胶片形式帮助医生诊断，大量的胶片不但存储是很大问题，也是不小的开支。实现医院数字化，不但方便医生诊断、医院管理，更能节约开支。因而医学图像三维可视化技术对于促进图像处理和可视化技术的发展和应用将起到积极的推动作用，研究并实现二维医学图像的三维可视化，具有重要的价值。医学数据可视化算法分为两大类，基于构造中间几何图元的面绘制算法和体绘制算法。面绘制需要先从体数据中重建物体的表面模型，然后利用传统图形学技术及硬件实现表面的绘制与显示；而体绘制则直接以体素为基本单元，应用视觉原理，通过对体数据重新采样来合成产生三维图像。面绘制是从3D数据场中抽取有意义和直观信息的一种重要手段。它实际上是把体数据转换成一种逼近面表示，从而可以进一步利用计算机图形学技术，甚至已有的硬件加速技术完成感兴趣信息的提取。由于它借助于面表不这样一个中间转换过程，而不是直接把体数据投向屏幕进行绘制，有时我们又称之为间接体绘制。这种三维几何表达的模型最常用的就是表面模型，一般以平面片特别是三角面片来逼近表示。面绘制重建的方法可分为基于轮廓的表面重建和基于体素的表面重建，基于体素的表面重建最具代表性的是Cuberille 算法，Marching Cubes 算法和 Dividing Cubes 算法。体绘制是以某种方式将整个数据场半透明地投影到2D屏幕上，并不借助中间几何图形。直接由三维数据场产生屏幕上的二维图像。这种算法能产生三维数据场的整体图像，包括每一个细节，并具有图像质量高、便于并行处理等优点。其主要问题是，计算量很大，且难以利用传统的图形硬件实现绘制，因而计算时间较长。目前应用于医学体数据的体绘制算法主要有基于图像空间的光线投射(Raycasting)算法,基于物体空间扫描的派射(Splatting)算法,错切-变形(shear_warp)算法，硬件辅助的三维纹理映射(hardware-assisted 3D texture mapping)算法。体绘制方法都可以通过积分函数来进行计算，光线投射算法是一种以图像为序的直接体绘制方法，最早由M.Levoy于1988年提出，也是体绘制中最基本和最灵活的算法。其基本思想及特点是:首先对于图像中要绘制的像素，根据视点方向，发射出一条射线，穿过像素中心向物体投射一单束光线，沿着射线，从遇到体数据开始，计算并整合遇到的体数据的光学属性。光线投射算法采用为每个体素分配透明度和颜色值的方法来合成图像，因此有利于保留图像的细节，绘制高品质的图像，特别适用于绘制区域特征模糊，体素特征相关性高的图像。但因为需要遍历每个体素，而且当观察方向发生变化时，数据中的采样点之间的前后关系也发生变化，这样就需要进行重新采样和计算，因此计算量较大。随着医学影像技术的发展，成像分辨率日益提升。光线投射法的每一次绘制都需要对整个体数据进行计算和处理，处理数据量巨大，而传统的腹部脏器三维可视化系统主要由基于CPU的编程语言(如c\c++等)实现，运行时间较慢，缺乏可以用于医学诊断的实时性，使得该系统在实际应用中有一定的限制。因此系统中可视化算法性能上的提升，对于整个医学诊断工作和计算机图形图像处理工作都有着极为重要的价值。在腹部脏器可视化过程中，存在着高度的并行性，因此适合用GPU来进行加速。与CPU相比，GPU具有并行的计算架构以及超长流水线，这种架构虽然在处理少量数据时速度不及频率更高的CPU，但在处理如医学图像等大量数据时速度远远超过CPU。同时GPU使用SIMD (Single Instruction Multiple Data,单指令多数据流)技术,对于需要并行计算的数据只需编写单个运算操作即可实现并行计算过程。综上所述基于GPU的腹部脏器三维可视化系统相对于传统基于CPU的该系统，运算效率更高，能够保证系统的实时性，具有更高的应用价值。GLSL(The OpenGL Shading Language)是 OpenGL 着色语言的简称，它是一种用于创建OpenGL着色器的语言。作为OpenGL 2.0标准的一部分，GLSL允许应用程序显式的指定处理顶点和处理片段时进行的操作。借助GLSL能够实现基于GPU的腹部脏器三维可视化方法，提升系统的实时性。
技术实现思路
针对现有技术存在的问题，本专利技术提供一种基于GPU的腹部脏器三维可视化方法。本专利技术的技术方案是:一种基于GPU的腹部脏器三维可视化的方法，包括以下步骤:步骤1:计算机X射线断层扫描设备对人体腹部进行CT扫描；步骤2:获取腹部CT图像，每组腹部CT图像包含若干张腹部断层扫描图像；步骤3:对获取的腹部CT图像进行三维可视化处理，得到三维可视化处理的结果图像；步骤3.1:对腹部CT图像进行对应点匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据；步骤3.2:从人在计算机内的虚拟视点经由光线投射初始空白图像中每一像素点向体数据中发射虚拟光线，进行光线投射，计算每一像素点的光线投射显示结果，所有像素点光线投射显示结果的集合构成腹部CT图像光线投射结果图像；步骤3.3:对体数据进行任意平面切割，将非感兴趣区域从体数据中切除，对切割后的体数据进行光线投射，得到切割后体数据的光线投射结果图像；步骤3.4:若需要同时显示多个体数据，则进行多体素光线投射，将不同体数据的光线投射结果图像合成至同一结果图像内作为三维可视化处理的结果图像，否则直接将切割后体数据的光线投射结果图像作为三维可视化处理的结果图像；步骤4:将腹部CT图像的三维可视化处理的结果图像在装有独立GPU的计算机屏幕上显示。所述步骤3.1对腹部CT图像进行对应点匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据，具体步骤如下:步骤3.1.1:对腹部CT图像进行匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像；步骤3.1.2:对每个新插值图像上的所有点，计算其相邻前后两张腹部断层扫描图本文档来自技高网...

【技术保护点】
一种基于GPU的腹部脏器三维可视化的方法，其特征在于：包括以下步骤：步骤1：计算机X射线断层扫描设备对人体腹部进行CT扫描；步骤2：获取腹部CT图像，每组腹部CT图像包含若干张腹部断层扫描图像；步骤3：对获取的腹部CT图像进行三维可视化处理，得到三维可视化处理的结果图像；步骤3.1：对腹部CT图像进行对应点匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据；步骤3.2：从人在计算机内的虚拟视点经由光线投射初始空白图像中每一像素点向体数据中发射虚拟光线，进行光线投射，计算每一像素点的光线投射显示结果，所有像素点光线投射显示结果的集合构成腹部CT图像光线投射结果图像；步骤3.3：对体数据进行任意平面切割，将非感兴趣区域从体数据中切除，对切割后的体数据进行光线投射，得到切割后体数据的光线投射结果图像；步骤3.4：若需要同时显示多个体数据，则进行多体素光线投射，将不同体数据的光线投射结果图像合成至同一结果图像内作为三维可视化处理的结果图像，否则直接将切割后体数据的光线投射结果图像作为三维可视化处理的结果图像；步骤4：将腹部CT图像的三维可视化处理的结果图像在装有独立GPU的计算机屏幕上显示。...

【技术特征摘要】
1.一种基于GPU的腹部脏器三维可视化的方法，其特征在于:包括以下步骤: 步骤1:计算机X射线断层扫描设备对人体腹部进行CT扫描； 步骤2:获取腹部CT图像，每组腹部CT图像包含若干张腹部断层扫描图像； 步骤3:对获取的腹部CT图像进行三维可视化处理，得到三维可视化处理的结果图像； 步骤3.1:对腹部CT图像进行对应点匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据； 步骤3.2:从人在计算机内的虚拟视点经由光线投射初始空白图像中每一像素点向体数据中发射虚拟光线，进行光线投射，计算每一像素点的光线投射显示结果，所有像素点光线投射显示结果的集合构成腹部CT图像光线投射结果图像； 步骤3.3:对体数据进行任意平面切割，将非感兴趣区域从体数据中切除，对切割后的体数据进行光线投射，得到切割后体数据的光线投射结果图像； 步骤3.4:若需要同时显示多个体数据，则进行多体素光线投射，将不同体数据的光线投射结果图像合成至同一结果图像内作为三维可视化处理的结果图像，否则直接将切割后体数据的光线投射结果图像作为三维可视化处理的结果图像； 步骤4:将腹部CT图像的三维可视化处理的结果图像在装有独立GPU的计算机屏幕上显不O2.根据权利要求1所述的基于GPU的腹部脏器三维可视化的方法，其特征在于:所述步骤3.1对腹部CT图像进行对应点匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像，由新插值图像与原腹部CT图像得到匹配插值后的腹部CT图像即体数据，具体步骤如下: 步骤3.1.1:对腹部CT图像进行匹配插值，在每相邻前后两张腹部断层扫描图像间得到新插值图像； 步骤3.1.2:对每个新插值图像上的所有点，计算其相邻前后两张腹部断层扫描图像中对应点的灰度差， 并设定灰度差阈值； 步骤3.1.3:判断计算出的灰度差是否大于灰度差阈值:是，则设定搜索窗口半径和灰度阈值，在相邻前后两张腹部断层扫描图像中分别以该两张图像的对应点为中心以搜索窗口半径为半径的匹配窗，在匹配窗内选取最佳匹配点以决定当前点的灰度值；否则，通过三次样条插值求当前点的灰度值； 步骤3.1.4:得到每个新插值图像上的所有点的灰度值； 步骤3.1.5:由原腹部CT图像灰度值与每个新插值图像上的所有点的灰度值，得到匹配插值后的腹部CT图像即体数据，体数据中的所有像素点构成一个三维的数据场。3.根据权利要求1所述的基于GPU的腹部脏器三维可视化的方法，其特征在于:所述步骤3.2从人在计算机内的虚拟视点经由光线投射初始空白图像中每一像素点向体数据中发射虚拟光线，进行光线投射，计算每一像素点的光线投射显示结果，所有像素点光线投射显示结果的集合构成腹部CT图像光线投射结果图像，具体步骤如下: 步骤3.2.1:将体数据作为三维纹理加载至GPU的缓存中； 步骤3.2.2:设定颜色传递函数与透明度传递函数，并将这两个传递函数作为一维纹理分别加载至GPU的缓存中； 步骤3.2.3:从人在计算机内的虚拟视点经由光线投射初始空白图像中每一像素点向体数据中发射虚拟光线，进行光线投射； 步骤3.2.4:在GPU内并行处理所有投射的光线，沿光线进行重采样，对采样结果进行合成，并根据设定的颜色传递函数与透明度传递函数计算每一像素点的颜色和不透明度，即得到每一像素点的光线投射显示结果； 步骤3.2.5:所有光线投射显示结果的集合构成腹部CT图像光线投射结果图像。4.根据权利要求1所述的基于GPU的腹部脏器三维可视化的方法，其特征在于:所述步骤3.3对体数据进行任意平面切割，将非感兴趣区域从体数据中切除，对切割后的体数据进行光线投射，得到切割后体数据的光线投...

【专利技术属性】
技术研发人员：姜慧研，项飞，邹坤，
申请(专利权)人：东北大学，
类型：发明
国别省市：辽宁;21