1、第一次试验pic = imread(1.jpg); data=rgb2gray(pic);R = data(:,:,1);G = data(:,:,2);B = data(:,:,3);figure,imshow(pic);figure,imshow(im);figure,imshow(R);figure,imshow(G);figure,imshow(B); figure,imhist(im);figure,imhist(im1);figure,imhist(im2); figure,imhist(im3);num=ones(1,256);for i=0:255temp=find(data=
2、i);num(i+1)=length(temp);endn=0:255;subplot(211)stem(n,num);axis(0 255 0 max(num)+100);subplot(212)imhist(data);title(自变程序与 imhist 函数产生灰度直方图比较 );00.20.40.60.811.21.41.61.82x 1040 50 100 150 200 2500200040006000800010000120001400016000180000 50 100 150 200 25000.20.40.60.811.21.41.61.82x 1040 50 100
3、150 200 2500 50 10 150 20 25002004006008001001200140000.511.52x 104 与与与与与imhist与与与与与与与与与与与0 50 10 150 20 250第二次试验pic=imread(1.jpg);data=rgb2gray(pic);figure,imshow(data);title(原始灰度图像,fontsize,20);m,n=size(data);data0=zeros(m,n);mob=1 1 1;1 1 1;1 1 1/9;for i=2:m-1for j=2:n-1part=data(i-1:i+1,j-1:j+1)
4、;part=double(part); temp=mob.*part;data0(i,j)=sum(sum(temp);endendfigure,imshow(uint8(data0);title(遍历所有的像素点对图像进行平滑,fontsize,20);w = fspecial(average,3 3);g= imfilter(data,w,conv,replicate,same);figure,imshow(g);title(使用 matlab 自带的函数对图像进行平滑处理,fontsize,20);与与与与与与与与与与与与与与与与与与与与与与与matlab与与与与与与与与与与与与与与第三
5、次试验pic=imread(1.jpg);data=rgb2gray(pic); %转化为灰度图figure,imshow(data);title(原始灰度图);m,n=size(data); %返回图像的大小data1=ones(m,n);data2=ones(m,n);data3=ones(m,n);data4=ones(m,n);B1=fspecial(sobel);B2=mhori;for i=2:m-1for j=2:n-1 %遍历元素part=data(i-1:i+1,j-1:j+1); %取图像的一部分part=double(part);temp1=part.*B1; temp2
6、=part.*B2; %图像的部分与模板对应元素相乘data1(i,j)=sum(sum(temp1);data2(i,j)=sum(sum(temp2); %求相乘的和endendfigure,imshow(uint8(data1);title(图像水平梯度图);figure,imshow(uint8(data2);title(图像垂直梯度图); data3=data1+data2; %水平与垂直叠加 figure,imshow(uint8(data3);title(图像完整的梯度图);data4=data3+double(data);figure,imshow(uint8(data4);t
7、itle(原始图像加图像梯度);w = fspecial(sobel);g= imfilter(data,w,conv,replicate,same);figure,imshow(g);与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与与第四次试验pic = imread(4.jpg); data=rgb2gray(pic);figure,imshow(data);title(原始灰度图像,fontsize,15);r,c=size(data);f1 =zeros(r,c);f2 =zeros(r,c);w=fspecial(sobel);g=imfilter(data,w
8、,conv,replicate,same);index1=find(g=5);index2=find(g=100);f1(index1)=1;f2(index2)=1;figure,imshow(f1);title(sobel 边缘检测阀值为 5,fontsize,15);figure,imshow(f2);title(sobel 边缘检测阀值为 100,fontsize,15);BW0=edge(data,canny);figure,imshow(BW0)title(canny 边缘检);BW1=imnoise(data,salt title(对图像添加椒盐噪音);BW2=edge(BW1,canny,0.3);figure,imshow(BW2);title(canny 边缘检,thresh 值为 0.3);BW3=edge(BW1,canny,0.5);figure,imshow(BW3);title(canny 边缘检，thresh 值为 0.5); sobel与与与与与与与5sobel与与与与与与与100canny与与与与与与与与与与与与canny与与与,thresh与与0.3canny与与与与thresh与与0.5
