专利名称：混凝土ct图像三维重构方法
技术领域：
本发明涉及混凝土 CT图像处理技术，尤其涉及一种根据对象的二维断层图像数据重构其三维图像的方法。
背景技术：
混凝土是一种工程应用相当广泛的不均匀材料。由粗骨料、细骨料、水泥水化物、 未水化水泥颗粒、孔隙及裂纹等组成的非均质复介材料。研究各种空间尺度缺陷对混凝土材料的变形、破坏和稳定性所起的作用，或者作用的大小，是目前混凝土材料研究所面临的最深刻而又最艰难的课题之一。近些年X射线混凝土 CT试验成为混凝土、岩石材料细观破裂过程的热点研究课题，X射线CT观察混凝土裂纹演化过程最大优点在于它的无损探测性能和较高分辨率。对于实验后的得到混凝土的CT图像数据分析，采用直接图像分析法。直接图像分析是指直接对混凝土受力后的结构变化通过图像分析。然而混凝土 CT图像中的许多萌生细微裂纹都无法直接从图像中直接准确的判断是否存在；另外由于CT机的扫描原理，所得到混凝土CT图像只是横切面图片，对于混凝土体柱受力后图像分析，无法直接进行，要想通过二维断层图像认识真实物体，就必须完成从二维图像到三维图像的恢复工作，需要对混凝土 CT图像进行三维重构。三维重建结果得到的混凝土真实的细观结构模型，可以与混凝土细观数值模拟结果进行试验比较，并校正数值分析模型。体绘制算法用于完成三维结果图像的生成和显示。光线投影算法是最常用的以图像空间为序的典型直接体绘制算法。利用光线投射法绘制出的图像质量很高、效果很好，但是此方法存在很大的问题， 那就是处理的数据量比较大，对内存的要求比较高，由于要对每条光线进行重采样计算，所以绘制速度比较慢。如何进行算法的改进以加快绘制速度已经成为目前光线投射算法研究中最需要解决的问题，也是体绘制技术中最主要的研究方向。目前，人们最常用的加速技术主要从两个方面考虑一方面从图像空间的相关性着手使发射出的光线数目尽量减少，另一方面从对象空间的相关性着手，省略灰度、梯度方向相同或相近体素的采样，从而使获得的采样点数目尽量减少。因此，提供一种在不影响绘制图像质量的情况下，从这两方面考虑对传统的光线投射算法进行改进，使绘制速度得到很大提高的方法成目前亟待解决的技术问题。
发明内容
针对现有光线投射算法中存在的问题和不足，本发明的目的在于提供一种混凝土图像重构方法，利用对象的二维断层图像数据重构其三维图像。运用三维重构可视化技术， 对混凝土 CT图像进行立体重构，真实反映出混凝土柱应力后结构的变化。为分析混凝土的裂纹扩展及破坏形态，揭示出混凝土材料的损伤破裂机理提供了一种重要的辅助分析手段和方法。本发明的技术方案是一种混凝土 CT图像三维重构的方法，包括步骤如下
I)采用医用CT或工业CT扫描混凝土试件断面，获得混凝土静力压缩CT 二维断层图像数据；2)运用三维重构可视化技术对混凝土静力压缩CT 二维断层图像进行立体重构， 生成混凝土 CT三维显示图像；3)对混凝土 CT三维显示图像进行三维优化，最终完成混凝土 CT图像三维重构。本发明进一步的特征在于所述步骤2)中，运用三维重构可视化技术对混凝土 CT 二维断层图像进行立体重构，采用光线投影算法按照设定视线的方向，从混凝土静力压缩CT 二维屏幕上的每个像素点出发，发出一条射线，沿着该射线按照一定的步长进行等距采样；对距离某个采样点最近的八个体素的属性值做三线性插值运算从而计算出，得到该采样点的属性值(颜色值和不透明度值)，对得到的每个采样点的属性值进行逐点合成，从而得到所对应的屏幕上该像素点的颜色值和不透明度值，最终生成混凝土 CT三维显示图像。上述对得到的每个采样点的属性值按照由后向前的图像合成和由前向后的图像合成的图像合成算法进行逐点合成。所述步骤3)中，图像进行三维重建化过程包括下述步骤①数据预处理优化使用包围盒技术对三维空间数据场进行裁剪，在平行光线的投影下，利用规则网格数据场特点，使用Graham方法实现用一个能把三维体数据场中的有用信息都包围起来的最小长方体包围盒，把目标体数据限定在其中；把包围盒投影到成像屏幕上，形成一个投影多边形；对所得到的最优投影多边形进行一次扫描转换实现投射光线的生成；②投射光线体素化及重采样优化通过图像空间坐标计算转换为物体空间坐标， 接着在物体空间坐标系OXYZ中利用Bresenham方法推广到三维，对投射光线进行体素化和求交处理。所述步骤①中，用Graham求解最小包围盒凸壳的算法过程如下a)先计算平面凸集中所有点的y坐标，记录y值最小点，之后把它与凸集中所有剩余点进行连接形成线段，接着计算所得到的线段与水平线之间的夹角；然后对这些点先按照夹角大小进行排序，如果夹角相同，就按照与y值最小点的距离排序，依次连接得到的序列点集pl，p2,…，pn,从而形成一个多边形；b)通过判断某两个点是不是在某一条指定线段的同侧，把p3到pn-Ι中不是凸壳顶点的点删除；c)最后把经过删除所得到的凸壳顶点进行输出。所述步骤②中，投射光线体素化及重采样优化，坐标转换是通过下述迭代关系式实现的ρ2' ·χ = p2. x+r0 ；p2' . y = p2. y+r3 ；p2' . z = p2. z+r6其中，(p2. x, p2. y，p2. z)是点pi在物体空间坐标系中对应的点pl'的坐标；rQ、 r3> r6分别表示像素在X，y, z轴上的物体空间位置。与现有技术相比，本发明的积极效果是本发明所生成的混凝土 CT三维图像提供立体的混凝土骨料、砂浆、孔洞和裂纹区域结构，为混凝土细观破裂研究提供分析研究观测平台。本发明方法实时性强、运算效率高，能清晰的分辨骨料、砂浆及孔洞。
本发明具有下述特点I)混凝土图像重构方法是在光线投射算法基础上提出了一种加速的光线投射算法，极大的提高了传统光线投射法图像绘制的速度，对图像的绘制质量几乎没有影响。并且，由于大部分的加速工作都是在预处理阶段完成的，从而使图像的交互速度也有了很大的提闻。2)混凝土三维重建结果得到了混凝土真实的细观结构，可以与混凝土细观数值模拟结果进行试验比较，校正数值分析模型。3)三维重建技术为混凝土 CT图像立体细观研究提供了一定的观测技术平台。4)混凝土真实细观结构的三维重建结果可以作为导入到大型工程软件(如 ANSYS)的前期技术基础，从而可以建立起能更加准确表征混凝土非均质特性的三维有限元模型，为土木工程学者进一步研究混凝土材料的力学特性和细观破坏机理提供更好的支持。
下面结合附图对本发明作进一步的详细说明。图I为本发明方法步骤流程图。图2(a)-图2(e)为本发明Graham方法求解投影多边形凸壳过程；图3为本发明坐标转换成像平面相邻点；图4为本发明排序两交点；图5为本发明光栅化直线方法Bresenham方法；图6为本发明体素化投射光线。
具体实施例方式本发明的混凝土图像重构方法，如图I所示，该方法包括第一步，采用医用CT或工业CT扫描混凝土试件断面，获得混凝土静力压缩CT图像二维断层图像数据；第二步，运用三维重构可视化技术对混凝土静力压缩CT 二维断层图像进行立体重构，生成三维显示图像；具体方法如下运用三维重构可视化技术对混凝土 CT 二维断层图像进行立体重构，采用光线投影算法按照设定视线的方向，从混凝土静力压缩CT 二维屏幕上的每个像素点出发，发出一条射线，沿着该射线按照一定的步长进行等距采样；对距离某个采样点最近的八个体素的属性值做三线性插值运算从而计算出，得到该采样点的属性值(颜色值和不透明度值)，对得到的每个采样点的属性值按照从后到前或从前到后的顺序进行逐点合成，从而得到所对应的屏幕上该像素点的颜色值和不透明度值，最终混凝土 CT生成三维显示图像。第三步，对混凝土 CT三维显示图像进行三维优化，从而完成混凝土 CT图像三维重构。①数据预处理优化。A.使用包围盒技术对三维空间数据场进行必要的裁剪，这样可以避免很多不必要的数据计算，从而提供绘制速度。
运用包围盒技术对其进行裁剪的过程为在平行光线的投影下，利用规则网格数据场的一些特点，用一个能把三维体数据场中的有用信息都包围起来的最小长方体包围盒，把目标体数据限定在其中。这样就可以利用包围盒对数据场进行有效的裁剪，从而极大的缩小了数据场的规模。我们假设OUVW是图像空间坐标系，OXYZ是物体空间坐标系。假如用Λχ，Δγ, Λ Z表示三维数据场中X，Y，Z三个坐标轴上每个体素的步长间隔，对体素进行遍历就可以确定对应的数据场的最小包围盒的长、宽、高，它们分别是(maxx_minx) * Δ X, (maxy_miny) * Δ y, (maxz_minz) * Δ z,如果体素在包围盒内部，就可以对其进行采样；如果它在包围盒外部就忽略，即认为该体素为空体元。要对体数据进行显示必须完成体视见变换，即把体数据从物体空间OXYZ转换到图像空间0UVW。但是有些体素的体视见变换是不必要的，为了提高算法速度，必须减少这些不必要的变换。把包围盒投影到成像屏幕上，形成一个投影多边形，只从投影多边形中的像素发出光线进行投射，就可以保证所发射的每条光线都能与三维空间的体数据相交。这种方法不仅能减少了一些不必要的体视见变换，同时也通过减少所包含的像素点使发出的投射光线数目与求交运算也大幅度的下降，并且还能有效的保证数据场中非空体元的遗漏。 这一过程使用Graham方法实现。如图2(a)-图2(e)所示,用Graham求解最小包围盒凸壳的算法过程如下(I)先计算平面凸集中所有点的y坐标，把y值最小的点记作pl，之后把它与凸集中所有剩余点进行连接形成线段，接着计算所得到的线段与水平线之间的夹角。然后对这些点进行词典式排序，即先按照夹角大小进行排序，如果夹角相同就按照与Pi的距离排序，依次连接得到的序列点集pl，p2，…，pn，从而形成一个多边形。因为Pl是凸壳边界的起点，且plp2和plpn与水平线的夹角分别最小和最大,所以p2和pn也肯定是凸壳的顶点。(2)通过判断某两个点是不是在某一条指定线段的同侧，把p3到pn-Ι中不是凸壳顶点的点进行删除。(3)最后把经过删除所得到的凸壳顶点进行输出。B.对所得到的最优投影多边形进行一次扫描转换实现投射光线的生成。进行扫描转换的目的就是用一些常用的转换算法按照显示器要求的模式进行图像的精确组织，从而可以用视频技术来显示所得到的光栅图像，采用活化边表和排序Y桶的有序边表算法来实现投影多边形的扫描转换，从而记录在图像空间OUVW中扫描到的有效像素点位置。其具体算法如下(I)数据准备采用位于y+1/2的中心扫描线，对于多边形的每条边求出与其相交的最高扫描线。把该多边形的边存入与这一扫描线相对应的y桶。把X交点的初值，多边形边穿过的扫描线条数Ay以及相邻扫描线之间的X增量Λχ存入链表。(2)数据的扫描转换对于每条扫描线，检查相应的y桶是否有新的边，把新的边移入活化边表。在活化边表中对应交点安X递增顺序排序，即当X1 < X2时，X1位于X2之前。从已X排序的表中
成对地取出交点。在扫描线y上激活X整数值以满足X1Sx + iS&的像素。活化边表中每条边的Ay减I。若Ay < O,从活化边表中删除这条边。计算新的X截距交点X新=X旧 + Λχ。对下一条扫描线重复上述过程。②投射光线体素化及重采样优化。这一阶段的优化是与上一节中投影多边形的扫描转换相结合一起完成的。通过投影多边形的扫描转换我们记录了在图像空间OUVW中扫描得到的有效像素点的位置，首先需要把这些坐标进行一下转换运算，即通过图像空间坐标计算转换为物体空间坐标，接着在物体空间坐标系OXYZ中利用Bresenham方法推广到三维，对投射光线进行体素化和求交处理。Α.坐标转换。用R表示旋转矩阵，T表示平移矩阵，则可以假设从图像空间坐标系到物体空间坐标系的转换矩阵是M = R*T。在坐标转换过程中可以利用坐标之间的相关性进行优化，由它们之间的相关性得到迭代关系式，使坐标转换之间的运算量减少，从而极大的提高了图像的绘制速度。在图3中可知，pl的坐标为(ul，vl，wl)，p2的坐标为(u2，v2，w2)，它们是相邻的两个像素点，都位于图像空间坐标系中的同一扫描线上，由于Pi和P2之间存在相关性，所以P2的坐标也可以表示为(pl.u+l，pl.v，pl.w)。这样就可以运用下面的迭代关系式快速的把P2的坐标转换为物体空间坐标系上的坐标ρ2' ·χ = p2. x+r0 ；p2' . y = p2. y+r3 ；p2' . z = p2. z+r6其中(p2. x，p2. y，p2. z)是点pl在物体空间坐标系中对应的点pl'的坐标。rQ、 r3> r6分别表示像素在X，y, z轴上的物体空间位置。经过这样的迭代处理，使进行坐标转换处理时的运算量比之前普通的矩阵运算有了大幅度的减少，极大的提高了绘制速度。对应处在同一扫描线上的所有像素点全都采用这样的迭代关系式来进行坐标转换的计算。按照视线方向从点P2发射一条光线，利用投射光线与平面族求交代替直接与体元求交算法，同时结合上一阶段找到的体数据场的包围盒就可以非常快速的求出此投射光线与三维体数据场的起始和终止交点，有效的剪裁了投射光线，这样就可以避免后面阶段的盲目采样，从而减少了无用采样点的数目。在求交的过程中，我们只需要考虑光线与长方体包围盒六个表面Ml，M2，……，M6的相交情况。因为入射光线方向是任意的，它与长方体包围盒相交的入点和出点有以下三种情况1.交点落在长方体包表面的内点处；2.交点落在长方体棱边的内点处；3.交点落在长方体顶点处。如果求的的两个交点重合，那么该光线对应的屏幕上的像素点颜色值设为背景颜色值；否则需要利用投射光线的方向对这两个交点进行排序，从而确定出入点。在此我们利用判定两向量同向的方法对这两个交点进行排序。在图4中，我们假设M(xm，yM，zM)、N(xn，yN，zN)是当前入射光线与长方体包围盒Box的两个交点，其中I = Oci,表示投射光线的方向，并且它不是零向量。现在对Μ，N进行排序，从而确定出入点。我们先不妨假设M是入点，N是出点，从入点到出点的向量与投射光线的向量f的关系如下I=又2(又>0)，即要断向量
是否与向量；^同向。本文采用符合函数sign(x)进行判断 I if X > O
权利要求
1.混凝土CT图像三维重构方法，其特征在于，该方法包括下述步骤第一步，采用医用CT或工业CT扫描混凝土试件断面，获得混凝土静力压缩CT 二维断层图像数据；第二步，运用三维重构可视化技术对混凝土静力压缩CT 二维断层图像进行立体重构， 生成混凝土 CT三维显示图像；第三步，对混凝土 CT三维显示图像进行三维优化，最终完成混凝土 CT图像三维重构。
2.根据权利要求I所述的混凝土CT图像三维重构方法，其特征在于，所述第二步骤中， 运用三维重构可视化技术对混凝土 CT 二维断层图像进行立体重构，采用光线投影算法按照设定视线的方向，从混凝土静力压缩CT 二维屏幕上的每个像素点出发，发出一条射线， 沿着该射线按照一定的步长进行等距采样；对距离某个采样点最近的八个体素的属性值做三线性插值运算，得到该采样点的属性值，对得到的每个采样点的属性值进行逐点合成，从而得到所对应的屏幕上该像素点的颜色值和不透明度值，最终生成混凝土 CT三维显示图像。
3.根据权利要求2所述的一种混凝土CT图像三维重构方法，其特征在于，所述对得到的每个采样点的属性值按照由后向前的图像合成和由前向后的图像合成的图像合成算法进行逐点合成。
4.根据权利要求I所述的混凝土CT图像三维重构方法，其特征在于，所述第三步骤中， 对混凝土 CT三维显示图像进行三维优化过程包括下述步骤①数据预处理优化使用包围盒技术对三维空间数据场进行裁剪，在平行光线的投影下，利用规则网格数据场特点，使用Graham方法实现用一个能把三维体数据场中的有用信息都包围起来的最小长方体包围盒，把目标体数据限定在其中；把包围盒投影到成像屏幕上，形成一个投影多边形；对所得到的最优投影多边形进行一次扫描转换实现投射光线的生成；②投射光线体素化及重采样优化通过图像空间坐标计算转换为物体空间坐标，接着在物体空间坐标系OXYZ中利用Bresenham方法推广到三维，对投射光线进行体素化和求交处理。
5.根据权利要求4所述的一种混凝土CT图像三维重构方法，其特征在于，所述步骤① 中,用Graham求解最小包围盒凸壳的算法过程如下a)先计算平面凸集中所有点的y坐标，记录y值最小点，之后把它与凸集中所有剩余点进行连接形成线段，接着计算所得到的线段与水平线之间的夹角；然后对这些点先按照夹角大小进行排序，如果夹角相同，就按照与y值最小点的距离排序，依次连接得到的序列点集pl,p2,…，pn,从而形成一个多边形；b)通过判断某两个点是不是在某一条指定线段的同侧，把p3到pn-1中不是凸壳顶点的点删除；c)最后把经过删除所得到的凸壳顶点进行输出。
6.根据权利要求4所述的一种混凝土CT图像三维重构方法，其特征在于，所述步骤② 中，投射光线体素化及重采样优化，坐标转换是通过下述迭代关系式实现的p2, . X = ρ2· x+r0 ;p2, · y = ρ2· y+r3 ;p2, · z = p2. z+r6其中，(p2. x，p2. y，p2. z)是点pi在物体空间坐标系中对应的点pi,的坐标；Ivivr6分别表示像素在X，1，z轴上的物体空间坐标位置。
全文摘要
本发明公开了一种混凝土CT图像三维重建的方法，首先用CT机对混凝土试样进行断层扫描；其次光线投影算法对混凝土CT图像进行三维重建，得到骨料、砂浆、孔洞的真实图像；最后利用可视化工具VTK对混凝土CT图像进行三维立体重构。本发明的积极效果是实时性强、运算效率高、能清晰可靠的分辨骨料、砂浆及孔洞。三维重建结果得到了混凝土真实的细观结构模型，可以与混凝土细观数值模拟结果进行试验比较，并校正数值分析模型。
文档编号G06T17/00GK102609980SQ20121001595
公开日2012年7月25日 申请日期2012年1月18日 优先权日2012年1月18日
发明者党发宁, 徐胜军, 李昌华, 王宏波, 赵亮, 陈登峰 申请人:西安万方工业自动化有限公司, 西安建筑科技大学