一个实验报告.docVIP

数字图像处理实验报告 实验名称：利用二阶导数进行边缘增强 班级：信科0801 姓名：董飞 学号：2003080131 指导老师：程浩 完成时间：2012.1.5 1.实验目的 基本要求：了解二阶导数与图像及图像边缘之间关系，掌握利用导数进行边缘增强的方法； 重点：导数与图像特性之间的联系； 难点：利用导数进行边缘增强时的系数选取。 2.实验内容 （一）Laplace算子与图像边界之间的关系 对于给定的数字图像，在像素的梯度为 观察边界附近的Laplace算子与图像其它部分的Laplace算子的区别； （二）边缘增强型滤波 我们利用一下公式对给定图像进行边缘增强： 其中A是选定的一个常数。 3.实验原理 （一）拉普拉斯算子 1.定义 　　拉普拉斯算子是n维欧几里德空间中的一个二阶微分算子，定义为梯度的散度。因此如果f是二阶可微的实函数，则f的拉普拉斯算子定义为： 　　(1) f的拉普拉斯算子也是笛卡儿坐标系xi中的所有非混合二阶偏导数： 　　(2) 作为一个二阶微分算子，拉普拉斯算子把C函数映射到C函数，对于k ≥ 2。表达式(1)（或(2)）定义了一个算子Δ : C(R) → C(R)，或更一般地，定义了一个算子Δ : C(Ω) → C(Ω)，对于任何开集Ω。 　　函数的拉普拉斯算子也是该函数的黑塞矩阵的迹 2. 推广 　　拉普拉斯算子可以用一定的方法推广到非欧几里德空间，这时它就有可能是椭圆型算子，双曲型算子，或超双曲型算子。 　　在闵可夫斯基空间中，拉普拉斯算子变为达朗贝尔算子： 达朗贝尔算子通常用来表达克莱因-高登方程以及四维波动方程。第四个项前面的符号是负号，而在欧几里德空间中则是正号。因子c是需要的，这是因为时间和空间通常用不同的单位来衡量；如果x方向用寸来衡量，y方向用厘米来衡量，也需要一个类似的因子。 拉普拉斯-贝尔特拉米算子 　　主条目：拉普拉斯-贝尔特拉米算子　　拉普拉斯算子也可以推广为定义在黎曼流形上的椭圆型算子，称为拉普拉斯-贝尔特拉米算子。达朗贝尔算子则推广为伪黎曼流形上的双曲型算子。拉普拉斯-贝尔特拉米算子还可以推广为运行于张量场上的算子（也称为拉普拉斯-贝尔特拉米算子）。　　另外一种把拉普拉斯算子推广到伪黎曼流形的方法，是通过拉普拉斯-德拉姆算子，它运行于微分形式。这便可以通过Weitzenbouml;ck恒等式来与拉普拉斯-贝尔特拉米算子联系起来。 3.拉普拉斯边缘检测 拉普拉斯算子是对二维函数进行运算的二阶导数标量算子。 由于拉普拉斯算子是一个二阶导致，它将在边缘处产生一个陡峭的零交叉。拉普拉斯算子是一个线性的、移不变算子，它的传递函数在频域空间的原点为零。因此，一个经拉普拉斯滤波过的图像具有零平均灰度。 function PicOut=Lap_edge(PicInput,thresh)% 本程序能够将BMP格式的黑白灰度图像用拉普拉斯算子进行边缘检测% 格式为 a=Lap_edge(PicInput,thresh) 或者嵌套为 Lap_edge(imread(rice.tif),15)% thresh参数可自选，对于rice.tif这张图来说最合适的值大约为14到18% 使用例子：PicInput=imread(rice.tif);% a=Lap_edge(PicInput,15);%－－－－－－－－－－－－－－－－ BEGIN CODE －－－－－－－－－－－－－－－－ % 作者;archiless lorder% 一，原图像预处理，读入黑白图片并确定长和宽[m,n]=size(PicInput); %确定图片的长和宽 % 二，拉普拉斯变换预处理，定义镜框矩阵和输出矩阵r=m+2; %把图片的长和宽各加2c=n+2;PicFrame=zeros(r,c); %定义二维数组PicFrame，长、宽比Input各多2，成为镜框的尺寸b=zeros(m,n); %定义滤波后的数组 % 三，拉普拉斯运算的三个矩阵Temp=zeros(3); %定义三阶方阵Temp，为临时矩阵op=[0 -1 0;-1 4 -1;0 -1 0]; %定义拉普拉斯算子Result=zeros(3); %定义三阶方阵Result，为运算结果矩阵 % 四，原图像矩阵处理，做一个像框PicFrame(2:m+1,2:n+1)=PicInput; %把原图的矩阵放到新的矩阵PicFrame中心，它的第一行、最后一行、第一列、 %最后一列都是0，即原图矩阵周围有一圈0的边缘，好像给图像加一个像框PicFrame(1,:)=PicFrame(2,:); %把第二行的值赋给第一行PicFrame(r,:)=PicFrame(r-1,:); %把倒数第二行的