利用VTK进行肠道三维重建的方法技术

本发明专利技术涉及一种利用VTK（VisualizationToolkit可视化工具包）进行肠道三维重建的方法。本方法直接基于VTK进行三维重建。首先，对含有噪声的原始人体CT扫描数据进行滤波和体数据的分割，得到肠道组织的CT数据；为了得到更好的重建效果，对数据进行插值，采用形态学的方法对体数据进行插值；然后根据MarchingCubes面重建理论，设计出高效解决面重建二义性的算法。在面重建时，根据肠道的特性，设置好颜色、散射光及散射光强度。实验结果表明，采用本发明专利技术的三维肠道重建算法能逼真的重建出三维肠道，且具有较快的重建速度。

【技术实现步骤摘要】

本专利技术涉及利用VTK (Visualization Toolkit可视化工具包)进行肠道三维重建的方法，尤其是涉及利用VTK对肠道的CT图像进行三维重建的方法。
技术介绍
三维肠道重建就是用一组肠道的CT断层图像，通过某种算法将其转化为具有直观立体效果的肠道图像，更好的展现肠道的真实三维形态。三维肠道重建是科学诊断肠道疾病的基础，也是肠道疾病诊断治疗需要解决的关键技术之一。大肠癌是我国常见的恶性肿瘤之一，患者平均5年生存率为20%_25%，我国肠癌发病率的增速是世界平均水平的两倍，达到年均4%。发病率仅次于肺癌和肝癌，列第三位，已经成为癌症中的首要杀手之一。然而，在医疗诊断中，观察病人的一组二维断层图像CT、MRI 是医生诊断病情的常规方式。但是，要准确地确定病变体的空间位置、大小、几何形状以及与周围生物组织之间的空间关系，医生仅通过二维断层图像，在头脑中想象真实的三维模型是十分困难的，从而影响了疾病诊治的精确性。因此，迫切需要一种行之有效的工具来完成对人体器官、软组织和病变的三维重建和显示。医学图像三维重建就是辅助医生对病变体的周围组织进行分析和显示的有效工具，它极大地提高了医疗诊断的准确性和科学性。
技术实现思路
本专利技术的目的在于提供一种利用VTK进行肠道三维重建的方法，实现利用CT图像进行逼真的三维肠道重建，提高医疗诊断的准确性和科学性。为达到上述专利技术目的，本专利技术的构思是首先对原始的腹部CT图像进行中值滤波，去除由于影像设备中各电子间的随机扰动和周围环境产生的噪声，再用Live-Wire和区域生长相结合的分割方法分割出肠道的CT图像，接着对肠道图像进行基于数学形态学的改进的最近邻域插值算法对肠道图像进行插值，在VTK平台上构建环境，以此重建三维肠道。根据上述专利技术构思，本专利技术采用下述技术方案一种利用VTK进行肠道三维重建的方法，其特征在于具体操作步骤如下(1)对含有噪声的CT数据进行预处理，对CT图像进行中值滤波，从而滤除由于影像设备中各电子间的随机扰动和周围环境产生的噪声；⑵对CT图像采用边界分割Live-Wire (火线)方法和区域生长相结合的分割方法分割出肠道部分的CT图像；(3)对分割得到的肠道CT图像进行数学形态学的改进的最近邻域插值，提高CT图像的层间分辨率；⑷设置VTK平台上的重建环境；(5)执行肠道面重建算法。所述步骤⑴中对CT图像进行中值滤波，其操作方法对于CT图像上的每一个点，选取一个含有奇数个像素的窗口，对该窗口所覆盖像素按大小排序，将处在灰度序列中间的那个灰度值作为中值，用它代替窗口中心对应的像素的原灰度值。所述步骤⑵中对CT图像采用边界分割Live-Wire算法和区域生长相结合的分割算法分割出肠道部分的CT图像，其具体步骤如下①用Live-Wire算法在CT图像上选取边界点，作为肠道的边界；②在边界线内选取种子点；③根据设定的条件进行区域生长，从而找到完整的肠道。所述步骤⑶中对分割得到的肠CT道图像进行数学形态学的改进的最近邻域插值，提高图像的层间分辨率，其具体步骤如下①利用数学形态学的方法，得到被插值图像区域的轮廓；②求出插值图像Sk上所有点在前一层图像Sk_i和后一层图像Sk+1上的对应点；③求出插值图像Sk中各点的灰度值。所述步骤⑷中设置VTK平台上的重建环境，其具体步骤如下设置肠道RGB颜色为1，0. 49，O. 25，反射系数O. 3，反射光强20，在VTK中的设置为 设置肠道颜色intestine->GetProperty O->SetDiffuseColor (1，O. 49，O. 25)； 设置反射系数intestine_>GetProperty O->SetSpecular (O. 3)；设置反射光强intestine_>GetProperty O->SetSpecularPower (20)。所述步骤(5)中执行肠道面重建算法，其具体步骤如下①逐点判断各个体素的面是否为二义性面，在计算得到肠道CT体数据与等值面的交点位置的基础上，判断是否为二义性面，包括体素与等值面交点坐标的计算和体素各个顶点的状态判断；②消除二义性面用对边交点的连线代替渐进线，计算连线交点O处的坐标，并判断交点处的状态值，从而判断连接方式，消除二义性；③拼接三角形面片，得到重建的三维模型，并进行平滑滤波从而得到表面光滑的三维肠道；④利用VTK的渲染机制，渲染三维肠道面重建结果。本专利技术与现有技术相比较，具有如下显而易见的突出实质性特征和显著优点先对CT图像进行分割，大大减小了冗余数据量；面重建阶段的二义性处理方法计算简单，极大地减小了计算量，提高了计算速度；对重建环境的设置使重建结果更接近真实肠道；在可视化工具VTK平台上实现三维重建，可充分发挥VTK的强大功能，而且有助于后续三维肠道内漫游和肠道疾病的检测。附图说明图I是本专利技术利用VTK进行肠道三维重建的方法的流程图。图2是最近邻域插值示意图。图3是重建出的逼真三维肠道效果图。具体实施例方式本专利技术的一个优选实例结合附图详细说明如下参见图1，本专利技术利用VTK进行肠道三维重建的方法的一个优选实例结合附图说明如下本专利技术的实施利用解放军海军总医院32排CT机中的一套889张人体胸腔CT扫描数据，分割出包含肠道部分的356张肠道CT 数据，图像大小512X512像素，层间距0.62mm。利用VTK进行肠道三维重建方法共分为五步，其具体操作步骤如下(I)对含有噪声的CT数据进行预处理，对CT图像进行中值滤波，从而滤除由于影像设备中各电子间的随机扰动和周围环境产生的噪声；其实施方法如下在图像上选取一个含有奇数个像素的窗口，对该窗口所覆盖像素按大小排序，处在灰度序列中间的那个灰度值为中值，用它代替窗口中心对应像素的原灰度值。(2)对CT图像采用边界分割Live-Wire方法和区域生长相结合的分割方法分割出肠道部分的CT图像，其具体如下人体腹部的CT扫描数据是由几百张CT图像组成的，要得到肠道的完整数据就要对这几百张图片整体分割，由于基于边界的分割不能用于体数据，而基于区域的分割虽然能用于体数据，但是它会造成过度分割，即将图像分割成过多的区域。本方法采用边界分割与区域分割相结合的分割方法，充分利用两种分割的优点更好地分割出肠道，具体方法如下 ①用Live-Wire算法在图像上选取边界点，作为肠道的边界。②在边界线内选取种子点，③根据设定的条件进行区域生长，从而找到完整的肠道。(3)对分割得到的肠道图像进行数学形态学的改进的最近邻域插值，提高图像的层间分辨率；其具体如下给定两幅断层图像Slri和Sk+1，为了通过插值求出它们之间的图像Sk，用集合S= (V，F) 来表示已知的所有断层图像，V为所有点的集合,其元素v=(x, y, Z)表示点在空间的坐标, 函数f :V — Y定义了任一点的灰度值。第一步，利用数学形态学的方法，得到被插值图像每个区域的轮廓。由步骤⑵得到2幅断层图像，接着要确定已知两幅断层图像上肠道的轮廓，分别用Clri和ck+1表示分割得到的肠道区域的边界。对任意位于图像Slri某个轮廓上的点6 = ( ,//, .-!)，即(b&u) € Cm来说，对应于图像Sk+1上的点5 = 本文档来自技高网...

【技术保护点】

【技术特征摘要】
...

【专利技术属性】
技术研发人员：万旺根，王蕊，余小清，王亚男，马艳茹，
申请(专利权)人：上海大学，
类型：发明
国别省市：