冠脉造影图像心血管边缘提取

一开题背景近年来，图像成像技术为医学诊断提供了丰富和准确的信息，对疾病的无创伤精确诊断起了积极的作用。然而，人们在获取和传输数字图像的同时，难免于图像数据被外界噪声所污染，妨碍了人们对图像信息的理解，由此，图像去噪技术应运而生。图像去噪，即在尽可能地不损失原图像细节的前提下，去除图像中无关的噪点，现有的图像去噪方法很多，常用的有1.均值滤波器2.低通滤波器。3中值滤波器。为了研究方便，我们的方法是人工地给原图像添加噪声，主要是随机噪声和椒盐噪声。另外，图像边缘反映了图像中最有价值的信息，对边缘的检测是图像处理和计算机视觉中最重要，最经典的课题之一。对医学图像而言，常需要通过边缘的提取来确定病灶的大小、脏器的运动情况、血球流动状态等。但由于多种因素的影响，医学图像边缘往往不也很清晰，这直接影响医务人员的诊断和治疗。因此，图像边缘和轮廓特征的检测和提取方法，一直是医学图像处理与分析技术中的研究热点。为方便医学影像工作者的使用，本文在对图像平滑（去噪）及对图像边缘检测技术进行回顾的基础上，分析讨论了各种方法的优势与缺陷。二课题目的在学习了医学成像系统和医学图像处理等课程的基础上，以冠脉造影图像为例，利用Matlab软件，学习医学图像的图像增强和边缘提取等。三课题研究的主要内容.分别利用中值滤波和高斯平滑去除图像中的椒盐噪声和随机噪声，分析不同的模板大小对中值滤波的影响，以及尺度（方差）对图像平滑的影响。.对造影图像分别用灰度梯度算子、Roberts交叉算子、Sobel算子、Kirsch算子、Laplace算子、Marr算子进行边缘提取，并对结果进行讨论，并研究阈值选取对边缘检测的影响。四原理和方法1图像平滑图像平滑的概念与分类。图像平滑主要是为了消除被污染图像中的噪声。图像的简单平滑是图像增强处理中最基本的方法之一，它利用卷积运算对图像邻域的像素灰度进行平均化，从而达到减少图像中杂点的影像，降低图像对比度的目的。在空间域方法中，主要有线性平滑，非线性平滑，自适应平滑；在频率域，由于噪声频谱通常在高频段，因此主要是选择合适的低通滤波器对其频谱成分进行调整，然后经逆傅里叶变换得到平滑图像（由于噪声频谱通常在高频段）。中值滤波⑴中值滤波是一种非线性的信号处理方法，与其对应的中值滤波器当然也就是一种非线性的滤波器。中值滤波器在1971年由J.w.Jukey首先提出并应用在一维信号处理技术中，后来被二维图像信号处理技术所引用。中值滤波在一定的条件下可以克服线性滤波器如均值滤波带来的图像细节模糊，而且对滤除脉冲干扰（椒盐噪声即脉冲干扰）及图像扫描噪声最为有效。但是对于一些细节多，特别是点、线、尖顶细节多的图像不宜采用中值滤波。另外，中值滤波很容易自适应化，从而可以进一步提高其滤波性能。因此，它非常适用于一些线性滤波器无法胜任的数字图像处理应用场合。中值滤波一般采用一个含有奇数个点的滑动窗口，将窗口中的各点灰度值的中值来替代指定点（一般是窗口的中心点）的灰度值。对于奇数个元素，中值是指按大小排序后，中间的数值；对于偶数个元素，中值是指排序后中间两个元素灰度值的平均值。高斯滤波高斯滤波器是一种线性滤波模板，模板各位置的权重是根据高斯分布函数确定的。由于高斯滤波器的二维可分性（X轴与Y轴方向进行高斯滤波互不干扰），我们可以用两个一维高斯滤波器的连续卷积来实现一个二维高斯滤波器：（1）对图像使用一维高斯卷积模板，在一个方向上进行滤波（例如水平方向）。（2）转置图像；（3）对转置以后的图像使用同一个高斯卷积模板，在垂直方向进行滤波。（4）将图像转置回原来的位置，我们就得到了经过二维滤波的图像。一维高斯卷积模板可以由二项式展开的系数来模拟a)3X3模板±X4X11b)5X5模板146412边缘检测边缘检测的定义图像的边缘是图像的最基本特征。所谓边缘就是指其周围像素灰度有阶跃变化或屋顶变化的那些像素的集合。边缘广泛存在于物体与背景之间、物体与物体之间、基元与基元之间。因此，它是图像分割所依赖的重要特征。物体的边缘是由灰度不连续性所反映的。经典的边缘提取方法是考察图像的每个像素在某个邻域内灰度的变化，利用边缘邻近一阶或二阶方向导数变化规律。用简单的方法检测边缘，这种方法称为边缘检测局部算子法。几种常用的边缘检测算子。边缘检测算子检查每个像素的邻域并对灰度变化率进行量化，也包括方向的确定。大多数使用基于方向导数掩没求卷积的方法。灰度梯度算子当图像微分运算后，每个像素处有一个梯度值反映该点的边缘强度。f(x,y).；:f(x,y).'、f(x,y)二।jexcy模板:-11-1模板:-11-1差分：Gx=f(x+1,y)-f(x,y);Gy=f(x,y+1)-f(x,y)为了检测边缘点，选取适当的阈值T,对梯度图像进行二值化。g(x,y)=1Grad(x,y)＞二T二0其他这样，形成了一副边缘二值图像g(x,y)。Roberts交叉算子Roberts边缘检测算子是一种利用局部差分算子寻找边缘的算子。if(x,y)尸max(f(x,y)-f(x1,y1)|,|f(x,y1)-f(x1,y)|)分fx'=|f(x1,y1)-f(x,y)|".fy'=|f(x1,y)-f(x,y1)|Sobel算子-101

