冠脉造影图像心血管边缘提取资料.docx

1、冠脉造影图像心血管边缘提取资料目录1、开题背景 32、课题目的 33、课题主要内容 34、原理和方法 34.1 边缘检测技术概述 34.2 灰度梯度算子 44.3 Roberts交叉算子 44.4 Sobel算子 44.5 Kirsch算子 44.6 Laplace算子 44.7 Marr算子 55、步骤 56、结果 66.1分别利用中值滤波和高斯平滑去除图像中椒盐噪声和随机噪声 66.2用灰度梯度算子进行边缘提取响 76.3 用Roberts交叉算子进行边缘提取 86.4 用Sobel算子进行边缘提取 96.5 用Laplace算子进行边缘提取 106.6 用Kirsch算子进行边缘提取 1

2、16.7 用Marr算子进行边缘提取 127、结论 138、讨论 148.1 利用中值滤波和高斯平滑去除图像噪声的比较 148.2 不同的模板大小对中值滤波的影响以及尺度(方差)对图像平滑的影响 148.3 利用不同算子对造影图像进行边缘提取的比较 148.4 阈值选取对边缘检测的影响 159、主要参考文献 151、开题背景边缘检测技术在图像处理与计算机视觉中占有特殊的位置，因为边缘是所要提取目标和背景的分界线，提取出边缘才能将目标和背景区分开来。在图像中，边界表明一个特征区域的终结和另一个特征区域的开始，边界所分开区域的内部特征或属性是一致的，而不同区域内部的特征或属性是不同的，边缘检测技术

3、正是利用物体和背景在某种图像特征上的差异来实现的，这些差异包括灰度、颜色或纹理特征1。边缘检测技术实际上就是检测图像特征发生变化的位置，它是底层视觉处理中最重要的环节之一，也是图像分割、目标区域识别、区域形状提取等图像分析方法的基础2。近年来，边缘检测技术已经逐渐被广泛地应用于医学、林业、遥感、工程技术中的零部件检查、故障诊断以及探伤等更多领域。但是在实际应用中，不断涌现的先进新技术、新方法给人们如何选择合适的、更佳的图像边缘检测方法提出了新的难题3。2、课题目的在学习了医学成像系统和医学图像处理等课程的基础上，以冠脉造影图像为例，利用Matlab软件，学习医学图像的图像增强和边缘提取等。3、

4、课题主要内容 1.分别利用中值滤波和高斯平滑去除图像中的椒盐噪声和随机噪声，分析不同的模板大小对中值滤波的影响，以及尺度(方差)对图像平滑的影响；2.对造影图像分别用灰度梯度算子、Roberts交叉算子、Sobel算子、Kirsch算子、Laplace算子、Marr算子进行边缘提取，并对结果进行讨论，并研究阈值选取对边缘检测的影响。4、原理和方法4.1 边缘检测技术概述边缘可定义为在局部区域内图像的差别，他表现为图像上的不连续性，包括灰度级的突变，纹理结构的突变，颜色的变化等，也即边缘是不同区域的分界线，指图像局部灰度显著变化的部分 。边缘提取首先要检出图像局部特征的不连续性，然后再将这些不连

5、续的边缘像素连成完备的边界。边缘的特征是沿边缘走向的像素变化平缓而垂直方向的像素变化剧烈。所以，提取边缘的算法就是检出符合边缘特性的边缘像素的数学算子。由于常常无法实现确定轮廓取向，因而挑选用于轮廓增强的微分算子时，必须选择那些不具备空间方向性的和具有旋转不变的线性微分算子 。在经典的边缘检测方法中，通常使用灰度梯度算子、Roberts交叉算子、Sobel算子、Kirsch算子、Laplace算子、Marr算子等。4.2 灰度梯度算子当图像微分运算后，每个像素处有一梯度值反应该点边缘强度。灰度梯度由幅值和方向两个因素构成。灰度梯度算子检测图像像素灰度在水平方向或垂直方向上的变化。因此，这种算子

6、对水平方向或垂直方向上的边缘较为敏感，仅计算相邻像素的灰度差，对噪声比较敏感，无法抑制噪声的影响。4.3 Roberts交叉算子 Roberts交叉边缘检测算子根据任意一对垂直方向上的差分可用来计算梯度的原理，采用对角线方向相邻两像素之差，实际时该点连续梯度的近似值，而不是所预期的点出的近似值。另外该算子提取的边缘较粗，对噪声很敏感，并且需要手动地指定一个阈值，对不同的图像使用不同的阈值后，边界提取效果差别很大，因此Roberts交叉算子一般使用自适应的阈值，并对提取的边缘进行细化处理，这样将能有效地改善Roberts交叉算子的边界提取效果。一般适用于低噪声图像6。4.4 Sobel算子Sob

