matlab学习_matlab

clear all; clc; M=6; N=9; %solarcell的行数与列数 I=imread('.\testpic\milkdrop.gif'); %I=imread('.\testpic\peppers.gif'); %I=imread('.\testpic\left.png'); %I = imread('.\testpic\bacteria.BMP'); %I = imread('.\use\IMAGE15.bmp'); %I = imread('.\IRPIC\IR_all_s_g.BMP'); [m,n]=size(I); imshow(I); %中值滤波 R=zeros(m,n); for i=2:m-1 for j=2:n-1 aa=[I(i-1,j-1) I(i-1,j) I(i-1,j+1) I(i,j-1) I(i,j) I(i,j+1) I(i+1,j-1) I(i+1,j) I(i+1,j+1)]; ii=0; for jj=1:9-ii for ii=1:8 if(aa(ii)>aa(ii+1)) cc=aa(ii); aa(ii)=aa(ii+1); aa(ii+1)=cc; end end end R(i,j)=aa(5); end end S=uint8(R); figure; imshow(S); %otsu 自动阈值计算 Y=zeros(1,256); thresholdValue=1; %阈值初值 for i=2:m-1 for j=2:n-1 k=S(i,j); Y(1,k+1)=Y(1,k+1)+1; end end sum=0; csum=0; n0=0; for k=1:256 sum=sum+k*double(Y(1,k)); n0=n0+Y(k); end fmax=-1; n1=0; for k=1:256 n1=n1+Y(k); if(n1==0) continue; end n2=n0-n1; if(n2==0) break; end csum=csum+k*double(Y(k)); m1=csum/n1; m2=(sum-csum)/n2; sb=n1*n2*(m1-m2)^2; if(sb>fmax) fmax=sb; thresholdValue=k; end end %图像二值化 U=zeros(m,n); for i=1:m for j=1:n if(S(i,j)>=thresholdValue) U(i,j)=0; else U(i,j)=1; end end end figure; imshow(U); %采用腐蚀的方法进行边缘检测 B=[1 1 1;1 1 1;1 1 1]; qsum=0; QQ=zeros(m,n); for i=2:m-1 for j=2:n-1 for ii=1:3 for jj=1:3 qsum=qsum+~xor(B(ii,jj),(1-U(i+ii-2,j+jj-2))); end end if(qsum==9) QQ(i,j)=0; else QQ(i,j)=1; end qsum=0; end end SS=abs(QQ-U); figure; imshow(SS); title('腐蚀处理后得到的边界图像'); SS1=SS; for i=1:m for j=1:n if(i<=5 | j<=5 | i>=m-5 | j>=n-5) SS1(i,j)=0; end end end figure; imshow(SS1); title('纯背景边界图像'); C=[0 1 0;1 1 1;0 1 0]; X=zeros(m,n); flag_x=0; for i=1:m if(flag_x==1) break; end for j=1:n if(SS1(i,j)==1) X(i,j)=1; flag_x=1; break; end end end figure; imshow(X); qsum=0; k=1; QQ=zeros(m,n); X2=zeros(m,n); ysum=1; while(ysum~=0) if(k>1) X=X2; end k=k+1; for i=2:m-1 for j=2:n-1 for ii=1:3 for jj=1:3 qsum=qsum+C(ii,jj)*X(i+ii-2,j+jj-2); end end if(qsum>=1) QQ(i,j)=1; else QQ(i,j)=0; end qsum=0; end end X2=QQ&SS1; Y2=(X2~=X); ysum=0; for iii=1:m for jjj=1:n ysum=ysum+Y2(iii,jjj); end end end figure; imshow(X2); figure; imshow(X) title('膨胀处理后的图像');