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

《遥感图像处理实验报告图像的空间域滤波处理》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、遥感图像处理实验报告（2013—2014学年第1学期）实验名称：图像的空间域滤波处理实验时间：实验地点：指导教师：专业班级：姓名：学号：一：实验目的掌握Matlab编程语言的基本语法掌握Matlab中imread()、imfilter()、imshow()、mat2gray()、等相关函数的用法及意义通过Matlab程序语言完成图像的空间域滤波处理，以达成边缘检测的目的二：实验内容在Matlab中，分别运用经典一阶梯度、Robert算子、Sobel算子、Prewitt算子及Laplacian算子进

2、行目标图像的边缘检测，互相对比实验结果，分析差异三：实验平台Windows7UltimateSP1Matlab2012b四：实验代码及Matlab使用心得（注释中）%清空工作空间clc;%清屏clearall;%关闭所有绘图窗口closeall;%定义一个矩阵I用来存储待处理图像,用imread()函数读取的数据格式为unit8I=imread('C:UsersAdministratorDesktopMatlabdatabarbara.png');%创建一个绘图窗口，并显示原图像fig

3、ure(1);imshow(I);title('原图');%把unit8格式数据转换为doulbe格式，以便滤波计算I=double(I);%%%%%%%%%%%%%%%%%%%%%%%%%%%%经典一阶梯度%%%%%%%%%%%%%%%yjtdx=[000;-110;000];yjtdy=[0-10;010;000];%imfilter()函数为一个滤波函数%其中有两个参数，外加一系列附加指令%如下中的I代表被处理图像%yjtdx代表滤波器%'replicate'是一个附加指令，因为处理过的数据

4、，可能会超出原图像数据大小的上下限%这个附加指令令超出上下限的数据为最接近上下限的数据ImageYjtdx=imfilter(I,yjtdx,'replicate');ImageYjtdy=imfilter(I,yjtdy,'replicate');%将x与y方向上的处理结果进行几何平均，作为处理结果ImageYjtd=sqrt(ImageYjtdx.^2+ImageYjtdy.^2);%figure()函数创建一个绘图窗口figure(2);%绘制处理后的图像%在用imshow()函数显示图像时

5、要注意，因为现在的处理结果是double类型的，而不是unit8类型%而imshow()在显示图像时，只有unit8类型的矩阵图像数据可以直接输入直接显示%而对于double类型的矩阵，imshow()认为其灰度空间是从0到1，而不是从0-256或其它%也就是说，对于double类型的矩阵，imshow()函数不知道矩阵所代表的图像的灰度级数是多少%这时就需要用mat2gray()函数对double型矩阵做归一化，所谓归一化，就是把原矩阵中的数据%映射至0-1的数据空间中，这时，经mat2gray

6、()函数处理过的矩阵，其中所有元素都在0-1范围内%再用imshow()函数进行显示，才可得到正常的显示结果imshow(mat2gray(ImageYjtd));title('经典一阶梯度处理结果');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Robert算子%%%%%%%%%%%%%%%%%Robert1=[-100;010000];Robert2=[0-10;100000];

7、ImageRobert1=imfilter(I,Robert1,'replicate');ImageRobert2=imfilter(I,Robert2,'replicate');ImageRobert=sqrt(ImageRobert1.^2+ImageRobert2.^2);figure(3);imshow(mat2gray(ImageRobert));title('Robert算子处理结果');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

8、%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Sobel算子%%%%%%%%%%%%%%%%%%Sobel1=[121;000;-1-2-1];Sobel2=[-101;-202;-101];ImageSobel1=imfilter(I,Sobel1,'replicate');ImageSobel2=imfilter(I,Sobel2,'replicate');ImageSobel=sqrt(ImageSobel1.^2+ImageSobel2.^2);figure