7、el算子是边缘检测器中最常用的算子之一，一般可以利用快速卷积函数，简单有效，因此应用很广泛。该算子有两个，一个是检测水平边缘的；另一个是检测垂直边缘的。水平与垂直相比，Sobel算子对于像素的位置的影响做了加权，因此效果较好。Sobel算子另一种形式是各向同性Sobel算子，也有两个，一个是检测水平边沿的，另一个是检测垂直边沿的。各向同性Sobel算子和普通Sobel算子相比，它的位置加权系数更为准确，在检测不同方向的边沿时梯度的幅度一致。但是该算子并没有将图像的主体和背景严格区分开来，没有严格地模拟人的视觉生理特性，所以提取的图像轮廓有时并不能令人满意。具有采用加权滤波，边缘较宽，剪短点较多

8、的特点，适用于灰度渐变、低噪声图像6。4.5 Kirsch算子对像素考察其8个相邻点。用3个相邻点的甲醛和减去剩下5个相邻点的加权和，并令3个相邻点环绕上述像素移位，多次计算，取其中差值最大者。Kirsch算子从八个方向对图像边缘信息进行提取，能获得较好的图像边缘。4.6 Laplace算子图像灰度值的强烈变化会导致其一阶导数在这个位置处出现一个尖峰脉冲，而在其他位置一阶导数都为0.对于二阶倒数来说，一阶导数形成的尖峰脉冲必然导致其二阶导数存在过零点，而过零点的位置恰好对应原图像中边缘的位置。这表明可以利用二阶导数的过零点检测图像中边缘的存在。Laplace算子正是基于以上原理的一种二阶导数算

9、子，在实际计算中一般采用Laplace算子模板与图像卷积来实现。74.7 Marr算子Marr算子是在Laplace算子的基础上实现的。Laplace算子对噪声比较敏感，为了减少噪声的影响，可以对图像先进性平滑然后运用Laplace算子。根据人眼视觉成像原理，图像中任一个像素点所对应的场景点周围的点对该点光强的贡献成正态分布。因此，平滑函数选择具有正态分布形式的Gauss函数。已经证明，Marr算子是按零交叉检测阶跃边缘的最佳算子。麻烦的是，Marr算子所用卷积模板一般较大，典型的是半径为8-32个像素。 Marr算子具有明显的优点：1.该滤波器中的高斯函数部分可以对图像进行平滑，消除图像中尺

10、度变化小于滤波函数参数的噪声或不必要的细节。而高斯函数在空域与频域具有相同的形式与性质，都是平滑的、定域的，基本上不会引入在原始图像中未出现的变化。2.采用拉普拉斯算子不仅减少了计算量，而且保证了各向同性。3.滤波参数可调，能够在任何需要的尺度上工作。大的滤波器可以用来检测图像的模糊边缘，小的滤波器可以用来检测聚焦良好的图像细节。在该算子中，滤波参数的选择很重要，参数小时边缘位置精度高，但边缘细节变化多；参数大时平滑作用大，但细节损失大，边缘点定位精度低。应根据噪声水平和边缘点定位精度要求适当选取参数。5、步骤 1、阅读并学习课题相关的文献及资料，学习掌握有关图像处理和边缘检测技术的基本原理和

11、应用。 2、复习并掌握不同边缘检测算子的计算公式，以及它们的检测特点。 3、分别利用中值滤波和高斯平滑去除图像中的椒盐噪声和随机噪声，分析不同的模板大小对中值滤波的影响，以及尺度(方差)对图像平滑的影响。 4、对造影图像分别用灰度梯度算子、Roberts交叉算子、Sobel算子、Kirsch算子、Laplace算子、Marr算子进行边缘提取，并研究阈值选取对边缘检测的影响。 5、分析结果，实验总结。6、结果6.1分别利用中值滤波和高斯平滑去除图像中椒盐噪声和随机噪声 从下图1和图2可以看出中值滤波的除噪效果比高斯平滑效果更好，但是前者处理后的图像对比度不如后者。中值滤波时，随着模板尺寸增大时，

12、对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少。且在滤除椒盐噪声时最为有效。对于高斯滤波，方差越大，对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少，且对随机噪声的滤除较差。 图1 不同模板中值滤波去除图像中椒盐噪声图2 不同尺度高斯平滑去除图像中随机噪声6.2用灰度梯度算子进行边缘提取响从图3可以看出：随着阈值的减小，图像边缘提取的效果有所增强，但是噪声影响也变大。因仅计算相邻像素的灰度差，对噪声比较敏感，无法抑制噪声的影响。图3 不同的阈值灰度梯度算子下进行边缘提取 6.3 用Roberts交叉算子进行边缘提取从图4可以看出：随着阈

