机器视觉实验报告..docx

------大学实 验 报 告作者所在学院： 作者所在专业： 作者所在班级： 作 者 姓 名 ： 实验报告1.1.1实验内容Averaging kernel(3*3 and 5*5)1.1.2算法实现clear all;clc;I=imread(..\PA1\House1.jpg);n=input(请输入滤波器模板大小\n);K1=avefilt(I,n);figure;subplot(1,2,1);imshow(I);title(原始图像);subplot(1,2,2);imshow(K1);title(均值滤波后图像);function d=avefilt(x,n)a=ones(n,n);p=size(x);x1=double(x);x2=x1;for i=1:p(1)-n+1 for j=1:p(2)-n+1 c=x1(i:i+(n-1),j:j+(n-1)).*a; s=sum(sum(c)); x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n); endendd=uint8(x2);1.1.3实验结果1.2.1实验内容Gaussian Kernel (σ =1,2,3 ) Use (2σ +1)×(2σ +1) as size of Kernel (You may write a separate function to generate Gaussian Kernels for different values of σ.)1.2.2算法实现clear all;clc;I=imread(..\PA1\House1.jpg);sigma=input(请输入方差的大小\n);n=2*sigma+1;K1=gaussfilt(sigma,n,I);figure;subplot(1,2,1);imshow(I);title(原始图像);subplot(1,2,2);imshow(K1);title(高斯滤波后图像);function d=gaussfilt(k,n,s) Img = double(s); n1=floor((n+1)/2);%计算图象中心 for i=1:n for j=1:n b(i,j) =exp(-((i-n1)^2+(j-n1)^2)/(4*k))/(4*pi*k); end end Img1=conv2(Img,b,same); d=uint8(Img1);1.2.3实验结果1.3.1实验内容Sobel Edge Operators and Prewitt Edge Operators1.3.2算法实现clcclear allI=imread(..\PA1\House1.jpg);J=imread(..\PA1\House2.jpg);K1= sobeledge(I);%调用sobel水平算子函数K2= sobeledgev(I);%调用sobel垂直算子函数K3= Prewittedge(I);%调用prewitt水平算子函数K4= Prewittedgev(I);%调用prewitt垂直算子函数K5= sobeledge(J);K6= sobeledgev(J);K7= Prewittedge(J); K8= Prewittedgev(J);figure;subplot(2,2,1);imshow(K1);title(sobel算子水平方向边缘检测);subplot(2,2,2);imshow(K2);title(sobel算子垂直方向边缘检测);subplot(2,2,3);imshow(K3);title(prewitt算子水平方向边缘检测);subplot(2,2,4);imshow(K4);title(prewitt算子垂直方向边缘检测);figure;subplot(2,2,1);imshow(K5);title(sobel算子水平方向边缘检测);subplot(2,2,2);imshow(K6);title(sobel算子垂直方向边缘检测);subplot(2,2,3);imshow(K7);title(prewitt算子水平方向边缘检测);subplot(2,2,4);imshow(K8);title(prewitt算子垂直方向边缘检测);sobel算子水平方向函数function d=sobleedge(I) I1=double(I);[m,n]=size(I1);u=zeros(1,9);for i=2:m-1 for j=2:n-1 u(1)=-1*I(i-1,j-1);u(2)=-2*I(i-1,j);u(3)=-1*I(i-1,j+1); u(4)=0*I(i,j-1);u(5)