专利名称：利用vtk进行肠道三维重建的方法
技术领域：
本发明涉及利用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 = (Hd1)只可能有下面3种情况
①点、= (WZjw)在轮廓Ck+1内，但不在Ck+1上，此时执行4^)05，其中壬力膨胀，B为结构算子。②点巧既不在轮廓(；+1上，也不在轮廓内，即该点和点不属于同一密度物质，此时，执行，其中表示腐蚀。③点Vi =在轮廓Ck+1上，即lU/Jfc+i) e Cm ,此时无需任何操作。综合上述3种情况，可以表示为根据式(I)、式(2)对任意的两幅断层图像进行操作，最终均能得到两幅完全一样的插值图像。此时，在断层图像Slri和Sk+1之间已经产生了一系列新的图像，找到与Sk位置对应的新图像，提取出它的每个密度物质轮廓Ck，这样就构成了插值图像Sk的轮廓。
权利要求
1.一种利用VTK (Visualization Toolkit可视化工具包)进行肠道三维重建的方法, 其特征在于具体操作步骤如下(1)对含有噪声的CT数据进行预处理，对CT图像进行中值滤波，从而滤除由于影像设备中各电子间的随机扰动和周围环境产生的噪声；⑵对CT图像采用边界分割Live-Wire (火线)算法和区域生长相结合的分割算法分割出肠道部分的CT图像；(3)对分割得到的肠道CT图像进行数学形态学的改进的最近邻域插值，提高CT图像的层间分辨率；⑷设置VTK平台上的重建环境；(5)执行肠道面重建算法。
2.根据权利要求I所述的利用VTK进行肠道三维重建的方法，其特征在于所述步骤 ⑴中对CT图像进行中值滤波，其操作方法对于CT图像上的每一个点，选取一个含有奇数个像素的窗口，对该窗口所覆盖像素按大小排序，将处在灰度序列中间的那个灰度值作为中值，用它代替窗口中心对应的像素的原灰度值。
3.根据权利要求I所述的利用VTK进行肠道三维重建的方法，其特征在于所述步骤 ⑵中对CT图像采用边界分割Live-Wire算法和区域生长相结合的分割算法分割出肠道部分的CT图像，其具体步骤如下①用Live-Wire算法在CT图像上选取边界点，作为肠道的边界；②在边界线内选取种子点；③根据设定的条件进行区域生长，从而找到完整的肠道。
4.根据权利要求I所述的利用VTK进行肠道三维重建的方法，其特征在于所述步骤 ⑶中对分割得到的肠CT道图像进行数学形态学的改进的最近邻域插值，提高图像的层间分辨率，其具体步骤如下①利用数学形态学的方法，得到被插值图像区域的轮廓；②求出插值图像Sk上所有点在前一层图像Sh和后一层图像Sk+1上的对应点；③求出插值图像Sk中各点的灰度值。
5.根据权利要求I所述的利用VTK进行肠道三维重建的方法，其特征在于所述步骤 ⑷中设置VTK平台上的重建环境，其具体步骤如下设置肠道RGB颜色为1，O. 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)。
6.根据权利要求I所述的利用VTK进行肠道三维重建的方法，其特征在于所述步骤(5)中执行肠道面重建算法，其具体步骤如下①逐点判断各个体素的面是否为二义性面，在计算得到肠道CT体数据与等值面的交点位置的基础上，判断是否为二义性面，包括体素与等值面交点坐标的计算和体素各个顶点的状态判断；②消除二义性面用对边交点的连线代替渐进线，计算连线交点O处的坐标，并判断交点处的状态值，从而判断连接方式，消除二义性；③拼接三角形面片，得到重建的三维模型，并进行平滑滤波从而得到表面光滑的三维肠道；④利用VTK的渲染机制，渲染三维肠道面重建结果。
全文摘要
本发明涉及一种利用VTK(VisualizationToolkit可视化工具包)进行肠道三维重建的方法。本方法直接基于VTK进行三维重建。首先，对含有噪声的原始人体CT扫描数据进行滤波和体数据的分割，得到肠道组织的CT数据；为了得到更好的重建效果，对数据进行插值，采用形态学的方法对体数据进行插值；然后根据MarchingCubes面重建理论，设计出高效解决面重建二义性的算法。在面重建时，根据肠道的特性，设置好颜色、散射光及散射光强度。实验结果表明，采用本发明的三维肠道重建算法能逼真的重建出三维肠道，且具有较快的重建速度。
文档编号G06T17/00GK102592311SQ20121000809
公开日2012年7月18日 申请日期2012年1月12日 优先权日2011年3月10日
发明者万旺根, 余小清, 王亚男, 王蕊, 马艳茹 申请人:上海大学