13、值的减小，边缘提取的效果增强，即边缘越清晰，但是噪声的影响增加，边缘较粗，但效果较梯度算子好，对噪声敏感。图4 不同的阈值Roberts交叉算子下进行边缘提取6.4 用Sobel算子进行边缘提取从下图5可以看出：随着阈值的减小，边缘提取的效果有所增强，但是噪声的影响变大。能进一步抑制噪声，但边缘定位精度不高，检测的边缘较宽，出现多像素宽度。其对这种渐变的图像边缘提取较好。图5 不同的阈值Sobel算子下进行边缘提取6.5 用Laplace算子进行边缘提取从下图6可以看出：Laplace算子是二阶导数算子，因此对噪声非常敏感，随着阈值减小，边缘提取效果虽稍有增强，但是噪声的影响太大抑制无法分辨边

14、缘。图6 不同的阈值Laplace算子下进行边缘提取6.6 用Kirsch算子进行边缘提取从下图7可以看出：随着阈值的减小，边缘提取的效果增强，但噪声的影响变大，但是其对噪声的抑制较强，较好地平滑了噪声信息，但边缘较粗。图7 不同的阈值Kirsch算子下进行边缘提取6.7 用Marr算子进行边缘提取从下图8可以看出：虽然该算子克服了Laplace算子抗噪声比较差的缺点，但是在抗噪声的同时，也将原有的比较尖锐的边缘也平滑掉了，造成部分边缘无法检测到。图8 不同的阈值Marr算子下进行边缘提取7、结论1 中值滤波的除噪效果比高斯平滑效果更好，但是前者处理后的图像对比度不如后者。中值滤波时，随着模板

15、尺寸增大时，对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少。且在滤除椒盐噪声时最为有效。对于高斯滤波，方差越大，对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少，且对随机噪声的滤除较差。2 灰度梯度算子仅计算相邻像素的灰度差，对噪声比较敏感，无法抑制噪声的影响。3Robert算子提取的边缘较粗，对噪声很敏感，并且需要手动地指定一个阈值，对不同的图像使用不同的阈值后，边界提取效果差别很大，因此Roberts交叉算子一般使用自适应的阈值，并对提取的边缘进行细化处理，这样将能有效地改善Roberts交叉算子的边界提取效果。4Sobel算子利用

16、垂直和水平方向模板可以得到不同方向的边缘，由于微分增强了噪音，平滑效果是Sobel算子特别引人注意的特征。5Kirsch算子的边缘化效果较好，能够较全面的提取边缘信息。6Laplace算子是二阶导数算子，因此对图像中的噪声非常敏感。另外，它常产生双像素宽的边缘，而且不能提供图像边缘的方向信息。所以很少直接用作边缘检测。7Marr算子中的高斯函数部分可以对图像进行平滑，消除图像尺度变化小于滤波参数的噪声或不必要的细节，采用拉普拉斯算子不仅减少了计算量而且保证了各向同性。8、讨论8.1 利用中值滤波和高斯平滑去除图像噪声的比较 中值滤波的除噪效果比高斯平滑效果更好，但是前者处理后的图像对比度不如后

17、者。8.2 不同的模板大小对中值滤波的影响以及尺度(方差)对图像平滑的影响中值滤波时，随着模板尺寸增大时，对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少。且在滤除椒盐噪声时最为有效。对于高斯滤波，方差越大，对噪声的消除效果有所增强，不过同时所得到图像变得更为模糊，可视的细节逐步减少，且对随机噪声的滤除较差。8.3 利用不同算子对造影图像进行边缘提取的比较1 灰度梯度算子仅计算相邻像素的灰度差，对噪声比较敏感，无法抑制噪声的影响。2Robert算子提取的边缘较粗，对噪声很敏感，并且需要手动地指定一个阈值，对不同的图像使用不同的阈值后，边界提取效果差别很大，因此Rob

18、erts交叉算子一般使用自适应的阈值，并对提取的边缘进行细化处理，这样将能有效地改善Roberts交叉算子的边界提取效果。3Sobel算子利用垂直和水平方向模板可以得到不同方向的边缘，由于微分增强了噪音，平滑效果是Sobel算子特别引人注意的特征。4Kirsch算子的边缘化效果较好，能够较全面的提取边缘信息。5Laplace算子是二阶导数算子，因此对图像中的噪声非常敏感。另外，它常产生双像素宽的边缘，而且不能提供图像边缘的方向信息。所以很少直接用作边缘检测。6Marr算子中的高斯函数部分可以对图像进行平滑，消除图像尺度变化小于滤波参数的噪声或不必要的细节，采用拉普拉斯算子不仅减少了计算量而且保证了各向同性。8.4 阈值选取对边缘检测的影响关于阈值的选择，对于不同的图像，不同的算子，有不同的最佳阈值。我们要根据不同的图像及算子进行选择，以获得最好的边缘提取。但是随着阈值的增大，提取的边缘越不清晰，这是一定的。9、主要参考文献