vcai93.blog.ir


%-- 5/25/15  2:04 PM --%
a=imread('tt.jpg');
imshow(a)
h=imhist(a(:,:,1));
plot(h)
b=a(:,:,1);
c=zeros(size(b));
c(b<100)=1;
figure
imshow(c)
p=h/(size(b,1)*size(b,2));
figure
plot(p)
clc
mean(b(:))
x=0:255;
size(x)
size(p)
x*p
p(p==0)=1;
p'*log2(p)
clc
b=a(:,:,2);
plot(imhist(b))
h=imhist(b);
p=h/(size(b,1)*size(b,2));
close all
plot(p)
T=zeros(size(p);
T=zeros(size(p));
for i=1:256
T(i)=sum(p(1:i));
end
figure
plot(T)
Tt=round(T*255);
n=b;
size(b)
for i=1:256
n(b==i)=Tt;
end
for i=1:256
n(b==i)=Tt(i);
end
hn=imhist(n/255);
figure
plot(hn)
imshow(n/255)
imshow(n)
hn=imhist(n);
plot(hn)
figure
plot(hn)
clc
f=zeros(1,256);
f=0:255;
f
clc
for i=1:256
g(i)=1/(2*pi*100)^.5*exp(-(i-163)^2/(2*100));
end
plot(g)
figure
plot(imhist(b))
figure
plot(Tt)
for i=1:256
tg(i)=sum(g(1:i));
end
figure
plot(tg)
figure
plot(Tt/(size(b,1)*size(b,2)))
clc
close all
clear
a=imread('tt.jpg');
b=a(:,:,2);
b=double(b);
size(b)
c=-1*b(1:767,:)+2*b(2:768,:);
imshow(b/255);
figure
imshow(c,[]);
c=-1*b(1:767,2:)+2*b(2:768,:);
clc
d=-1*b(1:766,2:1023)-1*b(2:767,1:1022)+4*b(2:767,2:1023)-1*b(2:767,3:1024)-1*b(3:768,2:1023);
figure
imshow(d,[]);
clc
clear
a=imread('tt.jpg');
b=a(:,:,2);
fb=fft2(b);
fb(1:5,1:5)
imshow(b)
figure
imshow(ft,[])
imshow(fb,[])
imshow(log10(abs(fb)+1),[])
imshow(log10(abs(fftshift(fb))+1),[])
bbf=fftshift(fb);
H=zeros(size(fb));
H(334:434,462:562);
figure
imshow(H)
H(334:434,462:562)=1;
imshow(H)
FH=H.*bbf;
figure
imshow(log10(abs((FH))+1),[])
ifb=fftshift(FH);
fg=ifft2(ifb);
figure
imshow((abs((fg))),[])
H=ones(size(fb));
H(334:434,462:562)=0;
FH=H.*bbf;
ifb=fftshift(FH);
fg=ifft2(ifb);
figure
imshow((abs((fg))),[])
fg(1:5,1:5)