遥感图像处理实验报告_图像的空间域滤波处理.docxVIP

遥感图像处理 实验报告 （2013 —2014 学年第1学期） 实验名称： 图像的空间域滤波处理 实验时间： 实验地点： 指导教师： 专业班级： 姓 名： 学 号：  一：实验目的 掌握Matlab编程语言的基本语法 掌握Matlab中imread()、imfilter()、imshow()、mat2gray()、等相关函数的用法及意义 通过Matlab程序语言完成图像的空间域滤波处理，以达成边缘检测的目的 二：实验内容 在Matlab中，分别运用经典一阶梯度、Robert算子、Sobel算子、Prewitt算子及Laplacian算子进行目标图像的边缘检测，互相对比实验结果，分析差异 三：实验平台 Windows 7 Ultimate SP1 Matlab 2012b 四：实验代码及Matlab使用心得（注释中） %清空工作空间 clc; %清屏 clear all; %关闭所有绘图窗口 close all; %定义一个矩阵 I 用来存储待处理图像,用imread()函数读取的数据格式为unit8 I = imread(C:\Users\Administrator\Desktop\Matlab\data\barbara.png); %创建一个绘图窗口，并显示原图像 figure(1); imshow(I); title(原图); %把unit8格式数据转换为doulbe格式，以便滤波计算 I = double(I); %%%%%%%%%%%%%%%%%%%%%%%%%%%%经典一阶梯度%%%%%%%%%%%%%%% yjtdx = [ 0 0 0; -1 1 0; 0 0 0]; yjtdy = [ 0 -1 0; 0 1 0; 0 0 0]; %imfilter()函数为一个滤波函数 %其中有两个参数，外加一系列附加指令 %如下中的 I 代表被处理图像 %yjtdx 代表滤波器 %replicate是一个附加指令，因为处理过的数据，可能会超出原图像数据大小的上下限 %这个附加指令 令超出上下限的数据为最接近上下限的数据 ImageYjtdx = imfilter(I, yjtdx, replicate); ImageYjtdy = imfilter(I, yjtdy, replicate); %将x与y方向上的处理结果进行几何平均，作为处理结果 ImageYjtd= sqrt(ImageYjtdx.^2 + ImageYjtdy.^2); %figure()函数创建一个绘图窗口 figure(2); %绘制处理后的图像 %在用imshow()函数显示图像时要注意，因为现在的处理结果是double类型的，而不是unit8类型 %而imshow()在显示图像时，只有unit8类型的矩阵图像数据可以直接输入直接显示 %而对于double类型的矩阵，imshow()认为其灰度空间是从0到1，而不是从0-256或其它 %也就是说，对于double类型的矩阵，imshow()函数不知道矩阵所代表的图像的灰度级数是多少 %这时就需要用mat2gray()函数对double型矩阵做归一化，所谓归一化，就是把原矩阵中的数据 %映射至0-1的数据空间中，这时，经mat2gray()函数处理过的矩阵，其中所有元素都在0-1范围内 %再用imshow()函数进行显示，才可得到正常的显示结果 imshow(mat2gray(ImageYjtd)); title(经典一阶梯度处理结果); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%Robert算子%%%%%%%%%%%%%%%%% Robert1 =[ -1 0 0; 0 1 0 0 0 0]; Robert2 =[ 0 -1 0; 1 0 0 0 0 0]; ImageRobert1 = imfilter(I, Robert1, replicate); ImageRobert2 = imfilter(I, Robert2, replicate); ImageRobert =