IVUS三维重建总结.docx

1.主要研究内容CAG图像和预处理，CAG二维信息的提取，血管中心线的三维重建，IVUS图像的图像预处理，内外膜分割，三维血管中心线与IVUS图像融合。1.1. CAG图像预处理冠脉造影图像属于高噪声，低信噪比图像，又因为下一阶段是要进行边缘检测，所以需要考虑减少噪声并增强边界的滤波器。图像预处理方法考虑的有高斯滤波，中值滤波，灰度拉伸，tophat变换。1.1.1.高斯滤波高斯滤波器是图像处理中最有效常用的去噪方法，能够有效的去除图像中的服从正太分布的随机噪声。高斯函数具有旋转对称性，因此处理后图像不会出现偏差，并且高斯函数有可分性，可以简化运算。因此高斯平滑对冠脉造影图像有一定的去噪效果。1.1.2.中值滤波中值滤波是一种非线性空间滤波器，其对处理椒盐噪声非常有效。因此首先选用中值滤波的方法来去除这种噪声。该滤波方法所选的模板越大则平滑去噪的效果越好。它的优点是一方面它的减少噪声和增强边缘效果比较明显，而另一方面相比其他形式滤波处理速度较快。缺点是去噪的同时也会平滑掉一些有用的细节。1.1.3.灰度拉伸由于造影图像在造影过程中会受到客观条件的限制，使得图像出现一些细节上的问题。灰度拉伸处理正是针对造影图像的这一缺点进行补偿，通过对造影图像的灰度进行拉伸，把灰度区间[a, b]内的像素映射到[c, d], 灰度拉伸的目的是增强图像的对比度和提高其质量。1.1.4.Tophat变换对结构噪声应用形态学TopHat变换消除图像背景的不均匀性。冠脉造影图像中存在多种噪声因此需要取出这些不同的背景噪声。除了背景噪声外，造影胶片还常常出现曝光不均匀的现象，使得背景不均匀。Top-Hat变换方法是形态学方法中的一种，在图像对比度增强、亮度调整中有很好的应用。1.2.CAG二维信息的提取在IVUS获得的图像检查时获得的造影图像，分为两种，一种是导管可见，一种是导管不可见。若导管不可见一般把血管路径作为导管回撤路径。二维信息提取方法有很多，目前为止有主要包括：Kristiansen和Wunderlich通过自动阈值法进行边缘提取；Yim和Cebral结合加权平均值和导数对边缘求解Ritchings和Colchestert采用模式识别法判断血管狭窄；Thackry和Nelson利用形态学法对血管进行分割；Liu和Sun利用自适应跟踪法对造影图像的血管轨迹提取；Park，Klein和Chakraboty采用动态模型方法，通过初始勾划血管中心线，并构造优化能量函数，在求解函数最小值过程中使动态轮廓模型逼进中心线，从而实现血管骨架提取。图1 导管与血管中心线示意图血管是一个弯曲的空间管状系统，所以有时候用骨架来表示血管中心线，如图1。血管二维骨架提取目前有直接法和间接法。直接法是利用造影图像中血管灰度特性直接提取血管骨架，人工参与比较多，而且有很多折线。间接法是图像的细化问题，获得感兴趣血管区域后，直接利用图像的细化获得血管骨架。导引丝提取是直接对未注入造影剂的造影图像利用图像细化提取，如果是提取血管骨架要注意前后冠脉图像差，剔除导引丝和导管，另外要注意两幅冠脉图像需要是同一时间的图像，即位于ECG信号的同一相位。不过对于血管骨架即理解为导管骨架的只用提取一个，但是在三维重建中误差会增大。有的文献采用snake模型获得血管中心线，首先获得血管大概中心位置和和形状后，利用适当的能量函数，使曲线的能量值达到最小，最终停留在最优位置获得血管中心线。另外snake模型方法还经常用于IVUS内外膜提取时。还有的利用动态规划和小波变换等方法进行提取的。1.3.血管中心线的三维重建1.3.1.任意点三维坐标求取三维重建一般利用两个不同角度的冠脉造影图像，在二维信息提取的基础上，寻找匹配的点对。在进行三维重建之前，需要了解X射线血管造影系统的成像构造，图2显示了造影机的几个运动方向：两个旋转和两个平移。在从两个角度分别拍摄CAG图像时，它们之间的成像示意图是图3。图2 造影机成像运动示意图图3 造影角度示意图单面冠脉造影系统在两个造影角度下的成像如图3所示。S1和S2是X射线源，投影面坐标系和所在平面分别代表两次造影过程中的图像平面分别代表两次造影过程中的图像平面，即图像A平面和图像B平面。和则分别为两个X射线源到其对应的平面的距离。和分别为两个不同平面所对应的中心。在中，可以设空间一点在A上的所对应的称之为投影点的坐标为；在平面坐标中，点P在图像B上的投影点的坐标为，和分别为LAO/RAO角度（即造影的左右角度）。和分别表示CRAN/CAUD角度。造影系统的四个坐标系统之间的几何变换关系为刚体运动，即可通过旋转运动和平移运动实现任意坐标系之间的变换，根据透视成影几何关系可以得到如下的关系式： （1） （2）实现坐标到坐标系的变换可以经过旋转和平移完成，使轴和重