matlab.rar_matlab 实例分割_去噪_图像分割_图像去噪_图像处理

clear; close all; clear all; f=imread('saturn.png'); f = rgb2gray(f); %f1=double(f); [MC NC]=size(f); [a1,b1,c1,d1]=dwt2(f,'haar'); [a,b,c,d]=dwt2(a1,'haar'); [m,n]=size(d); [m1 n1]=size(a1); aa=0; ah=0; av=0; for i=1:m for j=1:n aa=aa+a(i,j); ah=ah+abs(b(i,j)); av=av+abs(c(i,j)); end end aa=aa/MC^2; ah=ah/MC^2; av=av/MC^2; for i=1:m for j=1:n if((a(i,j)>1/4*aa)&(a(i,j)<2*aa)) wa(i,j)=1; else wa(i,j)=0; end if((abs(b(i,j))>3/4*ah)&(abs(b(i,j))<4*ah)) wb(i,j)=1; else wb(i,j)=0; end if((abs(c(i,j))>3/4*ah)&(abs(c(i,j))<4*av)) wc(i,j)=1; else wc(i,j)=0; end if(wa(i,j)+wb(i,j)+wc(i,j)==0) w(i,j)=0; else w(i,j)=1; end end end figure,imshow(w) title('提取出的图像')