下图所示的两个卷积-202点都用这两个核做卷积,-1~0-1一个对水平边缘响应最运算结果是一幅边缘幅度图像。下图所示的两个卷积-202点都用这两个核做卷积,-1~0-1一个对水平边缘响应最运算结果是一幅边缘幅度图像。模板：Krisch算子下图所示的8个卷积核组成了Krisch边缘算子。图像中的每个点都用8个掩模进行卷积，每个掩模都对某个特定边缘方向作出最大响应，所有8个方向中最大值作为边缘边缘幅度图像的输出。最大响应掩模的序号构成了边缘方向的编码。-31-3---3-31-3---30-3+5+5-+5+5-3-3+50-；3+5-3-3+5.+5-「5士口□&3-3□Laplace算子Laplace算子是对二维函数进行运算的二阶导数算子。通常使用的拉普拉斯由于拉普拉斯算子是一个二阶导数，它将在边缘处产生一个陡峭的零交叉，对噪声非常敏感，另外。，它常产生双像素宽的边缘，因而不能提供图像边缘的方向信息，所以很少直接用于边缘检测Marr算子Marr算子是在Laplace算子的基础上实现的。Lapace对噪声比较敏感，为了消除噪声的影响，可以对图像先进行平滑然后再运用Laplace算子。根据人眼视觉成像原理，图像中任一个像素点所对应的场景点周围的点对该点的光强的贡献成正态分布。因此，平滑函数选择具有正态分布的高斯函数h(x,y)oG(x,y)=f(x,y)*h(x,y)得到平?t图像G(x,y)后再用拉普拉斯算子处理。五实验步骤1流程图心脉造影图像心脉造影图像的影响分析不同模板的影响分析不同尺度(方差)的影响分析不同模板的影响分析不同尺度(方差)的影响图1实验流程图2步骤说明1对所给冠脉造影图像分别加入随机噪声和高斯噪声。2分别用3*3、7*7、13*13、19*19、25*25的模板对随机噪声图像和椒盐噪声图像